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We review recent developments in our understanding of how impurities influence the electronic 
states in the bulk of superconductors. Our focus is on the quasi-localized states in the vicinity of 
impurity sites in conventional and unconventional superconductors and our goal is to provide a 
unified framework for their description. The non-magnetic impurity resonances in unconventional 
superconductors are directly related to the Yu-Shiba-Rusinov states around magnetic impurities 
in conventional s-wave systems. We review the physics behind these states, including quantum 
phase transition between screened and unscreened impurity, and emphasize recent work on d-wave 
superconductors. The bound states are most spectacularly seen in scanning tunneling spectroscopy 
measurements on high- Tc euprates, which we describe in detail. We also discuss very recent 
progress on the states coupled to impurity sites which have their own dynamics, and impurity 
resonances in the presence of an order competing with superconductivity. Last part of the review 
is devoted to influence of local deviations of the impurity concentration from its average value 
on the density of states in s-wave superconductors. We review how these fluctuations affect the 
density of states and show that s-wave superconductors are, strictly speaking, gapless in the 
presence of an arbitrarily small concentration of magnetic impurities. 
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electronic correlations, and producing electronic states 
that are impossible in the bulk of a clean system. Un- 
til recently this idea has not been emphasized enough in 
the study of correlated electron systems, but in recent 
years more and more efforts are focused on understand- 
Average density of states in superconductors with impuBaitelab Ochanges produced by disorder in a wide variety of 

strongly interacting electronic matter. One of the most 
promising directions is the study of disorder near quan- 
tum critical points, where several types of order compete 
and exist in delicate balance that impurities have the 
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I. INTRODUCTION 

A. Aim and scope of this article 

Real materials are not pure. Sometimes excessive 
impurities hinder observations of beautiful physics that 
would be there in cleaner systems. Magnetic disorder 
destroys the coherence of the superconducting state. At 
the very least, in conventional metals, impurities lead to 
higher resistivity. It is therefore very tempting to treat 
impurities as unfortunate obstacles to our understanding 
of the true underlying physics of the systems we study, 
strive to make cleaner and better materials, and ignore 
imperfections whenever possible. 

Yet sometimes impurities directly lead to the desired 
physical properties. They are crucial in achieving func- 
tionality of doped semiconductors: undoped semicon- 
ductors are just band insulators and not useful for ap- 
plications in electronics. The entire multi-billion dollar 
semiconducting electronics industry is based on the pre- 
cise control and manipulation of electronic states due to 
dopant (impurity) states. 

Consequently, sensitivity of a physical system to dis- 
order can be a blessing in disguise. It can lead not only 
to achieving new applications but also to uncovering the 
nature of exotic ground states, elucidating properties of 



power to tip in favor of one of the orders (Milhs, 2003). 

This is a review of the impurity effects on the elec- 
tronic states in superconductors. The main purpose of 
our article is to give a reader an appreciation of recent 
developments, review the current understanding and out- 
line further questions on impurity effects in conventional, 
and especially unconventional superconductors. Super- 
conductors present probably the first example of a non- 
trivial many-electron system where effects of disorder on 
the electronic states were studied experimentally and the- 
oretically, and this review focuses on these effects. 

The subject of impurity effects in superconductors is 
well established and well covere d, see, for examp l e, ex- 
cellent textbooks and reviews ("Abrikos ov et IdL 119631 
Annetd 11990: Fetter and Walccka, 1 97l|; Ide Gennes , 



1989t ISchrieffeil Il96l ISigrist and Ue^ Il99lt lTink~ . 



1996|) . The main classical results, such as Abrikosov- 
Gor'kov theory of pairbreakin g by magnetic impurities 
([Abrikosov and Gorkovl ll96Cl|) . and Anderson theorem, 
that explains why non-magnetic impu rities do not de- 
stroy conventional superconductivity, l|AndersonL Il959|) 
are well known from the 60s, and are now taught in 
graduate school. The need to review the subject arose 
since a) there are many new results; b) the analyses of 
the classical papers have been substantially modified in 
applications to novel materials; c) the emphasis of the 
study of the impurity effects shifted from macroscopic to 
microscopic length scales. 

From the early days of superconductivity, impurity 
doping was one of the most important tools to identify 
the nature of the pairing state and microscopic proper- 
ties. A classical study of the role of magnetic impuri- 
ties in conventional s uperconductors was carried out by 
IWoolf and E'en 119651 and followed by many detailed in- 
vestigations. Both magnetic and nonmagnetic impuri- 
ties are pairbreakers in unconventional superconductors, 
and often impurity suppression of superconductivity is a 
strong early hint of the unconventional pairing state. For 
example, the rapid suppression of the transition temper- 
ature, Tc, in Al doped SrRu04 superconductor was the 
first an d strong indicatio n that i t is a p-wave supercon- 
ducto r ijMackenzie et all Il998t iMackenzie and Maenol 
I2003D . 

In the past two decades we have witnessed a tremen- 
dous growth of the number of novel superconductors. 
Many of them belong to the general class of strongly 
correlated electron systems, and, as a result of Coulomb 
interaction, the superconductivity is unconventional, see 
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below. Study of the effect of impurities on unconven- 
tional superconductors is a still developing field, yet it is 
mature enough to warrant an overview. 

Sometimes superconducting state emerges from com- 
petition between different phases, such as magnetically 
ordered and paramagnetic in high- temperature cuprates, 
organic materials and heavy fermion systems. Experi- 
mentally, superconductivity often is the strongest when 
the two competing states are nearly degenerate, near 
quantum critical points. This is the case for e xample for 
Ce based heavy fermion materi als ijSidorov et ai, 2002) 
and UGe2 |Saxena ei ali l2000j) . Study of impurity ef- 
fects in these materials allows (at least, in principle) to 
determine the nature of the superconducting state and 
reveal competing electronic correlations. 

This has driven in part the study of impurity effects 
in high-Tc superconductors. At present, despite much 
progress, there is no complete microscopic description 
and certainly no consensus in the community on the 
mechanism of superconductivity. Study of competing 
orders in the neighborhood of impurity atoms has the 
potential to reveal the nature and origin of the supercon- 
ducting state. 

The new states and structures that appear due to dis- 
order often are confined to micro- or mesoscopic length 
scales. They would remain in the realm of academic dis- 
cussion were it not for the development of new techniques 
and probes of disorder. At the time of classical work, ex- 
perimental interest lied solely in macroscopic properties 
of materials: transition temperature, Tc, specific heat, 
and the average density of states (obtained from planar 
junction tunneling measurements) were the experimen- 
tally measured quantities. With perfection of more local 
probes such as nuclear magnetic resonance (NMR), and 
especially with development of scanning tunneling mi- 
croscopy and spectroscopy (STM/STS), it became pos- 
sible to experimentally determine the structures on the 
atomic scales around the impurity sites. Therefore the 
emphasis of theoretical work also shifted to the study of 
these local properties. 

It is therefore timely and useful to review new results 
and ideas about impurity-generated states in supercon- 
ductors. 

We had to be selective about the topics that are in- 
cluded in this article. In the spirit of new approaches, 
our review primarily discusses the physics of the single 
impurity bound or quasi-bound states and the local elec- 
tronic effects in the vicinity of defects. We also discuss 
the physics on the mesoscopic scales, and the behavior 
of impurities in the presence of competing orders. In the 
specific case of high-Tc materials we discuss possible com- 
petition between superconducting state and some com- 
peting orders in the so called pseudogap state of these 
superconductors. 

In all our discussions we restrict ourselves to the study 
of the behavior of the density of states. A more compre- 
hensive review of all the effects that were studied experi- 
mentally and discussed theoretically is a much more diffi- 



cult task and would take substantially more space. We do 
not discuss the behavior of transport coefficients: while 
this is a subject of intense current interest and many im- 
portant results have been obtained there, it is beyond the 
scope of this article. 

To keep this review useful and relevant for people en- 
tering the field, we start with a simple Bardeen-Copper- 
Schrieffer (BCS) model for superconductivity, and use 
a modified version of this model throughout the arti- 
cle. We do not consider the corrections due to strong 
coupling that appear in the Eliashberg analysis; in the 
known cases of electron-phonon interaction th ese correc- 



tions are quantitative rather than qualitative llCarbotte . 
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. In many unconventional 



materials dynamical glue in the self-consistent theory is 
not known. For example, there is an ongoing debate on 
the very nature of the normal state in the high-Tc cuprate 
superconductors. Yet most people agree that the super- 
conducting state of cuprates is less anomalous then the 
normal state, and has the superconducting gap of d-wave 
symmetry. We take a view that at low energies it could 
be described for the purposes of our article by BCS with 
d-wave pairing state. 

At the same time, while this is a review of recent work 
on impurity effects in unconventional superconductors, it 
is emphatically not a comprehensive review of impurity 
effects in high-Tc cuprates. Nature of superconducting 
state, detailed microscopic description and competing or- 
ders in the cuprates are still a subject of intense debate at 
present. There is a number of excellent reviews of physics 
of cupr ates, including s canning tunneling microscopy 
(STM) ijFischer and et al . 2004 ), angle-resolved photoe- 
mission spectroscopv (A RPES) llCa,mnuza,no et al\ . \2004 
iDamascelh et allhOO^ and on nature of pseudogap state 
ijTimusk and Statti . ll999l) . Reader is referred to these re- 
views for the in depth discussion of the issues specific to 
high-Tc materials. 

B. Unconventional superconductivity 

Examples of exotic superconductors discovered in the 
last two decades include high-Tc, heavy fermion super- 
conductors, organic superconductors, SrRu04. The com- 
mon feature of all of them is that they are unconven- 
tional, i.e. the pairing symmetry is not s-wave, in con- 
trast to conventional materials, such as lead. 

Here any superconductor with the gap function that 
transforms according to a trivial representation of the 
point group of the crystal will be called an s-wave super- 
conductor. We call a superconducting order parameter 
unconventional if it transforms as a nontrivial represen- 
tation of the symmetry group. To be more precise, su- 
perconductivity is characterized by an order parameter, 
that describes pairing of fermions with time-reversed mo- 
menta, k and — k, 

4'(k)c,^ = (V'k,a'0-k,/3), (1-1) 
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where a, (3 are spin indices of the paired fermionic states. 
We distinguish between the spin singlet pairing (total 

spin of the pair S = 0), for which '^{\/i)a,i3 = *(k)(icr^'^j, 
where cr(2') is the Pauli matrix in spin space, and spin 
triplet state (5 = 1), when '^ap is a symmetric spinor 
in a, (3. Since the order parameter has to be antisym- 
metric with respect of permutation of fermion operators 
in Eq. the spatial part of ^(k)^^^ is even for spin 

singlet superconductors and odd in the spin-triplet case. 
Expanding in eigenfunctions of orbital momentum, it fol- 
lows that spin singlet pairing corresponds to even orbital 
function of momentum k and hence we call it s- (for 
Z = 0), d-wave (for Z = 2), etc. superconductor in analogy 
with the notation for the atomic states. For spin triplet 
superconductor, the orbital part is an odd function of 
k, and hence spin triplet superconductor can be p-wave 
{I = 1), f-wdve {I = 3) etc. More rigorously one would 
characterize pairing states by the irreducible representa- 
tion of the symmetry o f the cry stal lattice, inclu ding the 
spin-orbit interaction llBloiintl .198 5: Sigrist an d UedaL 
Volovik and Gor'kovl 119841 ). Characterization in 
terms of orbital moment is an oversimplification, and we 
will use this terminology with understanding that the 
correct symmetries are used for a given crystal structure. 
The above classification is given for BCS-like or even fre- 
quency superconductors. This classification is opposite 
for odd-frequency pairing, where, for example, spin sin- 
glet state has odd parity because pairing wave function 
is odd function o f time l)Balatskv and Abrahamsl Il992t 
iBerezinskiil Il974j) . We will focus on BCS like or even- 
frequency superconductors here. 

A reasonable definition of unconventional pairing state, 
that we adopt here, is that the order parameter average 
over the Fermi surface vanishes : 



(1.2) 



Hence superconductors with the constant or nearly con- 
stant order parameter on the Fermi surface are s-wave, 
while p-, d- or higher wave states, where Eq. (|1.2|) 
holds, are signatures of an unconventional superconduc- 
tor. There are several excellent recent reviews that ad- 
dress the unconventional nature of superconducting pair- 
ing states in specific compounds, such as p-wavc super- 
conductivity in SrRu04 (Mack enzie and M acno, 2003,) 
and d- wave state in high-T,. materials l|Annetti . Il990t 
iHarlingenL IT99I iTsuei and Kirtlevll2000(l . 



C. Outline 

We start with the general overview of BCS-like super- 
conductivity. To review the effects of impurities we need 
to discuss the properties of superconductors in general. 
In cuprates, as well as in some heavy fermion systems and 
other novel superconductors, there is some evidence for 
the existence of an order competing with superconductiv- 
ity on all or parts of the Fermi surface. The exact nature 



of the competing order parameter is only conjectured. A 
general feature of all such models is the enhancement of 
the competing order once superconductivity is destroyed, 
for example in the vicinity of a scattering center. It has 
been suggested that the reaction of the system to the 
introduction of impurities can be an important test of 
the order, or even growing correlations towards such an 
order, in the superconducting state. 

The prerequisite for such a test is the detailed un- 
derstanding of the behavior of "simple" superconductors 
with impurities. Work aimed at developing this under- 
standing spans a period of more than 40 years, and some 
of the very recent results continue to be fresh and un- 
expected. Therefore we devote a large fraction of this 
review to the discussion of the properties of supercon- 
ductors with impurities in the absence of any competing 
order. In this case, from a theoretical standpoint, before 
discussing the impurity effects we need to agree upon 
methods to describe the very phenomenon that makes the 
impurity effects so interesting: superconductivity. Even 
in the most exotic compounds investigated so far the su- 
perconducting state itself is not anomalous, in that it 
results from pairing of fermionic quasiparticles, and in 
that these Cooper pairs may be broken by interaction 
with impurities or external fields. 

Impurity effects in conventional superconductors were 
subject of the very ea rly studies by A nderson, so called 
"Anderson the orem". llAndersonl.ll959 ) and by Abrikosov 
and Gor'kov l)Abrikosov and Gorkov , Il960(i . hereafter 
AG. This pioneering work laid the foundation for our 
understanding of impurity effects in conventional and un- 
conventional superconductors, described in terms of elec- 
tron lifetime due to scattering on an ensemble of impuri- 
ties. AG predicted the existence of the gapless supercon- 
ductivity that was subsequently observed in experiments 
(,Woolf and Reif . 1965). The brief summary of the AG 
theory and its extensions to non- s-wave superconductiv- 
ity is given Table HI where effect of impurities on the su- 
perconducting state on average, or globally is listed. 

After intense interest in the early days of the BCS the- 
ory, the subject was considered "closed" in mid-60s, with 
most experimentally relevant problems solved. However, 
as often happens, recently there has been a revival of the 
interest in the studies of "traditional" low-temperature s- 
wave superconductors with magnetic and non-magnetic 
impurities, with many new theoretical and experimental 
results changing our perspective on this classical prob- 
lem. 

A special place in this review is devoted to the study of 
impurity induced local bound states or resonances. This 
is an old subject, going back to the 60's when the bound 
states near magnetic impurities in s-wave superconduc- 
tors were predicted in a pioneering work of Yu, Sh iba and 
Rusinov llRiisinovl Il969t IShibaL Il968t IyH ll96fiD . They 
considered pairbreaking by a single magnetic impurity 
in a superconductor, and found that there are quasipar- 
ticle states inside the energy gap that are localized in 
the vicinity of the impurity atom. The corresponding 
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gap suppression occurs locally and the concept of lifetime 
broadening is inapplicable. In general, in this situation 
it is more useful to focus on local quantities, such as lo- 
cal density of states (LDOS), local gap etc., rather than 
on average impurity effects (which vanish for the single 
impurity in the thermodynamic limit). Yet it is clear 
that this local physics at some finite concentration of im- 
purities suppresses superco nductivity completely. This 
connection was discussed in usinovl Il969t IShibal . 11968^ 
lYulll965|) . In particular, formation of the intragap bound 
state and impurity bands due to magnetic impurity leads 
to filling of the superco nducting gap, and there f ore co n- 
nects to the AG theory l)Abrikosov and Gorkoviri96Cl() . 

At the time there were no experimental techniques 
to directly observe single impurity states. As a result 
the entire subject was largely forgotten until the STM 
was applied to study the impurity states by Yazdani et 
al. l)Yazda,ni p.t all [imS). This reinvigorated the field 
and lead to a firm shift in the interest from global to lo- 
cal properties. Soon afterwards STM was used to observe 
local impurity states near vacancies and impurities in the 
high-re cuprates jHu d son e t ai, 2001, 1999; Pan et (3, 
l2000bt lYazdani et all Il999j) . These discoveries opened 
a new field of research where impurities open a window 
into the study of electronic properties of exotic materials 
with atomic spatial resolution. As a first test of theories 
this allowed a direct comparison of the local electronic 
features in tunneling characteristics with the theoretical 
predictions for the density of states. 

We start by briefly reviewing the BCS theory in Sec.lTTI 
Our main goal there is to review three approaches that 
will be used to analyze the impurity effects: direct di- 
agonalization of the hamiltonian via Bogoliubov-Valatin 
transformation, variational wave function of the original 
BCS paper, and the Green's function method which is 
well suited to the analysis of multiple impurity problems. 
Then we define different types of impurity scattering in 
Sec. mil We pay special attention to distinguishing be- 
tween magnetic and non-magnetic impurities, and dif- 
ferentiating between static and dynamic scatterers. The 
basic features of non-magnetic scattering in s-wave su- 
perconductors are outlined in Sec. IIVI 

To keep in tune with our intention to make the review 
readable by graduate students and researchers entering 
the field, we begin the discussion of the localized states 
by considering an example of an impurity bound state in 
a two-dimensional (2D) metal in Sec.0 Then we discuss 
the low-energy bound state in s- and c?-wave supercon- 
ductors in Sec. IVII and Sec. IVIII respectively. Changes 
in the ground state of a superconductor containing a 
classical spin as a function of the coupling strength be- 
tween the spin and conduction electrons are discussed in 
Sec.HXI 

We proceed to consider the situations when the impu- 
rities have their own dynamics, so that their effect on 
the electrons is complicated, see Sec. IXII and the com- 
bined influence of the collective modes and impurities. 
Sec. IXIII We briefly touch upon possible existence of 



impurity resonances in different models of the pseudo- 
gap state of the cuprates in Sec. I Villi and discuss recent 
STM measurements on both conventional and unconven- 
tional superconductors in Sec. IXIIII The final two parts 
of our review are devoted to the discussion of the effects 
on impurities on meso- and macroscopic scale. For com- 
pleteness, we briefly review the basics ideas of computing 
the average density of states for a macroscopic sample in 
Sec. IXIVI For lack of space we cannot do justice to this 
very rich subject and use it largely to discuss new results 
on the impurity effect on the scales small compared to 
the sample size, but large relative to the superconduct- 
ing coherence length. In that situation there are dramatic 
consequences of local impurity realizations that may be 
different from the average, and we overview the results 
for the density of states in Sec. IXVI We conclude with 
the summary in Sec. IXVII 



D. Other related work 

In focusing largely on the properties of impurities on 
atomic or mesoscopic scales, we cannot give due atten- 
tion within the confines of this review to several other 
questions that have been important in the studies of im- 
purities. One of these is how exactly does the impu- 
rity band grow out of bound states on individual impu- 
rity sites, i.e. what is the effect of interference between 
such sites is in real space. We briefly review some of re- 
cent work in Sec. IXIVI but do not discuss the subject in 
depth. The answer to this question is still somewhat un- 
settled ev en in general: while t he usual finite lifetime 
approach dCorkov and Kaluginl Il985l: iHirschfeld et all 
\198A ISchmitt-Rink et a/.Ul986j) gives a constant density 
of states at the Fermi level in a d-wave superconduc- 
tor, and even though the same result has been obtained 
in field theoretical models of Dirac fermions in two di- 
mcnsions, mimicking the d-wave superconductor ijZiegleil 
|1996; Zicgler et ai, 199 6^, it has a lso been argued that 
this DOS diverges (Pepi n and Led. 11998. 2001), or van- 
ishes. Vanishing DOS can occur with different power laws 
in energy depending on the approach one uses to treat 
disorder llNersesvan and T'svclik. 119971 iNersesvan et al 



1995 tlSenthil and Fisheil.1999.) (see also^Bhasee^T^I 



20011) '). The vanishing itself can be traced to level re- 



pulsion whe n the system is treated within random ma- 
trix theory l|Altland and Zirnbaueilll997.) . Detailed self- 
consistent numerical studies indicate, however, the the 
behavior of the DOS depends on the details of the impu- 
rity scattering and electronic structure ( A tkinson et al\ . 
I2OOOI: IZhu et ai . 2GGGb). In particular, the divergence 
only occurs in perfectly particle-hole symmetric systems, 
and generically Atkinson et al. find that there is a non- 
universal suppression of the density of states over a small 
energy scale close to the Fermi level. 

The interference be tween many impuritie s have been 
inves t igated recently (jAtkinson et all 120031: IZhu et all . 
I200a l2004hD with the eye on the importance of these 
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effects for the interpretation of the features in the 
STM data on the high-T^ cuprates collected over a 
large area of the sample. The interference is also 
responsible for the formation of the impurity bands 
and therefore is crucial for determining the trans- 
port properties, which we do not address in this re- 
view. Within the framework of the ^-matrix approxi- 



mation transport pr operties of unconventional supercon- 
ductors in gene r al llArfi et all flOsS 
Hirschfeld etdl Il98fll Il98fil Il988t [ Pethick and PinesL 



^ entional sup ercon- 
iGraf et alL Il99f : 



II. A BCS THEORY PRIMER 

We begin by reviewing the Bardeen-Cooper-Schrieffer 
(BCS) theory. This section only briefly summarizes 
the results pertinent to our discussion; many excel- 
lent textbooks provide an in-depth view on the theory 
llde Ge nneij'. '198^, 'Kette rson and Sond Il999l: ISchrieffeil 
1X9641: rTinkham. 1996). Consider a general hamiltonian 

Ti-Bcs = ^o(r) + H,nt, where 



in particular dPuffv al.. 2001: ICraf et all 
[ Hirschfeld and Goldenfeld. 199,'i; .Hirschfeld et all 



1995 
1994, 



iSchmitt-Rink et oil Il98(f~an d high-T. c u prates ijo(r) = V / d'^n/'t (r)[£(_iVr) - dV'a(r) (2.1) 

rticular jDuffv et al. 2001: ICraf et all Il995l 

is the band hamiltonian of quasiparticles with dispersion 
e(k), /i is the chemical potential, and the interaction part 

i^^nt = -^E / rf'r(i'^rVt(r)V,t(r')K.^^5(r,r')^^(r')V'5(r). 



: iQuinlan et all ll996L 11993) have been extensively 
discussed, and the experiments on both microwave, opti- 
cal, and thermal conductivity are us ed to extract prop- 
erties of impurity scattering, see ijTimusk and Stam 
Il999j) for a revie w as well as very recent results in 
both experim e nt ifCarr et al]. [2000; Chi ao et al., 20 0C ; 
Corson et all l200(]t [Hill et all \2004 iHosseini et al. 



Turner et «iJ. I200.1D and theorv ([Berlinskv et al. 



Chubukov et a2l 1 2003 
Howell et all \200 



1999HLee eit_a/.i 12004 Isegre et a/.[. |2002HTu et a;J .[200;: 

J, 200C; 



Hettler and Hirschfeld, 



200C 
199E 



■Nicol and Carbotte. .2003.). The 



quest ion of localiz ation in both s-wave dMa an^ Lee[ . 



1985 ) and d-wave f Atkinson and Hirschfeld, 2002; Lee, 
1993t ISenthil and Fisher ,2000.: Senthil et gl. 1998; 
Vishveshwara et al. L 120001: lYashenkin et a,ll l200ll) " con- 
tinues to be investigated. Some of these results have 
been summarized in recent r eviews on high-Tc sys- 
tems l|Timusk and Statl Il999) . We also do not touch 
upon the rich phenomena related to the surfaces play- 



the formation of the bound states (lAorili et all 


1998 
1981 


Blonder et all 


19821: 


Buchholtz and Zwicknaerl 


Covington et al. 


,ll997HFogelstrom et aZJ.Il997HHul 


1994 


Kashiwava and Tanaka 


.1200(1. 



By now there are also few re views availab le on the 
subject of impurity states. Joynt l[ Jovntl Il997[) reviewed 
early work on the impurity states within the t-matrix 
theory focusing on anomalous transport due to finite life- 
time of the quasibound states around impurities. By- 
ers, Platte and Scalapino , contributed substantially to 
studies of the detailed electronic str ucture of the reso- 
nance state and interference patterns ("Bvers et al!, '1993*, 
[Platte and Bycrs, 1997a b; Platte and B vers, 1998), and 
reviewed their and related work ijPlatte and Bversl 
I1999D . An excellent review of thermal and transport 
properties of low-energy quasiparticles i n nodal super- 
conductors was recently given by Hussey ([Hussevl 1200^ . 

The subject is so rich and well developed that it does 
not seem possible to do justice to both local quasiparticle 
properties around a single impurity site and the questions 
of interference and transport within the confines of a sin- 
gle paper. With this in mind we now are ready for a main 
discussion. 



(2.2) 

Here r is the real space coordinate, a and /3 are the spin 
indices, and and tp a-re the fermionic creation and an- 
nihilation operators respectively. The mean field approx- 
imation consists of decoupling the four-fermion interac- 
tion into a sum of all possible bilinear terms, so that 

H.nt - E/ o''rd'^r'|y„^(r,r')V'i(r)V/3(r') (2.3) 

+A„^(r, r')Vl(r)V^^(r') + A:^(r, r')Mr)Mr') 

The effective potential, ^^^(r, r') is the sum of the 
Hartree and Pock (exchange) terms, and the last two 
terms account for superconducting pairing. The pairing 
field. A, is determined self-consistently from 

A„0(r,r') = iK0^5(r,r')(V7(i-')V'5(r)). (2.4) 

The pairing occurs only below the transition tempera- 
ture, Tj,; above the average of the two annihilation op- 
erators in Eq. (|2.4|l vanishes, and therefore Aq,^ = 0. In 
contrast, Hartree and Pock terms are finite at all temper- 
atures, and can be incorporated in the quasiparticle dis- 
persion, e(k). These terms do change below Tc, upon en- 
tering the superconducting state. Their relative change, 
however, is of the order of the fraction of electrons par- 
ticipating in superconductivity, and therefore is small for 
weak coupling superconductors A/W <C 1, where W 
is the electron bandwidth). Therefore the effective po- 
tential, V , is not explicitly included in the following dis- 
cussion except where specified. 

Therefore we start with a reduced mean field BCS 
hamiltonian, 

■Hbcs = E / d''riji{r)Ho{r)i;a.{r) (2.5) 
+ Y.f d'*rdV|A„^(r,r')Vi(r)^^(r') + /i-c.|. 
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The spatial and spin structure of AQ,^(r, r') determines 
the type of superconducting pairing. In most of this re- 
view we consider singlet pairing, when A has only the 
off-diagonal matrix elements in spin space, and it is com- 
mon to write AQ^(r,r') = (icr^)Q^A(r, r'), where A is 
now a scalar function, see previous section. 

In a uniform superconductor the interaction depends 
only on the relative position of the electrons, so that 
V{r,r') — V{p = r — r'). Therefore in the absence 
of impurities, the structure of the order parameter in 
real space depends on the symmetry properties of V{p). 
These are easier to consider in momentum, rather than 
coordinate, space. In models with local attraction, when 
V{p) = VoS{p), the Fourier transform of the interaction 
is featureless, and A(k) = Ag; an example of an isotropic, 
or s-wave superconductor. 

In the remainder of this section we overview the main 
methods solving the BCS hamiltonian since the same 
methods are commonly applied to the studies of impu- 
rity effects in superconductors. The approaches that we 
consider are: a) direct diagonalization via Bogoliubov- 
Valatin transformation; b) variational determination of 
the ground state energy from the trial wave function; 
and c) Green's function formalism. 



Aafj, which can be obtained directly from Eq. (|2.4() to be 



A„^(r,r') = -iy„^^^(r,r')^ 



(2.9) 



Here the Fermi function f{E) = [exp{E/T) + 1]-^. 

In a uniform superconductor the Fourier transform of 
the Bogoliubov equations, Eqs. (|2.7|) - (|2.8|) . into the mo- 
mentum space gives 

(ek-i^k)"ka + A„;3(k)«k/3 = 0, (2.10) 
(a + £^k)«ka + A;^(-k)uk0 = 0, (2.11) 

where is the bare quasiparticle energy, measured with 
respect to the chemical potential, ^k — e(k) — /i. In a 
singlet superconductor 

(a - £^k)ukT + A(k)^;ki = 0, (2.12) 
(^k + Ek)vki - A'^(k)MkT = 0, (2.13) 

and recover the familiar energy spectrum E'k = 
■y^^ + |A(k)p, with the coefhcients u and v given by 



A. Bogoliubov transformation 

Since the effective hamiltonian of Eq. (|2.5|) is bilinear 
in fermion operators, V' a-nd V'^j it can be diagonalized 
by a canonical transformation of the form 



(r)7n + t^na(r)7l 



(2.6) 



subject to condition |u„Q(r)p -I- |'i;„Q(r)p = 1. The re- 
sulting equations on the coefficients u and v are 

Eu^{v) = i/o(r)ua(r) + j d'^r' A^p{v,v')vp{T'\^7) 
-Ev^ir) = H*{rW{r)+ [ dVA*^(r, r')«/5(r'02.8) 



Here we suppressed the label n for brevity. Clearly, when 
A = 0, coefficients u and v do not couple, and there is 
no particle-hole mixing. 

For each n there are four functions, 
M|(r),it|(r),i;|(r),U|(r) that need to be determined. 
However, for a singlet superconductor, for example, 
the matrix A^^ is off-diagonal in the spin indices, so 
that U| (u|) couples only to vi (w|), so that in practice 
only two of the equations are coupled. In the presence 
of the impurity potential, however, in general all four 
components become interdependent. 

Equations H2.7|) - (|2.8|l . are coupled integro-differential 
equations for the functions Una (r) and Vna (r) . They have 
to be complemented by the self-consistency equations on 



1 ± 



E^ 



(2.14) 



B. BCS variational wave function 

Superconductivity originates from the instability of the 
Fermi sea towards pairing of time-reversed quasiparticle 
states. Therefore a variational wave function approach, 
originating with the classic BCS paper, is to restrict the 
trial wave function to the subspace of either empty or 
doubly occupied states, 

l*(r))=n(«" + ^«4T4i)|0), (2.15) 



and to minimize the energy, Ebcs ~ This is 

a legitimate approximation at T = 0, and is a very good 
approach at low temperatures. In Eq. (|2.15() the vacuum 
state |0) denotes the filled Fermi sea, and c|^^ {'^ni) 
ates a quasiparticle with spin up (down) and with the 
wave function (/'„(r) ((/)*(r)) that is the eigenfunction of 
the single particle Hamiltonian. Normalization requires 
that |a„|2 + |6„|2 = 1. 

In the absence of impurities these eigenfunctions can 
be labeled by the same indices, k and a, as in the pre- 
vious section. Consequently, the variational approach is 
completely equivalent to the Bogoliubov analysis with 
the choice M„(r) = a„(/;>„(r), and w„(r) = 5„(^„(r). In 
general, however, interaction with impurities may lead to 
the appearance of the single particle states in the ground 
state wave function, see Sec. IIXI Moreover, it is worth 
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remembering that energy of the state described by the 
BCS wave function is greater or equal to that of the 
exact ground state obtained by solving the Bogoliubov 
equations. 



C. Green's functions 

The third approach that we will use in this work is 
the Green's function method, which originates with the 
work of Gor'kov. Following Nambu we introduce a 4- 
vector that is a spinor representation of the particle and 
hole states, 



vI/t(r) = (^|,^UT,V'i) 



(2.16) 



The matrix Green's function is defined as the imaginary- 
time ordered average 



(2.17) 



where the four-vector x — (r, r) combines the real space 
coordinate, r, and the imaginary time, r. The time 
evolution of the creation and annihilation operators in 
the Heisenberg approach is given by dip/dr = [Hbcs — 

For a singlet homogeneous superconductor the Hamil- 
tonian of Eq. I|2.5|l in the Nambu notation takes the form. 



Hbcs - J dr*t(r)(^(-iV)r3 + ATia2)*(r), (2.18) 

and we find l)Makill969|) 

Go\-k,Lu) = iu;„ ~ e(k)T3 - A(k)t72Ti. (2.19) 

Here LOn — nT{2n + 1) is the Matsubara frequency, ai 
are the Pauli matrices acting in spin space, are the 
Pauli matrices in the particle-hole space, and Tiaj de- 
notes a direct product of the matrices operating in the 
4-dimensional Nambu space. The self-consistency equa- 
tion for a single superconductor takes the form 

A(k) = -TJ2 f dk'l/(k,k')Tr[Tia2Go]. (2.20) 

In BCS the interaction is restricted to a thin shell of 
electrons near the Fermi surface, and therefore 

A{h) ^ -TNo I dh'V{h, h') Tr tiCTs [ d^wGo 

l^n L 

(2.21) 

where D, denotes a direction on the Fermi surface, and 
iVo is the normal state density of states. 

The off-diagonal component of Go, is often called 
the Gor'kov's anomalous F, (Gor'kov) Green's functions 
since it describes the pairing average 



In general Fai3{x,x') = gapF{x,x'), where g is the ma- 
trix describing the spin structure of the superconducting 
order. For the singlet pairing g — ia^y\ where ct*^^-* is the 
Pauli matrix. Therefore in a singlet spatially uniform su- 
perconductor normal and anomalous components of Gq 
are 



G(w„,k) = 
F(c^„,k) = 



A(k) 



(zc.„)2 -^l- |A(k) 



(2.23) 
(2.24) 



The connection with the Bogoliubov's transformation 
is provided by rewriting the Green's functions as 



G(c^„,k) 



F(cj„,k) = ui^vl 



1 1 



iujn - El, iujri + Eu 



(2.25) 
,(2.26) 



where Wk and Wk are given by Eq. (|2.14() . 

The three approaches discussed above are complemen- 
tary and equivalent in the case of homogeneous super- 
conductors. However, some of them are better suited for 
addressing specific questions in the presence of impuri- 
ties. In particular, we will see that the Green's function 
method is sometimes advantageous for determining the 
thermodynamic properties of a material and averaging 
over many impurity configurations. For inhomogeneous 
problems, where we are interested in the spatial varia- 
tions of the superconducting order and electron density, 
both Bogoliubov equations and Green's functions are of- 
ten used. In general, the choice of a specific methods 
depends on the type of question asked in the presence of 
impurities, and we briefly describe the basic models and 
issues related to impurity scattering in superconductors 
below. 



III. IMPURITIES IN SUPERCONDUCTORS 
A. Single impurity potential 

If we are to address theoretically the question of what 
defects do to superconductivity, we must describe the 
defects and superconductivity in the same framework. 
Grain and surface boundaries, twinning planes, and other 
structural inhomogeneities scatter conduction electrons 
and therefore affect the resulting order parameters. How- 
ever, here we focus on only one type of imperfection: im- 
purity atoms. 



a. Potential scattering. First and foremost an impurity 
atom has a different electronic configuration than the 
host solid, and therefore interacts with the density of 
conduction electrons via a Coulomb potential. 



Faf3{x,x') = -{Trll^a{x)lp/3{x')). 



(2.22) 



Hi, 



E 



drV'l(r)C/pot(r)-0Q(r 



(3.1) 
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In good metals the Coulomb interaction is screened at 
the length scales comparable to the lattice spacing, and 
therefore the resulting scattering potential is often as- 
sumed to be completely local, C/pot(r) — UQS{r — Tq), 
with the impurity at Tq. The resulting scattering occurs 
only in the isotropic, s-wave, angular momentum chan- 
nel. If finite range of the interaction is relevant, scatter- 
ing in / 7^ channels needs to be considered. In that case 
the treatment is similar to that of magnetic scattering in 
conventional superconductors, see Sec. IVII and was ap- 
plied to unconventional superconductors in, for exa mple 
l|Balatskv eF^ Il994 Ikampf and DevereauxL Il997|) . 

In the 4- vector notation of the previous section the po- 
tential scattering has to have the same matrix structure 
as the chemical potential, or e(k) in Eq. (|2.19() . so that 



Hi, 



ydrvI't(r){/p„t(r)T3M'(r), (3.2) 



or, in Nambu notation, 

Upot = UoT36{r - ro) 



(3.3) 



b. Magnetic scattering. In addition to the electrostatic 
interactions, if the impurity atom has a magnetic mo- 
ment, there is an exchange interaction between the local 
spin on the impurity site and the conduction electrons. 



Himp = X! / '^rJ(r)V4(i")S ■ crQ/3-0/3(r). 



(3.4) 



The range of interaction here is determined by the quan- 
tum mechanical structure of the electron cloud asso- 
ciated with the localized spin. Again, in reality we 
often consider a simplified exchange hamiltonian with 
J(r) = JoS{r — ro), which captures the essential physics 
of the problem. In the 4- vector notations of the previous 
section the electron spin operator becomes 



1 

"=2 



(1 + T3)cr + (1 - T3)a3(Ta3 



Therefore 



or, in Nambu notation, 

Umag = J(r)S • «. 



(3.5) 



(3.6) 



(3.7) 



c. Anderson impurity. However, even if the ground state 
of an isolated impurity has an electron spin, the result of 
putting such an impurity into a host matrix may modify 
the spin configuration as the impurity electrons couple to 
the conduction band. Therefore a realistic model for an 



impurity site is based on the Anderson model, with the 
Hamiltonian 

Ha = ^£;o4c?a + t^"dT'^di + ^.sd) (3-8) 

a 

Hsd = 51 ^''d^.a'^" + /i.e. (3.9) 

k,a 

Here Eq is the position of the impurity level relative to 
the Fermi energy, and d operate on the impurity site, 
U is the Coulomb repulsion for the electrons localized on 
the impurity site, and cj^, Ck create and annihilate the 
conduction electrons. This Hamiltonian allows the elec- 
trons to hop on and off the impurity site, resulting in a 
finite width of the impurity level, F = 7r|ysdpA^o- The 
model describes the potential scattering, when J7 <C F. 
On the other hand, when Eq ^ Ep, Eq + U '3> Ep, and 
[/ ^ F, we expect the local levels to remain split, so 
that the impurity state is singly occupied and has a local 
spin. Therefore it allows a natural interpolation between 
potential and magnetic scattering, as well as the study 
of the mixed valence regime. The price to pay for such a 
rich behavior of the Anderson impurities is the difficulty 
of studying them analytically, and therefore in practice 
many results have been obtained in the simplified mod- 
els above, although a number of very thorough numerical 
renormalization group studies of Anderson impurities in 
superconductors exist. We will review some of them for 
completeness, but will not focus on those extensively. 

B. Many impurities 

In all of our discussions we assume noninteracting im- 
purities, so that the net impurity potential is simply 



^imp(^) — ^ ^ Ujmp ^i) 

dr'p,mpir')Uimp{r - r' 



(3.10) 
(3.11) 



Here U denotes the matrix structure of the potential in 
both spin and particle-hole space, and we introduced the 
impurity density. 



(3.12) 



We also assume the dilute impurity limit of the average 
impurity concentration <C 1, where 



(3.13) 



For magnetic scatterers it was explicitly shown that 
the effect of the RKKY interaction between scatter- 
ing centers on the superconducting properties is small 
ijGalitskii and L arkin. 2002; Larkin e t aZ.l . ll97i(l . 

If we now compute a local physical quantity, such as 
the density of states measured at the position r by the 
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STM, it will depend on the distance from the nearby im- 
purities, and therefore will be different for different real- 
ization of impurity distributions. In contrast, thermody- 
namic quantities, such as the density of states measured 
in planar junctions, or the specific heat, average the den- 
sity of states over many random local configurations of 
impurities. Therefore in computing their values we av- 
erage over all impurity configurations (^Abrikosov et all . 
Il963l) . so that, for example. 



G'(w„,k)=J| — / drjG(w„,k, ri,...,rjv. 



(3.14) 



Here a bar denotes such an impurity average. 

By definition pimp = n-i. We also assume an uncorre- 
lated, or random, impurity distribution, which means 



Ni 



p(r)p(r') 



= ni5{Y-Y')+nl. (3.15) 

Since the impurities are dilute, <C rii, and we neglect 
the second term compared to the first. In Sec. IXIVl we 
implement this impurity averaging procedure to deter- 
mine the average density of states. 



C. The self-energy and the T-matrix approximation 

In practice to compute the Green's function in the pres- 
ence of impurities we will often employ the T-matrix ap- 
proximation. This method is described in detail in man; 



original articles and re views ("Hirschfcld and Goldenfeld', 



1993; 



2002; 



Hirschfeld et alV 11986. .1988- Hotta. .1993; Hussevi. 
Mahanll200(i|r and we only briefly summarize it. 



For a single impurity with the scattering potential 
f k,k' in the momentum space (given by one of the models 
discussed at the beginning of this chapter), the T-matrix 
accounts exactly for multiple scattering off of one impu- 
rity. In the language of Feynman diagrams, the corre- 
sponding process is shown in Fig. ^ Here, and through- 
out the review, the hat over a letter means that it denotes 
a matrix in Nambu space. Therefore the full Green's 
function is 

G(k,k') = Go(k) + Go(k)&k,k,Go(k') (3.16) 
+ ^Go(k)&k,k"Go(k")C/k",k'Go(k') + . . . . 

k" 

Here we suppressed the frequency index in the Green's 
function as the scattering is elastic. The series can be 
summed to write (see Fig. ^ 

G(k, k') = Go(k) + Go(k)fk,k'Go(k'), (3.17) 
where the T-matrix is given by 

fk,k' = f/k,k'+5It^k,k"Go(k")^/k",k' + ...(3.18) 

k" 

- t/k,k'+^C/k,k"Go(k")fk",k'. (3.19) 



This equation needs to be solved for T. If the impu- 
rity scattering is purely local, C/(r — r'), the scattering 
is isotropic, ?7k,k' = U , greatly simplifying the process of 
solving the equation for the T-matrix, as T depends only 
on frequency. 

Notice that we could draw the set of diagrams in Fig.Q] 
in real space, and write the corresponding set of equa- 
tions for the T-matrix and Green's function G(r, r') in 
complete analogy with Eq. H3.19|l . The main observation 
here is that, in the vicinity of the impurity, the trans- 
lational invariance is broken, and the Green's function 
depends on two momenta, k and k'. 

G(r, r'; w) = Go(r, r'; w) + Go(r, Tq; w)f (w)Go(ro, r'; lo) 

(3.20) 

The T-matrix lends itself easily to describe the effect of 
an ensemble of impurities. The so called self-consistent 
T-matrix approach (SCTM) considers multiple scattering 
on a single site of an electron tha t has already been scat- 
tered on all other impurity sites ijHirschfeld et al\ . Il986l 
Il988|) . This results in replacing the bare Green's func- 
tion in Eq. (|3.19|) by its impurity-averaged counterpart, 
G(k, w). Notice that after averaging over the random 
impurity distribution the translational invariance is re- 
stored, so that the Green's function depends on a single 
momentum k. The combined effect of impurities is given 
by the self energy, I](k, w) = riimpTk^k, so that 



G-^{k,uj) = Go \k,tj) - E(k,tj). 



(3.21) 



In contrast to the single impurity case where Eq. (|3.17|) 
with the T-matrix given by Eq. (|3.19|) is the exact so- 
lution of the problem, the Green's function given above 
is an approximation, and much recent research is moti- 
vated by questions about how accurately it describes the 
properties of nodal superconductors with impurities. 



D. Static and dynamic impurities 

So far we only discussed the static impurities, and 
most of our review addressed such a situation. However, 
even for purely potential scattering a situation is possi- 
ble when a vibrational mode leads to a time-dependent 
modulation of the charge on an impurity site, and, as a 
result, Upot acquires a characteristic frequency. Such a 
mode can be extended, as a phonon, or local. Influence 
of the dynamical impurity on the local properties of a su- 
perconductor is a relatively new subject of research and 
we summarize recent results in Sec. IXII 

For magnetic scattering the situation is more complex 
even in a normal metal. The degeneracy between the 
spin-up and spin-down states on the impurity site and 
the non-trivial commutation relations between different 
spin components ensure that quantum dynamics of the 
impurity is generated even if the exchange constant is 
purely static. In simple words, if the scattering process. 
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which flips both the spin of the conduction electron and 
the impurity spin, is relevant, the dynamics of the lo- 
cal spin flips becomes essential. This dynamics leads to 
Kondo screening of the impurity spin in a metal, and in 
Sec. IXII we briefly discuss the current status of the yet 
not fully understood problem of Kondo effect in a super- 
conductor. 

In the limit of large impurity spin, however, the change 
of the impurity spin by 1 during the spin flip scattering is 
not relevant, and its dynamics does not play a major role. 
In this limit of classical spin the static local and global 
density of states is discussed in Sec. ^3 and Sec. IXIVI 
respectively. Such a spin acquires dynamics only when 
placed in an external magnetic field, which also affects 
the superconducting state. 



IV. NON-MAGNETIC IMPURITIES AND ANDERSON'S 
THEOREM 

One of the most important early results was the ro- 
bustness of the conventional superconductivity to small 
concentrations of non-magnetic impurities. Theoretical 
unde rpinning of this r esult is known as Anderson's theo- 
rem l)Andersonlll959(l . Anderson's observation was that, 
since superconductivity is due to the instability of the 
Fermi surface to pairing of time-reversed quasiparticle 
states, any perturbation that does not lift the Kramers 
degeneracy of these states does not affect the mean field 
superconducting transition temperature. 

This is most clearly seen from the BCS analysis, which 
we carry out foUowing lMa and LeelHoSSl We consider an 
isotropic pairing potential, yo,/3^5(r, r') = VS{r — r'). In 
the absence of a magnetic field the coefficients a„ = sin On 
and bn = cos On can be taken real without loss of general- 
ity, so that the self-consistency condition, Eq. (|2.9|l reads 



equivalent to that of a pure superconductor provided the 
density of states 



(£2 _^ ^2 \ 



d'^r02(r)0^(r). (4.1) 



Here 



A„ = d'^rA(r)0^(r). 



(4.2) 



As noted above, in the BCS approach (/)'s are the eigen- 
functions of the single particle hamiltonian. In the ab- 
sence of impurities the system is translationally invariant, 
so that A(r) = A„ = Aq. The most important assump- 
tion underlying Anderson's theorem is that the super- 
conducting order parameter can be taken to be uniform, 
A(r) = Ai, even in the presence of impurities. In that 
case the individual eigenfunctions of the single particle 
hamiltonian including impurities are rather complicated. 
However, the gap equation, Eq. H4.1|l . takes the form 



iV(6,r) = ^</)2„(r)<5(e 



(4.4) 



V J yi^TAf' 



(4.3) 



is unchanged compared to the pure metal, A^(e,r) w po- 
ll this condition is satisfied, the solution Ai of the gap 
equation Eq. H4.3|) must be identical to that of the BCS 
equation in the absence of impurities, and therefore the 
transition temperature and the gap are insensitive to the 
impurity scattering at the mean field level. 

Anderson's theorem helped explain why superconduc- 
tivity was robust to impurities in many early experi- 
ments. It is important to realize however that it is an 
approximate statement about the thermodynamic aver- 
ages of the system. Beginning with the next section we 
will analyze in more detail the changes that impurities 
create in superconductors in their immediate surrounding 
as well as on average. We will see that even purely poten- 
tial scattering does induce changes in the local properties 
of superconductors, albeit the corresponding change in 
the transition temperature remains minimal. Anderson's 
theorem brings to the fore the need to separate the study 
of impurity effects on different lengths scales, from lattice 
spacing to the coherence length, to sample size. 

Before we proceed to study the local properties we dis- 
cuss the extensions of the Anderson's treatment of im- 
purities. In weakly disordered systems the density of 
states remains nea rly constant as a function of disor- 
der. Ma and Lee l|Ma and Led [198 5*) argued that An- 
derson's theorem remains valid in the form above even 
in a strongly disordered system provided the localization 
length, L ^ (pqAq)"^/'^. In that case there is a large 
number of states localized within energy Ag of the Fermi 
surface, and these state form a local superconducting 
patch. The Josephson interaction between the patches 
then leads to global phase coherence at T = 0. More- 
over, they argued that the theorem holds all the way to 
the limit of site localization. 

It is important to note that the supcrfluid stiffness, i.e. 
the ability of the superconductor to screen out the mag- 
netic field, is affected by disorder. In particular, when 
the quasiparticle lifetime, r, becomes sufficiently short, 
Aqt <C 1, the superfiuid density ps ~ Aqt. Conse- 
quently the superconductor is sensitive to the local phase 
fluctuations of the order parameter, and the experimen- 
tally observed transition temperature may be severely 
suppressed compared to the mean field Tc, as it is, for 
example, in granular superconductors. Approaches ex- 
tending beyond the mean field picture are largely outside 
the scope of this review. 

Therefore for dilute impurities Anderson's theorem 
is valid provided the superconducting order parameter 
can be taken to be nearly uniform. Since the "heal- 
ing length" of A(r) over which it can change apprecia- 
bly is the coherence length, — ^i^f/Aq, where vp is 
the Fermi velocity, while the Coulomb screening length 
for the charged impurities in metals is of the order of 
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the lattice spacing, a, for ^ a the order parame- 
ter remains essentially uniform, and Anderson's theo- 
rem holds. Much work has been done recently on the 
effect of the ultrashort coherence length on the impu- 
rity scattering in superconductors. In particular, it has 
been shown that when the superconducting pairing is of 
the order of the electron bandwidth, Anderson's theorem 
is violated JCh osal et al, 1998; Moradian et al., 2001; 
iTanaka and Marsigliol I2OO0I) . 

Ghosal et al. ("Ghosal et al}, |199^, and Xiang and 
Wheatley (Xiang and Wheatlevi. Il99,'j) have explored in 
detail the discrepancy between the single particle exci- 
tation gap and the superconducting order parameter as 
a function of disorder in these circumstances. Beyond 
Anderson's regime of the constant density of states, both 
quantities decrease at first, since the disorder depletes 
the density of states. Then, however, the spectral gap 
persists while the supercond ucting order vanishes. As 
pointed out by Ma and Lee l)Ma and Lee . 1985) in the 
limit of strong disorder the models with on-site pair- 
ing, such as those studied by Tanaka and Marsiglio, and 
Ghosal et al., show the on-site spectral gap (so-called 
Anderson negative-C/ glass) without the off-diagonal long 
range order and without symmetry breaking. 

In most experimentally relevant situations, however, 
the corrections to the main statement of Anderson's the- 
orem are quantitative rather than qualitative. This is 
generally true of most results pertaining to the impurity 
scattering in superconductors, and therefore it is very 
instructive to consider this problem in BCS-like systems. 



V. SINGLE IMPURITY BOUND STATE IN 
TWO-DIMENSIONAL METALS 

Before we proceed to calculate the effect of impurity 
in a d-wave superconductor it is instructive to review 
a simpler problem of a single impurity in a metal. We 
show here a T-matrix calculation for finding the bound 
states due to a single impurity in d dimensions with an 
attractive potential Uo < 0. The Hamiltonian for the 
problem is 



H 



Ck', 



(5.1) 



k,k' 



the Uq term describes the on-site energy change of elec- 
tron density n{r) in external potential U{r) = Uod{r). 
We consider a single particle (/i = 0), although the cal- 
culation for the normal metal with a finite density of 
states follows simply by replacing e(k) — > ^{\bfk) in the 
following. 

The bare Green's function for a free particle is 



(5.2) 



Since the vertex of the impurity interaction, Uq is mo- 
mentum independent, the equation for the T-matrix is 



particularly simple and follows from Eq. (|3.19|) . 



Uo 



l-C^oEkGo(c^,k) 



(5.3) 



Summation over momentum in is easily performed 
using the density of states 



iV(e)=5]^(6-6(k))^r,6^-\ 



(5.4) 



where is a constant dependent on dimension. There- 
fore 



90 



{uj) = V Go (t^, k) - / ^^^^^ ~ -TdUJ^ , (5.5) 
"'0 



deN{e) 
10 w - e 



for d 2, where W is the electron half bandwidth. In 
two dimensions go ~ — r2 ln(T4^/|a;|). Consequently, the 
T-matrix for d 7^ 2 is given by 



T : 



Uo 



1 - .gd^" 



(5.6) 



where gd — —Uo^d is the effective coupling constant, and 
by the same expression with the obvious substitution of 
ln{W/uj) for d = 2. 

Since the Green's functions in the presence of impurity 
scattering is G = Go + GqTGo, see Eq. 1)3.17(1 . poles of 
the T-matrix are the new poles of G that are not poles 
of Go, signifying the appearance of new states. We can 
find this pole, ujo, from Eq. 1(5. 6() for different dimension 
d. The bound state {ujo < 0, see Fig. is formed for an 
arbitrarily small potential \Uo\ in d = 1,2, but requires 
a critical coupling for d = 3. The energy of this state is 
given by 



1 



ifd=l; (5.7) 
ifd=2; (5.8) 

^0 ^ [gi - g|c], 93 > 93c, if d = 3, (5.9) 



ujo = W^exp( ), 

92 



where the d = 3 critical coupling g^c ~ W~^^'^. 

We focus in more detail on the two-dimensional case, 
when 52 = r2|C^o| and T2 = is the electron density 

of states. The bandwidth, W ~ „ is the ultravio- 
let cutoff corresponding to the lattice parameter a for 
free particle. This result can be compared to the solu- 
tion of the Schrodinger's eq uation for the parti cle in the 
2D attractive potential Uo l|Landau and Lifshitj . 12000). 
Ch. 45. For an arbitrary potential U{r) the solution 
obtained using the T-matrix is asymptotically correct 
if the scattering length is greater than a. For shallow 
potential the bound state energy — wq is small, and the 
characteristic extent of the bound state wave function is 
lo = { ^mt^ Y^"^ ^ Therefore in this limit we can safely 
approximate U{v) = UoS{r), where Uo = J f/(r)dr. 
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Finding the energy of the bound state, Eq. H5.8|l . is 
only one part to the solution. We also want to determine 
the corrections to the local Density of States due to bound 
state. We write the equation for the Green's function in 
real space, Eq. (|3.20() . 

G(r, r';u) = Go(r, r'; lu) + Go(r, 0; lo)T{lu)Go{0, r'; lo) 

and read off the position dependent Density of states 
(DOS) 



N{r,uj) 



1 



/mG(r, r; u) 



Noir,Lj)+SN(r,uj). (5.10) 



Here the first term is the DOS of a clean system, and 
the second is the correction due to the bound state. We 
focus on the energy range close to the bound state energy, 
w ~ wq. Since the bound state is below the bottom of 
the band, the unperturbed Green's function G° has no 
imaginary part in this range (A^o = ^Im.9o(^ = 0)/7r). 
Therefore the only contribution to the imaginary part 
of the full Green's function, Eq. (|5.1U|) comes from the 
T-matrix : 



ImT(a;) 



Im 



1 



= Im In 



1/52 " log[W/{~Lu)] 
(jj + iS. 



TrS{uj - loq), 



(5.11) 



and the correction to the DOS of a clean system is : 



6N{r,iu) = \Goir,uo)\^Siuj-iUo) 



(5.12) 



with Go(r,cj) = NoJf){kpr) ln[^] is the real part of 
Green's finction in 2D systems. Equations (|5.9() and 
(|5.12|) are the main results of this section. They establish 
the logic we will adhere to in finding impurity induced 
bound states: a) find the poles of the T matrix in the and 
the poles of the dressed Green's function Eq. (|5.9() . b) find 
the inhomogeneous DOS due to impurity induced state, 
Eq. (|5.12l) . One should keep in mind that this approach 
is just one of many one can implement in a search for 
scattering induced bound states. Alternatively one can 
use for example the exact numerical solution of a finite 
system. As we will argue for superconducting case the 
self-consistency condition can not be implemented an- 
alytically and the numerical solution remains the only 
method available. 



VI. LOW-ENERGY STATES IN s-WAVE 
SUPERCONDUCTORS 

A. Potential scattering 

Even though the potential scattering does not change 
the bulk properties of isotropic s uperconductors, it does 
affect the local density of states l)Machida and ShibataL 



119721: IShibal 119731) . Let us consider the Anderson im- 
purity model, Eqs. (|3.8|) - (|3.9|l in the hmit U = (res- 
onance scattering). As discussed above the localized 
state acquires a finite width, T = TT\Vsd\'^No, due to hy- 
bridization with the conduction band. The Green's func- 
tion of the conduction electrons can be evaluated in the 
T-matrix approach, with the result at real frequencies 
ifMachida and ShibataK .1972; .Shiba. .1973,1 



T{lu) = \Vsd\^T3 



'-EoT3-\Vsd\^T3Y,GoiK^)r3 



T3- 



(6.1) 

The poles of the T-matrix determine the location of the 
bound states 



2r 



(6.2) 



In most physical situations F 3> A, so that the bound 
states are located at 



(6.3) 



where Nd{0) = 7r~'^r / (T^ + E^) is the density of states of 
the resonant impurity level. For typical densities of states 
ANd{0) ^ 10^"^, so that the bound states lies essentially 
at the gap edge. Shiba considered a finite but small value 
of the Coulomb repulsion and allowed for the induced 
pairing on the impurity site (jShiba, 1973). He concluded 
that, even though there may be a shift of the bound 
state to lower energies, it still lies within lO^'^A of the 
mean field gap edge, and therefore can be neglected in 
the discussions of physical properties. 



B. Classical spins 

If the substitution atoms have a magnetic moment, the 
time-reversal symmetry is violated, and therefore super- 
conductivity will be suppressed. We consider the mag- 
netic scattering, Eq. (|3.6|l . which we rewrite in the mo- 
mentum space as 



Elpx — 



(6.4) 



k,k' 



We first review a simplified version of this problem, where 
we do not need to consider Kondo screening. We re- 
view scattering on classical spins first studied indepen- 
den tly at ab out the same time by Shiba. Rusinov, and 
Yu l)R,usinovl [19681 119691: IShibaL 119681: IYuL 119651) . Quan- 
tum mechanical properties of spin can be neglected when 
S' — > oo, and we simultaneously take J — > so that the 
product JS = const. In this limit the localized spin acts 
as a local magnetic field. 

Therefore we study the effect of the impurity with the 



potential U{y) = Uo + Uex, or Himp 



Hi, 
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BCS s-wave superconductor with the unperturbed hamil- 
tonian of the form 

Ha = Y^ £k4 „Cka + Ao Xl't^kT^-ki + c-kiCktl- (6-5) 

k 



ka 



This problem serves as a starting point for aU subsequent 
analysis of the resonance states in superconductors. 

To find a localized state with energy < i? < Aq near 
a single parama gnetic impurity we perform a Bogoliubov 
transformation usmo ^ 119681 IYuL [1961 to find 

Eua{r) = e(k)u„(r) + iAcr^^W0(r) + J7a;3(r)w/3(r)(6.6) 
Ev^{r) = -e(k)w„(r) - iAal^up{r) - [/„^(r)w^(r)(6.7) 

This system is solved by Fourier transforming the equa- 
tions and expanding the impurity potential in spherical 
harmonics, and has solutions with energies 

^ l + jTrNoVif -iTrNoJiS/2f 

Ao V[l + i^NoVi)^ - {ttNoJiS/2)^]^ + 4(7riVo J/5/2)2 ' 

(6.8) 

where A^o is, again, the density of states at the Fermi 
energy in the normal state. This result can be written in 
a more elegant form if we introduce the phase shifts. Si, 
of scattering for up (+) and down (-) electrons, in each 
angular channel. 



tan(5± = {7:No){Vi±JiS/2). 
Then the energies of the states in the gap become 



(6.9) 



— = cos{6l 



(6.10) 



Clearly, for purely potential scattering {S^ = S^) the 
spectrum begins at the gap edge, and there are no in- 
tragap states. However, as the magnetic scattering in- 
creases, a series of low-energy states below the gap edge 
appear. Purely magnetic scattering corresponds to 6'^ = 
— i5j~, and strong scattering (large phase shift) yields a 
localized state deep in the gap, while weak scattering 
(small phase shift) results in the bound state very close 
to the gap edge. 

The same result ca n be obt ained us ing the Gr een's 
function formulation l)Rusinovl .1969: ShibaL Il968ll and 
solving the single impurity problem using the T-matrix 
method described above. With the impurity hamiltonian 
of Eq. H3.6|) in the Nambu notations the matrix Green's 
function for the system is 

G(k,k';tj) = Go(k,t^)5(k-k')+Go(k,w)f (k,k')Go(k',w) 

(6.11) 

Here the T-matrix is computed as in Sec. lIIII and we sum 
over the indices of the matrix a in each vertex. The l-th 
angular component of the T-matrix satisfies the matrix 
equation (for a spherical Fermi surface and isotropic gap) 



3.12) 



Tiiuj) ^Ui + Ui deGo{k,uj)Ti{u;). 



The full expressions for the T-matrix for both poten- 
tial and magnetic scatteri ng in all angul ar channels is 
straightforward to obtain l)R,usinovl Il969l) but is rather 
cumbersome, so that we don't give it here. Even the case 
of only spherically symmetric scattering (I — 0) with both 
t/p 7^ and J the T-matrix is simple yet lengthy 
l)Okabe and Naeil Il983tl . The main results for the en- 
ergy of the Shiba states remains the same, of course as 

Eq. inini- 

In the particular case of purely magnetic spherically 
symmetric exchange, J(k — k') — J, only I = compo- 
nents are non-vanis hing and the T-matrix has a partic- 
ularly simple form ijShibal Il968() . The diagonal in spin 
indices component is. 



T(i)(a;) 



1 (JSmoiu) 



N I -{JSgo{oj)l2Y- 
Here is the local matrix Green's function. 



-nNo- 



VAf^ 



The bound state energy 



eo = 



Ao 1 + {JSTTNo/2y 



(6.13) 



(6.14) 



(6.15) 



The wave functions of the bound states at Ei can be 
computed using the Bogoliubov equations above. In the 
simplest case of isotropic sc attering at distances r ^ pl^^ , 
both u{r) and v{r) vary as l|FetteillT965: Rusinov. 19691) 



sin(p_Fr - Sq) 



Ppr 



exp(-r/^o|sin(5o+-5o-)|, (6.16) 



that is, the state is localized near the impurity site at 
distances 



ro 



Co 



6 



|sin(J+-<5o)| yr 



(6.17) 



The square of these coefficients gives the spatial depen- 
dence of the amplitude of the particle and hole com- 

onents of the densi ty of states at a given position r 

lYazdani eroZI. Il997|) . 

The analysis above was carried out under the assump- 
tion that the variation of the superconducting order pa- 
rameter, A, around the impurity site does not change 
the position of the resonance low energy state. There 
are several characteristic length scales for this variation, 
(5A(r). Far away from the impurity, r 3> ^O: at tempera- 
tures close to Tc, where this variation can be determined 
perturbatively , (5A(r)/Ao ~ l/{pFr) (Heinrichs,, 19(^ 
lRusinovLll968() . This power law is insensitive to the phase 
shifts of scattering on the impurity. At low temperatures 
a fully self-consistent treatment is required, which leads 
to SA{r) decaying as {pFr)~^ and oscillating on the scale 
of ^qAq/loo, where the Debye temperature ujo sets the 
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scale for the interaction between electrons l)SclilottmannL 
[l97i). 

In the immediate vicinity of impurity, vf/ojd ^ <C 
^0, the variation of the order parameter is (5A (r)/An c± 
l£ipj^r)'^ in the linear response approximation <)Rusinovl 
[l968). In the fully self-consistent treatment at distances 
r <C S,qujd/ Ep, this dependence was found to acquire an 
oscillating factor sin^p^^r (,,Schlottma nn. 1976). 

In all these cases, since the suppression of the order 
parameter is determined by the Fermi wavelength, the 
effect is negligible in determination of the position of the 
bound state. 



VII. IMPURITY-INDUCED VIRTUAL BOUND STATES 
IN d-WAVE SUPERCONDUCTORS 

We are now ready to extend our discussion to impu- 
rity induced states in d-wave superconductors. Scalar 
(non-magnetic) impurities are pair-breakers for "higher- 
orbital-momentum" states, such as a d-wave pairing 
state. This occurs because change of the momentum of 
particles in the Cooper pair upon scattering disrupts the 
phase assig nment for p articular momen ta in a nontriv- 
ial pairing llAndersonl . ,1959 : Markowitz and Kadanofl 
iTsunetol 1 1962(1 . More rigorously this follows from 
the analysis of the normal and anomalous self-energies 
due to scatterin g within the Abrikosov-Gorkov theory 
l)Abrikosov et all . fl963(l . We also note that one of the 
first arguments about pairbrea king effects o f potential 
scattering was given by Larkin l|Larkinl Il965f) . 

As we have emphasized, for pairbreaking impurities 
the local properties of the superconductor near an im- 
purity site, such as the local density of states and the 
gap amplitude, will be modified dramatically. To cap- 
ture these modifications, we use a variation of the Yu- 
Shiba-Rusinov approach llR,usinovl[T968tlShibal.lT965lYl 
119651) . which treats magnetic impurities in the strong 
scattering limit, see Sec. I VII We restrict our consider- 
ation to the s-wave scatterers with the phase shift close 
to the unitarity limit, — 7''/2, when the bound state 
has energy away from the gap edge. In contrast to the 
s-wave superconductors, in d-wave systems the density 
of states below the gap maximum is non-zero, and varies 
linearly with energy in a pure system. Consequently, the 
overlap with the particle-hole continuum only allows the 
formation of virtual bound states with a finite lifetime. 

We focus in this section on point-like defects, and 
use the T-matrix approach. A closely related method 
uses quasiclassical approximation and picture of Andreev 
scattering i deas to rep r oduce the results of T-matrix 
calculatio n llChen et all Il999t IChoi and Muzikaii 1199(1 
IShnirma n et a^J . ll999j) . Even more interesting results are 
obtained within the quasiclassical formalism for extended 
defects. For example, it has been shown that index theo- 
rem d ictates the existence of the low energy quasi-bound 
state l|Adagideli et all \199^ . 

Zn substitutions in cuprates are one example of non- 



magnetic atoms that are predominantly potential scat- 
terers in high-Tc superconductors. Although Zn ions 
are nominally non-magnetic, Tc is strong ly suppressed 
by Zn subst i tution of Cu in the planes l|Hottal Il993t 
llshida et all ll99lD . Therefore, it is reasonable to as- 
sume that Zn ions are non-magnetic unitary scatterers, 
see below. 

We analyze virtual impurity-bound states in a d- 
wave superconductor and, within this framework, ex- 
plore possible implications of the assumption that the 
pairing in cuprates is in the dx2_y2 channel. We model 
cuprates as a 2D d-wave superconductor, based on strong 
anisotropy of electronic transport. Our results, can 
be easily extended for any nontrivial pairing state and 
may be relevant, e.g. for heavy-fermion superconductors 
with impurities . Here we closely follow the references 
llBalatskv efdl Il995l iBuchholtz and Zwicknael Il98lt 
'Salko la et gi.L Il996l ll997HStamr]lll987ir 

Main results of this section are as follows: (i) A 
strongly-scattering scalar impurity is a requirement for 
a localized, virtual or virtually bound state ( or reso- 
nance) to exist in a d-wave superconductor. It is intu- 
itively obvious that any strong enough pair-breaking im- 
purity — magnetic or non-magnetic — will induce such 
a state. Indeed, the low-lying quasiparticle states close 
to the nodes in the energy gap will be strongly influ- 
enced even by a non-magnetic impurity potential, result- 
ing in a virtual bound state in the unitary limit. (ii)This 
should be compared to the fact that, in s-wave super- 
conductors, both magnetic and resonant non-magnetic 
impurities produce bound st ates inside the energy gap 
l|Machida and ShibataL 11972^ . The energy fl' and the 
decay rate fl" of this state are given by 



n = n' + in" 



-Ao 



■Kc/2 



log(8/^c) 



ITT 



2 log(8/7rc) 



(7.1) 



where c — cot(5o. These results are computed assuming 
logarithmic accuracy is sufficient, with log(8/7rc) ^ 1. 
In the unitary limit, defined as 5q — > 7r/2 (c — ^ 0), the 
virtual bound state becomes a resonance at f2 — > with 
0. In the opposite case of weak scattering with 
c < 1, the energy of the virtual bound state formally 
approaches fi' ~ Aq and the state is ill-defined because 
Q." ^ Q! (see Fig. ^| The wave function of the bound 
state is found to decay as a power law: ^'(r) ~ l/r and 
is not normalizable. This is consistent with the virtually 
bound state being not really a bound state. Wave func- 
tion is localized along the directions of the vanishing gap, 
so called nodal directions. 



A. Single potential impurity problem 

Consider the single scalar impurity problem with 

i^int - t^OcLck'a (7.2) 

kk' cr 
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where Uq is the strength of the scalar impurity potential 
at r = 0, resulting in s-wave phase shift Sq. 

The scattering is described by a T-matrix, T{uj), which 
is independent of wavevector. The Green's function in 
the presence of an impurity is 

G(k, k'; uj) = Go(k, L^)4k' + Go(k, uj)f{u;)Go{k' , uj){7.3) 

where Ak — Aocos2(^ is the gap function of dx2_y2- 
symmetry, 



Hirschfeld and Goldenfeldl 


Il99,'4 iHirschfeld ei al. 


1988; Led Il993l 


iPethick and Pinesl. Il986h 


Schmitt-Rink et alV 119861 


Shiballl968HStamnlll987|). it 



is known that T = Tofo + T^t^ for s-wave scattering. 
The T-matrix takes the form 



r(w)ii = l/[c-5ii(w)], 



(7.4) 



where ffiiM = ^-i^ Tr G^") (k, c^)(fo + T3)n. The 
quasi bound states in the single-impurity problem are 
given by the poles of the T-matrix 



c = .gii(f^), 



(7.5) 



which is an implicit equation for impurity resonance fJo 
as a function of c, the strength of impurity scattering. 
Choosing the gap function at the Fermi surface so that 
A(iy9) — Aocos2(^, one finds for particle-hole symmetric 

case gii = go{uj) = luj / /S.{if>Y — lu'^) , where the 

\ / F S 

angular brackets denote averaging over the Fermi surface; 
for simplicity, we take {•)fs = / • dip/2iT ^. For |w| <C 
Ao, one finds 



2uj 



log 



4Ao 



ITT 

y 



(7.6) 



In Fig. 131 we illustrate a solution of the Eq. (|7.5() . 

In principle, the solution of Eq. IjT.Sfl is complex, indi- 
cating a resonant nature of the quasiparticle state, bet- 
ter described as a virtual state. This is easily seen form 
Eq. (|7.1() , which solves Eq. (|7.5() to logarithmic accuracy. 
However, as c — > 0, the resonance can be made arbitrar- 
ily sharp. For c = 0, the virtual state becomes a sharp 
reson ance state bound to the impurity (,Balatskv et al\ . 
Il995j) . As c 1^, ri' and VL" increase without bound so 
that ri"/r2' 1^, and the solution becomes unphysical. 
For c > 1, no solution has been found for fJ. ^ 



^ We assume that the energy gap has Une nodes in three dimensions 
with weak quasiparticle dispersion along the z axis; an extension 
to a general three-dimensional case is straightforward. 

^ The related model of the Anderson impurity in an unconventional 
superconductor has been considered by L. Borkowskii and P. 
Hirschfeld, Phys. Rev. B 46, 9274 (1992). The results found 
here for pure potential scattering require the generalization of 
the Anderson model to include the impurity potential phase shift, 
independent of the Kondo temperature. This aspect of impurity 
scattering has not been addressed previously. 



To properly solve a single impurity state one has to 
retain both components of the T-matrix. Assumptions 
about particle-hole symmetry alone are not sufficient to 
leave the T3 contribution out while computing the den- 
sity of states around the impurity site. Full solution, 
Eq. H7.4|) , leads to the definite sign of the resonance en- 
ergy depending on the sign of c. Indeed, if we assume 
that T3 is zero we find that 



Tn 



.90 (t^) 



(7.7) 



and To has now two poles at c = ±go{uj). This means 
that there would be two solutions to the " poles" of the 
T-matrix for each c, one on positive and one on nega- 
tive frequency. This would clearly contradict the obvious 
particle hole asymmetry introduced by impurity. One 
can not get a symmetric density of states if we have only 
repulsive (C/q > 0) or attractive {Uq < 0) impurity po- 
tential. This argument shows that one is not allowed 
to ignore T3 contributions, because T3 is not a smooth 
function of energy near the pole. 

Now we turn to the physical implications of these vir- 
tual bound states in a d-wave superconductor. Consider 
the most interesting case of unitary impurities in the di- 
lute limit, separated by a distance greater than the co- 
herence length ^. Before averaging over impurities, these 
bound states are nearly localized close to the impurity 
sites (see below) and can substantially modify the local 
characteristics of the superconductor: for example, the 
local density of states, observed in STM and the local 
NMR relaxation rates of atoms close to the impurities. 

Consider a local density of states, defined as 

N{r,uj) ^ --lmgn{r,r-uj + iO+) (7.8) 

TT 

with the total Green's function in the presence of the 
impurity 

G(r,r';c^) = Go{r-v',uj) 

+Gair,u;)f{u;)Go{r',Lu). (7.9) 

We find two terms in the local density of states 

7V(r,c^) = iV(^) + iVi,„p(r,c^) . (7.10) 

The first term originates from the bulk quasiparticles, 
which are described by plane- wave eigenstates with E]^ — 
Va + M, g^'HO,^;) ^j:>i/i^-E^) + vl/{u; + Ey,)], 
where itk and Uk are the standard Bogoliubov factors. 
The bulk density of states is constant in the system with 
N{iLj)/No — uj/Aq, for lj <C Aq. The second term, 

iVu„p(r,c^) = -ilm[Go(r,c^)f (u;)Go(r,c^)]n (7.11) 

TT 

originates from the virtual bound state created at the 
impurity. This impurity state will have a form of a 
cross with long tails extended along the gap nodes see 
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Fig. 0] As an clarifying example, consider the limit of 
unitary scattering for which the resonant state is formed 
at -Eimp,n = ^ 0, see also previous Sec. Because 
ImG(")(r,tj = 0) = -'kN{uj = 0) = 0, only the imagi- 
nary part of the T-matrix contributes to Mmp and the 
bound-state probability density is found to decay as the 
inverse second power of the distance from the impurity, 
mostly localized along the nodes of the gap function, 



M„p(r,cj = 0) = Re[G(")(r,c^ = 0)] 



^cxr^^ (7.12) 



and similarly, but with smaller amplitude, in the vicinity 
of the extrema of the gap function. 



A^i,„p(f,w = 0) cx -fr- 



(7.13) 



In addition to the power law decaying large dis- 
tance asymptotics there is also an additional exponen- 
tially decaying piece that decays with ^((^), the angle- 
dependent coherence length of the superconductor, de- 
fined as ^(1^9) = hvp I\1^{lp)\. Exponentially decaying part 
is important to compare the induced DOS to the observed 
in STM near impurity site although it does not change 
asymptotic behavior at large distances. In practice the 
intensity near impurity is mapped out only within few 
lattice sites. 

Gap nodes lead to the power law decay of the wave 
function along all directions at large distances r >> ^. 
This follows from the power counting of the d-wave 
propagator we estimate: G{Y,bj ^ 0) ^ J (i^fcexp(ik • 
r)G(k,w ^ 0) ~ / kdkexp{ik ■ r)^ - 1/r. The fact 
that the impurity state is virtually bound is reflected 
in the logarithmically divergent normalization. This di- 
vergence should be cut off at an average distance be- 
tween impurities at any finite concentration. More gen- 
erally, for an arbitrary position of the resonance, taking 
into account that only one state has been produced with 
E- 



Q' + we find 



iVi,np(r, w) 



\v{r - ri)p 



(7.14) 



where we have introduced the sum over different impuri- 
ties, located at r^, and u(r — r^), ^(r — r^) are the eigen- 
function of the Bogoliubov-de Gennes equation at the 
impurity level. 

The local effects of impurities are best revealed by local 
probes. NMR experiments on Cu in Zn-doped cuprates 
are quite useful in this regard. From Eq. H7.12(l and be- 
low, one concludes immediately that the local NMR sig- 
nal would show two distinct relaxation rates (or even the 
hierarchy of rates): one coming from the Cu sites, far 
away from the impurities, and another from the sites, 
close to the impurities. The Cu sites near the impuri- 
ties will be sensitive to the higher local density of states 



and will have a higher relaxation rate at low temper- 
atures. At finite impurity concentration (~ 2%), the 
volume-averaged density of states will have a finite limit 
at w ^ 0, as follows from Eq. (|7.12l) . The relaxation 
rates of Cu atoms close to and away from an impu- 
rity will, therefore, have the same temperature depen- 
dence (TiT)~^ = const, but will be of a different mag- 
nitude. Precisely this be havior has been obse rved exper- 
imentally: Ishida et al. ijishida et all Il99lj) have mea- 
sured two NMR relaxation rates for Cu in Zn-doped 
YBa2Cu307_5. The second NMR signal with higher re- 
laxation was inferred arising from the near-impurity Cu 
sites. Alloul and collaborators have pointed out that 
the NMR signal coming from the sites close to impurties 
shows a distribution of relaxation times and reflect local 
electroni c and magnetic disto rtions produced by impuri- 
ties, see llBobroff et qiJ . l2n0lD and references therein. 

More direct evidence for the impurity induced reso- 
nances in high-Tc is coming from STM experiments. Lo- 
cal variations of the density of states can be probed using 
scanning-tunneling microscopy. These experiments were 
crucial in establishing the existence of the impurity in- 
duced resonances in cuprates and their anisotropic nature 
l)Hudson et a;J . ll999tlPan et a/.l.l2000b|) . see section lXIIll 

We would like to contrast our picture of the dilute 
limit of strongly scattering centers to the usual approach 
of averaging over impurities at finite concentration. If 
one considers averaging over impurities, two NMR relax- 
ation rates, arising from unequivalent sites, cannot be 
resolved; similarly local inhomogeneous aspect of the lo- 
calized states will be lost after averaging over impurity 
positions. 

For practical purposes the distinction between the a 
true bound states and continuum in our case is not well 
defined, as it is in s-wave superconductors. Any finite 
temperature will produce a finite lifetime for these bound 
states, and they will be hybridized with the continuum 
of low-energy quasiparticles as they are not separated by 
a well defined gap. 



B. Single magnetic impurity problem 

The similar analysis for the magnetic impurity is more 
involved. For a quantum spin s = 1/2 in d-wave super- 
conductor one needs to address the Kondo effect. It is 
discussed in more details in Sec. IXll For a classical spin 
S ^ I the analysis within t he mean field i s simil ar to the 
one in the previous section fealk ola et gIlll997D . 

The main effect is that exchange coupling between the 
local spin S and electron spin leads to the renormaliza- 
tion of the effective scattering potential. Namely for elec- 
trons of two spin polarization the net impurity potential 
is Uq ± J, where Uq is the potential scattering strength 
and J is the exchange coupling to impurity spin, see 
Sec. IXIIII There are two virtual bound states, one for 
each electron spin orientation. STM data on Ni-doped 
Bi2212 are fit well using this simple formula. This mean 
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field approach does not address the dynamics of the large 
spin S. More analysis is required to address this problem. 



C. Self-consistent gap solution near impurity 

Impurity scattering will produce local modifications of 
the order parameter. We already addressed some of the 
effects in s-wave superconductors in Sec. lIXI Here we will 
discuss the self-consistently determined gap in d-wave su- 
perconductors. 

To address these effects at small distances one would 
need to use an numerically determined exact spectra near 
impurity and solve self-consistent gap equation. We de- 
fine: 



A{i,i + S) = 



■^[un{i + (5X(«) + <(i)w„(i + S)] 



(7.15) 



Numerical solution of this problem was presented in 
jFranz et a[.^ 19 96; Salkola et ai, 1997: Tsuchiura et all . 
I2OOO; Zh u e"<: ayj , i20 00a). The main result is that impurity 
scattering suppresses the gap magnitude. Suppression is 
strongest on the impurity site and quickly gap recovers a 
bulk value, although there are always oscillating tails at 
far distance due to 2kp oscillations, Fig.|Sl 

In practice, the difference between self-consistent and 
non self-consistent solution is not important at distance 
beyond few lattice sites away from impurity. Local gap 
suppression is clearly seen in STM data, see Sec. IXIIII 



D. Spin-orbit scattering impurities 

Spin-orbit coupling in impurity scattering in supercon- 
ductors is the least discussed among all other kinds of 
impurity scattering. Standard spin orbit scattering is of 
the form: 

Hsoimp = ^ AsoCk „i?Q/3 • (k X k')ck'/3 (7.16) 

k,k' 

where Xso is the strength of Spin Orbit scattering. This 
kind of SO scattering would be present even for non- 
magnetic impurities, it will be a pairbreaker and it will 
produce the the quasi-bound states inside the gap, al- 
though the detailed calculation has not been done to our 
knowledge. 

Another kind of SO scattering impurities scattering 
in d-wave from the magnetic impurity can be consid- 
ered as well. It was initially motivated by e xperi- 
ments on Ni doped Bi2212 (Movshovich et ai, Il998t 
iNeils and HarlingenL 12002^) . In this approach impu- 
rity spin is coupled to the orbita l motion of the con- 



duction electrons (Balatskvl Il998t iBarash et al. I. Il997t 
iGraf el J . 1200(1 lOrimaldl 1 19991) : 



where 7so is the strength of coupling and S is the im- 
purity spin. This term is the SO coupling Hso,imp — 
7soL • S written in second quantized notation. Pre- 
dominantly in plane motion of electrons as is the case 
in Bi2212 will couple Lz to Sz, Hiz = ihd^ is the angular 
momentum operator Lz with respect to the impurity site. 
The net effect of this term is twofold. This scattering 
term is a pairbreaker, locally gap is suppressed and reso- 
nance is formed. However the more interesting and non- 
trivial is the distortion of the d-wave order parameter in 
the vicinity of impurity, which results from the nontrivial 
orbital structure of the d-wave order. The d-wave state is 
a linear combination of the state with / = 2 and I = —2, 
A(0) = Aocos(2(/)) cx exp(2z0) -I- exp(— 2i0) ^ — y"^. 
The orbital angular momentum components will be af- 
fected differently as a result of scattering. In the first or- 
der perturbation theory in Hsoimp one generates the cor- 
rection to the order parameter A = jAo7so sin(20) ~ 
xy. There is a finite amplitude for incoming d-wave pair 
\in) cx jx^ — y^) to scatter into \out) cx i\xy) channel: 

\out) — i75oAoL • Sjin) = 1/17^0 Aq sin(20) . (7.18) 

Therefore as is the case for SO scattering, sometimes im- 
purity scattering can produce nontrivial distortions of 
the initial order parameter, aside from trivial suppres- 
sion. For more details, see ^Balatskvl IT998t iGraf all 
2000; Zhu and Balatskv, 2002). Similarly, magnetic field, 
that acts similarly to the Sz term in Eq. (|7.18l) . not 
only suppresses the d-wave order parameter but also pro- 
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HsO,^mp = XI ^So4,<tS ' (k X k')ck'a (7.17) 
k,k' 



E. Effect of doppler shift and magnetic field 

Here we are focusing on the orbital effect of magnetic 
field and thus the problem is closely related to the effect 
of Doppler shift on impurity resonance. 

In the presence of a superflow with velocity vs propa- 
gators are modified: G'(k, lu) G(k, w — k • V5) for a pla- 
nar wave state at vector k and similar change for F func- 
tion. Hence the rest of the calculation for the impurity 
state goes through as before. Since only local propaga- 
tors enter into solution for impurity resonance Eq. (|7.1|l . 
the modifications will arise as a result of changes in den- 
sity of states due to Doppler shift. 

Interesting effect of the superflow produced by the 
screening currents on th e impurity state was studie d 
by Samokhin and Walker ijSamokhin and Walkeill200l|) . 
They pointed out that Doppler shift will result in the 
broadening of the impurity induced resonance. This is 
a consequence of the local scattering nature of impurity 
that means summing over all momenta to obtain local 
Green's function Go(w). One would need to compare the 
typical value of the Doppler shift vskp and the energy of 
the resonance f2'. In the case when Doppler shift is small 
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effect is negligible. In the opposite limit of vskf ^ fl' 
superflow produces broadening of the resonance but not 
the energy shift of the resonance. 

F. Sensitivity of impurity state to details of band structure 

In the above discussion we were using the single band 
model with particle-hole symmetric structure as a sim- 
plest example to prove the existence of the impurity in- 
duced resonance. The effect of asymmetri c band abou t 
the gap midpoint was considered by Joynt ijJovntl Il997|) , 
by assuming a constant DOS with different energy ranges 
outside the gap edge. To make a comparison with the real 
experimental data on impurity resonances, see Sec lXIlH 
one has to understand the details of realistic band struc- 
ture. Microscopically, relevant bands in Cu-0 planes are 
Cu dx2-y2 and O px,y bands. In the above analysis we 
have assumed that upon the reduction of the complicated 
band structure of high- Tc to a single band model, one can 
still describe nonmagnetic impurity by a single parame- 
ter, i.e. the on-site potential Uq. We do not have a proof 
for this and assume that the major physical effects, such 
as that the impurity induced resonance will be properly 
captured in a simplified model. 

Within the framework of this simplified one band 
model one can still investigate some effects beyond the 
simplest assumptions. One of the most obvious ques- 
tions is the position of the impurity-induced resonance 
with respect to the Fermi energy. We find the reso- 
nance depends on the sign of the impurity potential, the 
electron occupation, and the band structure. We have 
performed numerical exact diagonalization for the t-t'- 
V model with nearest-neighbor hopping next nearest 
neighbor hopping t' , and a negative V , that describes 
the nearest neighbor attraction to produce effectively a 
d-wave pairing. The single particle energy dispersion for 
the normal state is given by: 

^fc — — 2t(cos kx -\- cos ky) — At' cos k^ cos ky — fi , (7.19) 

with fi the chemical potential. Impurity was modeled as 
an on-site potential Uq. We have looked at three possibil- 
ities: (i) t — l,t' — 0, fj, = (the filling factor n — 1.0), 
with band particle- hole symmetry present. Fig. El (ii) 
t = l,t' = -0.2,^ = -0.784(n = 0.84), with no band 
particle-hole symmetry, Fig.[7| (iii) t = l,t' = —0.3,fi — 
— 1.0(n — 0.85), again with band particle-hole symmetry 
absent, Fig.|Sl Here we talk about the band particle-hole 
symmetry because as long as an impurity is introduced, 
the local particle-hole symmetry is always broken. 

As shown in these figures, for the cases (i) and (ii), the 
band DOS has two coherent peaks. Also for the case (ii), 
the DOS is asymmetric with respect to the zero energy 
point. In these two cases, a repulsive potential Uo > 
leads to a negative energy impurity state fig < 0. This 
position is manifested as a resonance peak appearing be- 
low the Fermi energy in the LDOS directly impurity site 
but appearing above the Fermi energy in the LDOS at its 



four nearest neighbors. An attractive impurity potential 
?7o < induces a positive energy impurity state fig > 0, 
as reflected as a resonance peak above the Fermi energy 
in the LDOS directly on the impurity site but below the 
Fermi energy in the LDOS at its nearest neighbors. 

For the case (iii) , in addition to the two coherent peaks, 
there are also two van Hove singularity peaks (the one 
on the negative energy side being more pronounced and 
the other on the positive energy side being faint). Then 
for a repulsive impurity, the on-site resonance peak does 
shift from the negative energy side across the zero energy. 
This phenomenon is absent for the cases of (i) and (ii). 
For Uo < 0, the result is similar to the cases (i) and (ii). 
We point that in the STM data, the van Hove singular- 
ity peaks are absent. Here we have chosen for (ii) and 
(iii) only in the optimal doping regime. In other doping 
regime, all possibilities uncovered above could appear. 
More detailed analysis, especially realistic band struc- 
ture calculations will allow us better address the details 
of impurity states in high- Tc materials. 

As for the sign of the impurity potential from the Zn 
and Ni atoms in cuprates, it is still an unsettled issue. It 
is believed that these atoms substitute Cu in the Cu-0 
plane, and do not change the hole doping. In case of Zn 
that has Sd^^As^ electrons, Zn++ is in d^'^ configuration. 
The third ionization energy should be a rough measure of 
the energy level for d-orbital of Cu and Zn/Ni ions though 
the electrons from the d-orbital of Cu form a band. By 
comparing the energies of Cu atom Ec^++ = — 36.83ey, 
and Zn atom Ezn++ = — 39.722el^, we estimate energy 
to be on the order of Uq ~ — 2.89el^. Therefore, Zn atom 
plays the role of a strong attractive potential. 

In case of Ni doping, Ni++ has a 3d^ configuration 
and a spin S = 1 ground state is formed. Therefore, Ni 
impurity will produce both potential scattering Uq and 
magnetic scattering J. We can estimate the energy Uq 
again by taking the difference between atomic energies 
for Ecu++ = -36.83ey and £^1++ = -35.17ey. We es- 
timate Uq ~ 1.66ey for Ni. It provides a weaker repulsive 
potential for the non-magnetic scattering part. 



VIII. SINGLE IMPURITY BOUND STATE IN A 
PSEUDOGAP STATE OF TWO-DIMENSIONAL METALS 

In the high-tempe ra ture cu prates, many exper- 
iments llLoram et ali l2000t iNorman et~al\ . Il998t 
iRenner et all 1 1998(1 show that the electronic density 
of states near the Fermi surface is suppressed within 
the range of Apc above the superconducting phase 
transition temperature Tc but below a characteristic 
temperature T*. So far, the mechanism for the PC 
state is still hotly debated. F or a review see (jTimuskl 
l2003t iTimusk and StattL Il999|) . The typical competing 
scenarios for this anomalous phenomenon, including 
mainly the pre-formed pa i r wit h phase-fluctuation 
model (|Emerv and Kivelsonl Il995(l. the Bose- Ein stein 
condensation of Cooper pairs ijChen et all . Il998(l . the 
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time-re versal-symmetr y- breaking circulating current 
model (IVarmal Il999^■ and t he rf- density-wave model 
(DDW) l|Chakravartv e.t all 120011). which is typical 
kind of st aggered flux state ('Affleck and Marstonl llQSSt 
iHsu et a?i . .l991; Marston and Affleck. 1989). In the first 
scenario, the normal state contains preformed Cooper 
pairs, and the phase fluctuation of the pairing field 
destroys the long range order, that is, the superconduc- 
tivity. Since the pairing field has at the onset a d-wave 
symmetry in the momentum space for the relative 
motion degrees of freedom, a d-wave-like pesudogap 
follows naturally. In the second category of scenarios, 
it is speculated that the PG state comes as a result of 
a new order not inheriting from the superconducting 
pairing. 

In this Section, we are not going to discuss the origin 
for the PG. Instead we are particularly interested in the 
consequences of the electronic property around a single 
impurity in the PG state. One can consider the temper- 
ature evolution of the impurity state as the temperature 
increases and eventually becomes larger than Tc- Then 
there are two possibilities for the evolution of impurity 
resonance at T > Tc- a) the impurity resonance grad- 
ually broadens until the superconducting gap vanishes, 
at which point the impurity resonance totally disappears 
and b) the resonance gets broader but survives above 
Tc- Which of the possibilities is realized depends on 
the normal state ph ase the superconductor evolves into . 
It has been argued ijKrasnov et alV l2000t iLoram et alV 
[2000) that in the underdoped regime the superconduct- 
ing gap opens up in addition to the PG present well above 
Tc- Hence, we find that the impurity resonance survives 
above Tc in the PG state of high-Tc materials. The po- 
sition and the width of the resonance are determined by 
the impurity scattering strength and PG scale. In the 
absence of PG above Tc the impurity state disappears. 

Specifically we calculate the resonant state generated 
by the substitution of one Cu atom with a Zn atom us- 
ing the T-matrix approach. We rely on the fact that the 
density of states (DOS) is depleted at the Fermi energy 
in the PG regime. We argue that the mere fact that the 
DOS is depleted at the Fermi energy is sufficient to pro- 
duce a resonance near the nonmagnetic impurity, such as 
Zn. Before we consider the impurity effect within spe- 
cific scenarios, we give a general analysis, which is valid 
in the PG state with no superconducting phase or ampli- 
tude fluctuations above Tc, as long as there are interac- 
tions that lead to the PG state. The approach we take is 
similar to the previous analysis of the nonmagnetic impu - 
rity in the superconducting state ijBalatskv et a/.Lfl995|) . 
See also Fig.|^ The states generated by the impurity are 
still given by the poles of the T matrix: 



.1) 



This is an implicit equation for Q as a, function of C/q- 
This solution can be complex, indicating the resonant 
nature of the virtual state. To solve this equation, we 



need to know the form of nonperturbed Green's function 
on the impurity site, go{E). To do so, we split go into 
its imaginary and real part go = g'^ + ig^ with Gq{uj) = 
— 7rA'o(w), where iVo(a;) is the density of states. 

Mea surements on t he ele ctronic specific heat by Loram 
et al. l|Loram et a.l\ . l2000fl show that the normal state 
pseudo gap opens abruptly in the underdoped region be- 
low a hole doping equal to pcrit ^ 0.19 holes/Cu02. In- 
spired by these data, we will assume that around the 
pseudogap region, states are partly depleted and the den- 
sity of states is linear, that is Nq{ijj) ~ A^oI^I/Apg for 
\lo\ < ApG and No{u}) = Nq for Apc < \uj\ < W/2 with 
W the bandwidth. This density of states is depicted in 
Fig. [rUT al . As it is obvious from the solution of Eq. (|8.1|l , 
the precise position and the width of the resonance will 
depend on the specific form of the PG. We will use this 
linearly vanishing PG DOS. Results for other forms of 
N{u!) like a fully gapped DOS or a DOS with a quadratic 
dependent gap, can be obtained in the same way and lead 
essentially to similar expressions. 

Using this DOS for and invoking the Kramer-Kronig 
relation ifMahan. 2000) 
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with P the Cauchy's principle value, one can obtain the 
real part g'^ as 
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(8.3) 



This function is plotted in Fig. llOf b) together with I/C/q- 
If 2UqNo > 1, one can see from this figure that equation 
(|8.1|) has four solutions. Since the width of a resonance 
state is proportional to \Q\, the only state with sharp 
width is the solution with \Q\ close to zero and we will 
only consider this solution. After expansion in lu of equa- 
tion (|8.3I) we arrive at an expression for this solution f2 
of Eq. l|0) : 
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(8.4) 



^ We argue that the appearance of the intragap impurity 
state is a robust feature of any depleted DOS around the 
Fermi surface. We also considered the model DOS with 
No{iii) = No [a + {1 — a)iii^ / A^q] which leads essentially to simi- 
lar results as a function of the impurity strength with a resonant 
state at U = -Apa(i + iT^aNoUo) / (iNoUoil - a - Apc /W)) 
^ -Apa{l + i-!TaNoUo)/{4:NoUo{l - a)) when Apc/W is 
small. But also a fully gapped DOS equal to N{u)) = Nq for 
|aj| g [ApG, W/2] and zero otherwise gives rise to a comparable 
expression with a resonant state at H = —Apq/{2UoNo). 
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This equation can be solved exactly in terms of Lam- 
bertW functions, which, to logarithmic accuracy with 
ln|2;7oiVo| > 1, gives: ^ 



n = 



n' + m" 



1 



2t/oA^oln|2C/oiVo| 

1 in sign([/o) 



1 - 



ln|2i7oiVo| 2\n\2UoNo\ 



.5) 



where il' is the energy position and fl" the decay rate. 

Using formula (|8.5|l . and taking A'o = 1 state/eV, 
Apc ~ 300K ~ 30meV and the scattering potential 
Uq ~ ±2eV, we estimat e 1] ^ zfc2meV ~ ifc 20K as was 
found by Loram et al. l)Loram et al\ . l200Cl|) . This en- 
ergy is close to the Zn resonanc e energy u>n = — 1%K, 
seen in the superconducting state ijPan et adl2000bj) . By 
combin ing these results w ith the band-structure argu- 
ments ijMartin et a;J.l2002() . we come to conclusion that 
the Zn impurity in Bi2212 is strongly attractive, with 
[/q ~ — 2eV. This result, as we will now see, may be 
modified due to the particle-hole asymmetry characteris- 
tic of doped cuprates. 

In the absence of particle-hole symmetry, a similar cal- 
culation can be done. The simplest way to introduce the 
asymmetry is by making the upper and lower cutoffs in 
the DOS unequal. This situation corresponds to a chemi- 
cal potential /i, located away from the center of the band. 
Keeping the DOS otherwise unchanged, with the pseu- 
dogap centered at the chemical potential, results only 
in the following change in the first logarithmic term of 
Eq. 
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Neglecting the frequency lo relative to chemical potential 
/J, and assuming that ^ is small relative to the bandwidth, 
we obtain that the results for the asymmetric case can 
be obtained from the symmetric ones by the substitution 
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The effect of the asymmetry term can be estimated for 
the superconducting cuprates. For 20% hole doping, 

^l (l/5)W^/2 = -W/IQ. Hence, the modified value 

for the Zn impurity strength in Bi2212 can be obtained 



from the symmetric result, l/U* = l/Uo+ANofi/W. The 
new value is C/* ~ — 1 eV, which is a strongly attractive 
potential, as is expected from the band structure argu- 
ments. 

The solution of the impurity state deep in the super- 
conducting regime involves two aspects: the energy po- 
sition and the width of the resonance and secondly, the 
real space shape of the impurity state. We have dis- 
cussed the energy of the impurity state above. Great 
advantage of the on-site impurity solution for the local 
potential Uq is that only on-site propagator goi^) enters 
into calculation. Hence the knowledge of the DOS was 
sufficient to calculate the impurity state. On the other 
hand, to calculate the real space image of impurity in- 
duced resonance, one would require more detailed knowl- 
edge of the Green's functions in the PG regime. Quite 
generally, one would expect for a c?-wave-like PG with 
nearly nodal points along the (±7r/2, ±7r/2) directions, 
that the impurity resonance in the pseudogap regime 
would be four-fold symmetric, si milar to superconducting 
solutions l)Balatskv et all Il995f) . This calculation would 
require a specific model for the PG state, some of which 
are considered below. 

Impurity state in a DDW system. The model Hamil- 
tonian for a clean DDW system is written as: 



Here {cia) is the creation and annihilation of an elec- 
tron with spin projection a at the ith site. The DDW 
order parameter Wij has the value Wi^i±x = W±x = ^ 
and Wi^i±y — W±y = — while is zero otherwise. No- 
tice that the prefactor i = a/— T before Wij . It indicates 
that the DDW state breaks the time reversal symmetry. 
The quantity is the chemical potential. Without loss of 
generality, we have assumed that the impurity is located 
at the origin. In the momentum space, the Hamiltonian 
is given by: 



k,(T 



(8.9) 

Here cj,^ {cka) are the creation and (annihilation) oper- 
ators of an electron at the wave vetor k and with spin 
projection a. The single particle energy dispersion for 
the normal state is given by Eq. H7.19|l . The DDW gap 
is given by: 



The exact solution in terms of a LambertW function, Lw{—1, x), 
is f2 = — Ap(3Sgn(C/o) exp{Lui(— 1, — sgn(l/o) exp(j7r/2 — 
l)/{2NoUo)) -I- 1 - iTT/2}, where Lw{x) is such that 
Lw(x) exp[Lw{x)] = x. 

The simplest model for thermal broadening is to assign the tem- 
perature dependent width: Thermal broadening at high temper- 
atures T > Tc substantially broadens the impurity resonance 
peak Q"(T) = ^/{n"{T = 0))2 -|-T2. 



Wo 



(cos kx 



iky) . 



(8.10) 



Here the wave vector k^ and ky are defined in the first 
Brillouin zone. 

In view of the fact that the DDW state breaks the 
translational symmetry with lattice constant but con- 
serves that by \/2a along the diagonals of the square 
lattice, it is convenient to rewrite the Hamiltonian in the 
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reduced Brillouin zone with its area one half of the orig- 
inal by introducing a two-component electron operator 



with Q — (tt, tt) 



Ho = 



E 

kerBZ,a 



ka 



-2iWk £.k+Q 



*fc. , (8.11) 



where rBZ represents the reduced Brillouin zone. 

Accordingly, we can introduce the following Green 
functions in the clean limit: 



( ^1? 



cio) 

C(0) 
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with 



g[l\k;r) 
{k-r) 



-(rr[cfc.(r)4(0)]) , 
-{Tr[cka{r)cl^Q,M\) 
-{Tr[ck+QAr)clM]) 



C(0) 
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-{Tr[ck+Q,a{T] 



k+Q,a 



(8.13a) 
(8.13b) 
(8.13c) 
(0)]> .(8.13d) 



Here the factor Tr is a r-ordering operator as usual, and 
Cfecr(r) = e^^ Ckae~^'^ is the operator in the Heisenberg 
representation. Given the Hamiltonian Eq. H8.11|l . with 
aid of the equation of motion for the field operator Cka [t) 
and c\^{t), and by performing a Fourier transform with 
respect to t, 



Q{k] t) = kBT'^Q{k] iujn)e'' 



with ujn = (2n 

g^''\k-iLo^) 



ViTTkuT, one can arrive at: 



2iWk 



-2iWk 



(8.14) 



(8.15) 



As will be obvious immediately, we also need to calcu- 
late the Green function in real space, which through the 
Fourier transform is given by 



The existence of the resonant states will directly manifest 
in the local density of states: 



Ni(uj) = Im CJ(i, i]LU + i5) 



(8.19) 



The results are displayed in Figs. El and El As 
shown in these figures, the electronic excitation spectrum 
around the impurity in the DDW state is very sensitive to 
the parameter values, which control the band structure. 
For t' = and at the half filling (/i = 0), the electron 
density of states in the clean limit is vanishingly small 
around the Fermi energy, which leads to resonance states 
near the Fermi energy. With t' = but the system doped 
away from the half filling, the resonant peak in the LDOS 
is shifted away from the Fermi energy. This is because 
the energy at which the band DOS vanishes no longer 
coincides with the Fermi energy. For a set of more real- 
istic parameter values, the density of states in the clean 
limit shows negligible reduction at low energies, which 
makes the local density of states near the impurity not 
exhibit any signature of a resonance state. T hese results 
were indepen dently obtain ed by Zhu et al. llZhu et all 
120011) . Wang (|mni|200l, and Morr llMoril 120021) . The 
quasiparticle states in the DDW state with a finite con- 
centration of non- magnetic impurities was investigated 
by Ghosal and Kee ijGhosal and Ke(j.l2004|) . 

Phase- fluctuation scenario: We now devote to a discus- 
sion on the impurity state in a pha s e-fluc tuating pairing 
field. For more details, see IWand iooi The effective 
mean field Hamiltonian in a two-dimensional square lat- 
tice for a d-wave superconductor can be written as: 



't, 



-A* 



{~tj, - ^%) 



, (8.20) 



k-Ri 



\G^^'' (k- iuj ) -\- G^'''' (k- iuj ) 

ivii \ , n) -r ^22 ^n) fluctuating. For d-wave superconductivity, one can write 
^^"^^^ the pairing field as: 

-'Q-^^^f2HA:;za;„) + e'Q-^-g^;Hfc;zc.„)] , 

(8.16) = -^e'^^= = Ay-e*'^'^- , (8.21) 



7(0), 



where = (c|| , Cii ) is the two-component spinor op- 
erator in the Nambu space, the other notations are the 
standard except the pairing field A^j is fluctuating. In 
the phase fluctuation scenario, the amplitude of the su- 
perconducting order parameter is fixed while its phase 



where are the lattice vectors and Rf 



Ri — R, 



Specially, from Eqs. H8.15|l and (|8.16|) . one can obtain: 



7(0) 



1 

TV 



y 



2iw„ — ik+Q — ik 



Cfe)(iw„ - iu+q) 

(8.17) 

In the presence of a single non-magnetic impurity in 
the DDW state, the T-matrix analysis just shows again 
that the resonance state is determined by the poles of 
T{iLOn —fLO + i0+), i.e., the following equation 



with riij = 1 (— 1) for x (y) direction bonds and the phase 
Lpij = [ipi + ipj)/2. The spatial variation of the phase will 
give rise to a superfluid flow associated with the Cooper 
pairs. By performing a gauge transformation. 



i<PiO-3/2^j, 



.22) 



with (T3 the third component of the Pauli matrix, one can 
remove the phase factor on the pairing field such that 



g(")(0,0;a;-f zO+) 



1 



(8.18) 
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E*l 



-A*- 



-A,,- 



, (8.23) 
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where iij = tije~'^^'^'~'^^^^'^ . In type-II superconductors, 
the length scale of the phase variation (the London pene- 
tration depth) is much larger than the Fermi wavelength. 
Hence one can specify the phase for the Cooper pair 
ipi = 2qs • Ri, where is the average momentum per 
electron in the superfluid state. This ansatz leads to the 
Green function in the clean limit: 

(8.24) 

where At = (cos fc^ — cos fcj, ) and the energy dispersion 
is still given by Eq. H7.19|) . 

In the presence of a non-magnetic single-site impurity 
at site i = (0, 0) in the 2D system, the Green function 
for the impurity system can be obtained, within the T- 
matrix approximation, 

Q{i,j;qs;iuJn) = ^"^"^ (i, j; q^; icj„) + (i, 0; g^; i^n) 
xr(g,;ic^„)gW(0,j;g,;zL^„) , (8.25) 

with the T-matrix given by 

f-i(iw„;q,) =r3/[/o-a^"^(0,0;g,;iw„) . (8.26) 

With a fixed qs, the local density of states (LDOS) at 
site i is given by: 

N{i;qs;uj) ^ ~{2/TT)lmgii{i,i;qs;cu + iO+) . (8.27) 

For the fluctuating phase, one needs to take 
average over q^. If the fluctuation is ther- 
mal, the statistical distribution is Gaussian, as 
gove rned by the Koster l itz-Thouless (KT) the- 
ory l)Kosterlitz and Thoulessl Il973l Il974t ISheehv et all . 
120011) . that is piq,) = e-«'/2"''/Eg, e-'?'/2"„^ .^l^ere 

riy = exp[-y/aTc/{T -Tc)] within the KT theory. In 
the continuum limit, \/J(ff) = ^/n^- The averaged LDOS 
is caculated as N{i;uj) = {N {i; qs] uj)) . 

The results for various values of n„ are shown in 
Fig. 1131 For small n^, the resonance peak is very sharp 
and similar to that in the superconducting state {riy — 0). 
When riy is increased, the peak is broadened with de- 
creasing peak height. Finally, the spectrum at low ener- 
gies becomes featureless. 

Finally, in the phase fluctuation scenario of the PG 
state, the electron excitation spectrum around the im- 
purity is very sensitive to how far the temperature is 
away from the superconducting transition temperature 
Tc- However, in the normal-state ordering scenario, the 
resonance states are not sensitive to the temperature up 
to the PG critical temperature. Another piece of differ- 
ence of the electronic states around the impurity between 
two scenarios is the energy position of the resonance state 
in the phase fluctuation state is not sensitive to the dop- 
ing while that in the state with a normal ordering shifts 
with the doping. More generally, if the superconducting 
fluctuations are present, then an additional satellite peak 



should appear at the opposite bias due to the particle- 
hole nature of the Bogoliubov quasiparticles. The rela- 
tive magnitude of the particle and the hole parts of the 
impurity spectrum can be used to determine the extent to 
which the PG is governed by the superconducting fluctu- 
ations. In the case of fully non-superconducting PG (e.g., 
the DDW state), there should be no observable counter- 
part state^_Combined with other experimen tal propos- 
als l|Janko et a/.l. Il999l: iMartin and Balatskvl 120001. the 
impurity state can help to better understand the myste- 
rious PG state. 



IX. QUANTUM PHASE TRANSITION IN S-WAVE 
SUPERCONDUCTOR WITH MAGNETIC IMPURITY 

A. Introduction 

In this Chapter we revisit the problem of a localized 
classical magnetic moment in a superconductor. A re- 
markable aspect of this interaction we will focus on here 
is the first- order zero temperature transition which takes 
place in an s-wave superconductor as a function of the 
magnetic moment, JqS, where S is the local impurity 
spin and Jq is the exchange coupling. In this transition, 
the spin quantum number s of the electronic ground state 
of the superconductor \'^o) changes from zero for a sub- 
critical moment Jq < Jcrit to 1/2 for Jq > Jcrit- The 
total spin becomes S' ± 1/2 depending on the sign of the 
exchange interaction Jp . The first to p oint out the phase 
transition was Sakurai ljSakurail Il970^ it corresponds to 
a level crossing between two ground states as a function 
of the exchange coupling Jq. In a singlet superconductor 
level crossing occurs between states where the impurity 
spin is either unscreened or partially screened. The states 
have different spin quantum numbers and level crossings 
are generally allowed. These quantum phase transitions 
are of the first order and hence do not have divergent 
time or length scale associated with them. 

We address the above problem at zero temperature by 
using the mean-field approximation within the T-matrix 
formulation and the self-consistent approach, which takes 
into account a local gap-function relaxation. Local 
Coulomb interaction U which breaks particle-hole sym- 
metry and leads to an asymmetric spectral density for 
the impurity-induced quasiparticle states. Figure 1141 il- 
lustrates the local effect of a magnetic moment on the 
low-energy spectral density in an s-wave superconduc- 
tor. Since we limit our considerations to a classical spin, 
S" 3> 1, the impurity moment cannot be screened com- 
pletely by the quasiparticles. We show that the gross 
features of the impurity-induced quasiparticle states in s- 
and c?-wave superconductors can be qualitatively under- 
stood within the non-self-consistent T-matrix formalism. 

The transition is not unique to the classical spin. Sim- 
ilar effect is found in a Kondo model of a quantum spin, 
see Sec.lXl 
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B. Quantum phase transition as a level crossing 

The physical picture of the quantum transition fol- 
lows from the behavior of the impurity-induced bound 
state. This transition is a consequence of the instability 
of the spin unpolarized ground state, because, for a large 
enough Jq , the energy of the impurity-induced quasipar- 
ticle excitation would fall below the chemical potential. 

As we have discussed it in the context of Yu-Shiba- 
Rusinov solution for a classical spin, see Sec. I VII impurity 
state has energy that is always below the gap threshold : 



1 ~ inJoSNo)^ 
1 + {nJoSNoy 



(9.1) 



and has particle (u_i) and hole (w_i) amplitudes at posi- 
tive and negative energies. One can see the level crossing 
and change of the ground state already from this result. 
Ignoring the self-consistency, and using Eq. H9.HI , we find 
that the singularity occurs at 



^0 = Jcrit = I/t^NqS 



(9.2) 



For a weak coupling Jo < Jcrit , the ground state of the 
superconductor is a paired state of time-reversed single- 
particle states in the presence of the impurity scattering, 
with ground state wave function being the BCS like : 



1* 



n^.K + wnViV'lJIo) = |*o) (9.3) 



Here, since the translational symmetry is broken by im- 
purity, and we consider the eigen-states of the scattering 
problem in the presence of impurity, states are labeled 
by a discrete index n = 1, ...oo corresponding to a dis- 
crete scattering states that are the basis of the Bogoli- 
ubov Hamiltonian with impurity. The n = 1 state would 
correspond to an impurity bound state, localized on im- 
purity site. Index —n correspond to a time reversal state, 
i.e. the localized state with opposite spin. The first ex- 
cited state at Jq < Jcrit would correspond to a single 
quasiparticle state where one excitation is present. The 
lowest excited state corresponds to a intragap Yu-Shiba- 
Rusinov state at energy Qq, see Fig. [T5l 
The excited state can be written as: 

|*-l)j0<Je„. ^7lll*0) = |<I'-l) 

I'i'-i) - i'-i n + ^n^i^Lm (9.4) 

ri>l 

with standard quasiparticle definitions of 71 = uiipi 



with ui + 



uitp[ — Vitp. 



wiV'li + viipi, etc, 



for aU quasiparticles: e.g. 7±i n„>o[""+^"V'i^l„]) = 0. 
^ This state is a true spin singlet (^'o|SeZectron|^'o) — 0. 
To avoid confusion with impurity spin S we explicitly 
indicate that ^eiectrori is the net spin of conduction elec- 
trons. Hence if JoKJa^t = l*o) is a ground state, the 
total spin of electrons is zero, and only the spin of im- 
purity counts. The first excited state at energy fio has a 
spin 1/2 quasiparticle in it: -l\Sll^^^^^J^ -i) = -1/2. 

Upon increasing coupling constant Jq one reaches the 
regime where these two states cross, Fig. At that 
point the excited and ground states changes the roles. 



l^o)jo>Jc„. = = l*-l> 

|*-l)jo>J.„* = 1*0) 



(9.6) 



Another way to see this quantum phase transition is 
to examine energy levels as a function of Jo / Jcrit ■ For 
variational wavefunctions |5'o,-i) we define the respec- 
tive energies as expectation vales of the Hamiltonian: 



Eo,-l{Jo/Jcr^t) = (%_l|iJ|%-l) 



(9.7) 



Energy of the first excitation, the impurity bound state 
is then simply 



^o{Jo/ Jcrit) ~ E-l — Eq, Jq < Jcrit 
^o{Jo/ Jcrit) = Eq — ^0 > Jcrit 



(9.8) 



There are several implications of this result. Firstly, 
the ground state of superconductor with a magnetic im- 
purity in the strong coupling limit is a non-BCS state. 
The ground state is a pair condensate with a single un- 
paired quasiparticle. Similar result was observed for a 
Kondo screening in superconductor l)Sakai et fl/J . ll993[) . 
One can easily understand the result by going to a strong 
couphng hmit J0-/V0 ^ 1- In this case long before any su- 
percondcuting correlation are established the single elec- 
tron state bound to the impurity site will form. This is 
the strong coupling limit of a Kondo screening problem. 
In our case bound electron will partially screen the large 
impurity spin. In a case of a spin 3=1/2 the screening 
will be complete and th e net spin of the state will be sin- 
glet ijSakai et a/.l . ll993|) . Therefore state would evolve to 
a superconducting state with one unpaired electron that 
is bound to impurity. Ground state in a strong coupling 
limit is a non-BCS state with only pairs present. Ulti- 
mately this state is formed by the energy balance between 
superconducting and magnetic energies. Single electron 



'j'i =1. For future use we introduce 



1*0) = \{[Un+Vni^ii^l, 



(9.5) 



n>l 



Then = Vl^;^|*o). The state 7j|^'o) does not ap- 

pear inside the superconducting gap and hence is not rel- 
evant for this discussion. Note that I^I^o) is a true vacuum 



" Here the spin of the state n = 1 is determined by the sign of 
exchange coupUng Jq. We will assume it to be antiferromagnetic. 
So the electronic spin of the state n = — 1 in Eq. 19. 4i is opposite 
to the local spin S assumed to be up without loss of generality. 
Case of ferromagnetic coupling is similar. Indeed classical spin 
solution Eq. 19.11 is symmetric between Jo ^ — Jo as it contains 
only even powers of exchange. 
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state bound to a local spin provides a much stronger en- 
ergy gain ~ Jq compared to the gain due to pairing Ap. 

Crossing point and related quantities are shown in the 
Fig. 1171 This crossing point corresponds to a true quan- 
tum phase transition. 

Secondly, the crossing point occurs exactly at critical 
point in Eq. H9.1|) only in a non-self-consistent treatment 
where single particle levels provide the only contribu- 
tion to total energy. In practice the true phase tran- 
sition occurs slightly earlier. The gap suppression and 
quasiparticle interaction also contribute to free energy 
and in self-consistent mean-field approximation, we find 
that order-parameter relaxation shifts Jcrit downwards 
and the energy of the impurity-induced bound state does 
not reach zero when a first-order transition between the 
two ground states occurs. In practice the result are qual- 
itatively similar and are typically within 10 percent of 
the numerical results obtained in a selfconsistent treat- 
ment (Salkola et at, 1997). In contrast, a d-wave super- 
conductor has no quantum transition for any value of 
the magnetic moment when its quasiparticle spectrum in 
the normal state has particle-hole symmetry. The ab- 
sence of the transition follows from the behavior of the 
impurity-induced quasiparticle states which are pinned at 
the chemical potential for an arbitrarily large magnetic 
moment, see Sec. IVIII However, if particle-hole sym- 
metry is broken or if the pairing state acquires a small 
s-wave component, the transition is again possible for 
a large enough moment. The impurity moment induces 
two virtual-bound states which have four-fold symmetry 
and extend along the nodal directions of the energy gap. 

C. Particle and hole component of impurity bound state 

In this section we show that the excited states inside 
the gap in superconducting state appear in pairs at pos- 
itive and negative energies. This is a direct consequence 
of the fact that natural excitations are Bogoliubov exci- 
tations. Particle and hole coefficients of the excited state 
IvP-i) jn< j^^.j are given by the u and v components of 
the quasiparticle operators 7„, see Sec.m To be specific 
we confine subsequent discussion to s-wave case, however 
the results are applicable to any superconducting state. 

Consider two independent processes: a) electron at en- 
ergy r^o and spin down, n = —I and b) hole with spin 
up, n — 1 injected in superconductor with the same en- 
ergy f2o- Hole creation means electron with spin up is 
extracted from superconductor. This can be achieved 
by reversing the bias of the STM tip for example and 
it would correspond to the negative energy axis. Varia- 
tional wave functions that would describe these processes 
are 

l/'lll*o)jo<Jc«t = 

V'll^'o) =fi|$-i> (9.9) 

Here, to be specific we consider the case of Jq ^ Jcrit- 
This illustrates the point that in BCS like ground state 



the particle excitation with energy Hq and hole excitation 
with negative energy — Jloj aside from irrelevant normal- 
ization factors, is the same excited state, namely l^*-!). 
Therefore the poles in the density of states are always 
coming in pairs at positive and negative energies. The 
true quasiparticles in superconducting state are Bogoli- 
ubov excitations 7„ that have a finite component of par- 
ticle and hole with amplitudes w„ and u„. The strength 
of the electron absorption and emission process is con- 
trolled by these coherence factors. This is generally true 
for BCS superconductor even without impurities. For the 
case at hand, impurity states are well distinguished from 
the continuum. The two poles at ifio are part of the 
same physical excitation. We can write impurity state 
local spectral function Ai{r,uj) = — ImG'ii(r, Li;)/7r as 

Ai{Lo) = Z+S{Lo-Qo) + Z-S{io + no) (9.10) 

and the relative strength is Z+ ~ u'^-^ and ^ v'^i, 
so that the net strength of the poles Z+ -I- Z~ >= 1 as 
it should for a physical excitati on. For more details and 
references reader is referred to l|Salkola et a/.l . ll997|l . 

Analysis for Jq > Jcrit is more involved. The ground 
state wave function is now I'l'-i). Injection of the elec- 
tron with spin opposite to the spin of the bound state 
and extraction of electron with spin down will produce 
respectively 

4\<^-i) = ijUjl,\¥o) , v^_i|a>_i) = I*;) (9.11) 

with complementary annihilated states = and 

~ 0. Although the two states written in 
Eq. (|9.11l) are different, the only difference is that the 
number of Cooper pairs in these two states differ by 
one. For a macroscopically large system with number of 
Cooper pairs » 1 this produces negligible difference 
in the energies and matrix elements. Therefore, again, 
the injection of electron with spin up (in our convention) 
and extraction of electron of spin down physically will 
produce the same state. This state will have a parti- 
cle and hole projection just as we discussed in case of 

Jq ^ Jcrit- 

Similar quantum phase transition occurs in a d-wave 
superconductor even for a nonmagnetic impurity. In the 
case of particle-hole symmetric band the unitary scatter- 
ing produces a zero energy state, see Sec. IVIII Eq. ()7.1|l . 
However for the particle hole asymmetric band impurity 
state will reach zero energy and eventually will change the 
sign as a function of impurity strength. This transition 
occurs at Uq > Ucrit ~ Ml M being the chemical potential 
that leads to a particle-hole nonsymmetric band. It is 
known that single quasiparticle bound state will form at 
Uq > Ucrit and the ground state wavefunction will have 
unp aired single quasiparticle, apart from the BCS pairs, 
see l|Salkola et ad Il996l Il997j) . 
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D. Intrinsic tt phase shift for Jo > Jcrit coupling 

Here we would like to point a little noticed by impor- 
tant fact about the impurity induced modifications of the 
order parameter at the impurity site, namely a tt phase 
shift of the order parameter near impurity. 

^From the above Fig. El it follows that at Jq > Jcrit 
the self-consistent solution indicates that the phase of the 
order parameter on the impurity site is shifted by tt with 
respect to the phase in the bulk, see Fig. 1181 

Spatial extent of the tt shifted region in numerical cal- 
culation was found to be few atomic sites. Such a sharp 
change in the phase of the order parameter would cost 
a superconducting condensate energy and would not be 
preferential under normal circumstances. In the case at 
hand however, condensate energy is secondary to the 
magnetic exchange at strong coupling limit near impu- 
rity site. Physics is driven by magnetic interactions in a 
strong coupling limit. Even thought eh phase difference 
is TT the phase shift does not lead to any time reversal vi- 
olating observable effects. One can see that for the phase 
difference tt there are no superconducting currents near 
the impurity: / ~ IcSm4> — 0- These results were ob- 
tained in the self-consistent treatment within a negative 
U mo del that allows for the on-site pairing! Salkola et al., 
[l£i3)- 

We are not aware of a simple explanation of this effect. 
It appears to be general and not related to a particu- 
lar model. It is connected to the tt junctions discussed 
in the context of tunneling barrier. The tt phase shift 
is preferential in the junction in the presence of mag- 
netic impurity or ferromagnetic layer. This subject is 
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X. KONDO IMPURITY 

In all of previous discussions, we have concentrated on 
static impurities. The next two Sections are devoted to 
the examples when impurity atoms have their own in- 
ternal degrees of freedom, the dynamics of which is cou- 
pled to the scattering of conduction electrons. This dy- 
namical behavior often leads to qualitatively new results. 
Quantum dynamics is particularly important for study of 
magnetic impurities as the earliest and the best known 
example of non-trivial consequences of this dynamics is 
the Kondo effect (.Kondo. 1964) . 

Dilute magnetic impurities doped into an otherwise 
nonmagnetic metallic host have dramatic effects on the 
low temperature resistivity, susceptibility, and specific 
heat. These anomalies are associated with screening of 
the impurity spin by conduction electrons. For S = ^ 
a global singlet is formed by coupling an electron state 
to the impurity site; dynamics of spin flips is crucial 
for the formation of the singlet. For a single magnetic 



impurity in the metallic host, this is manifested as a 
crossover from the Curie susceptibility at high temper- 
atures, X = C/T with C = A^i%S{S + l)/3fcB, where 
S the magnitude of the spin and /i the Bohr magne- 
ton, to singlet behavior below a characteristic Kondo 
temperature Tk — W^exp(— l/2JiVo), where W is the 
electron half band width and J is the exchange con- 
stant. Two important notes are that a) Kondo screen- 
ing occurs only for antiferromagnetic exchange constant 
J > 0; and b) the process is non-perturbative, as is 
clear from the non-analytic dependence of the Kondo 
temperature on the exchange constant. Understanding 
of the single impurity Kondo proble in a met al required 
concerted use of the renorm ahzation group l|Andersonl 
1197(1 lAnderson et olV 1197(11). numerical renormalizat on 
group (NRG) I Wilsor . I1975 K exact solu ti ons v ia the 
Bethe ansatz ( Andrei |_|198( , Wiegmaim , |198C ) , and 



large- expansions ('Colemaii Il984 Il985l iReadi Il985t 
iRea d and Ncwns. 1983a b). M any important results 
are summarized in r ecent reviews (jCox and Zawadowskil 
ll998HHewscHlll993D . 

Kondo effect depends on the existence of the host elec- 
tronic excitations at the Fermi energy. In metals, the 
density of states is constant, which simplifies the analysis. 
If the density of states varies in the immediate vicinity 
of the Fermi surface, the effect is realized rather differ- 
ently (or not at all). Kondo effect in an insulator, with 
a har d band gap, was investigated in iQgura and Sasd 
I1993I In superconductors, however, the gap in the single 
particle spectrum only arises below the transition tem- 
perature for the Cooper instability driven by the finite 
DOS. Consequently, the two effects compete. 

Studies of magnetic impurities in superconductors be- 
gan, and largely continued, with the investigations of the 
properties of classical spin, 5* ^ 00, for which no re- 
duct ion in ina.Knitude due to Kondo screen ing is possi- 
ble l)R,usinovl 119681: IShibaL Il968t F^. |l963). The ques- 
tion of which conclusion of their analysis are robust for 
small spin values and in gapless superconductors remains 
of intense current interest. 



A. Kondo effect in fully gapped superconductors 

In Shiba-Rusinov analysis the sign of the exchange 
interaction between the conduction electrons and local- 
ized impurity spins is irrelevant. As discussed above in 
real metals antiferromagnetic exchange leads to a com- 
plete screening of the impurity below Kondo tempera- 
ture Tk , while ferromagnetic exchange does not produce 
resonance states. Consequently, treatment of quantum 
impurity spins has to bring out the differences between 
two signs of J . 

For J > opening of the superconducting gap com- 
petes with Kondo screening as both instabilities are 
driven by a finite DOS at zero energy. Qualitatively, 
it is clear that if Tk ^ T^, the impurity is completely 
screened by the time of the onset of superconductivity. 
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In contrast, for Tr- Kondo screening is suppressed 

by a decreased density of states upon opening of the su- 
perconducting gap. 

Kondo screening can be viewed as growth and diver- 
gence of the effective exchange couphng as we look at the 
properties of the system at lower and lower energies. This 
is the essence of the the NRG method, and the under- 
lying principle of successful analytical work. Therefore 
the effective exchange constant, and, with it, the phase 
shift of scattering on the impurity depends on the en- 
ergy of the incoming electron. Consequently, the effect 
of scattering varies with temperature. 

We immediately conjecture that, while a localized 
Shiba-like state exists in the gap of a superconductor, 
its energy in the Kondo limit changes upon lowering the 
temperature. Therefore the "final" , T = 0, position of 
the impurity resonance is a complex function of Tc and 
Tk- Moreover, as results of the previous section imply, 
the nature of the ground state, i.e. its spin and degen- 
eracy, depend on the relation between Tc and Tk- As 
a result, there has been increased interest in determin- 
ing the properties of the ground state, and the localized 
excited (Shiba-Rusinov) state of a BCS superconductor 
with a quantum impurity spin. 



1. Ferromagnetic exchange 

Early analytical attempt s were car- 
ried out^ jMlillcr - HartmannL Il973t 

IZittartz and Miiller-Hartmannt Il97(l in the frame- 
work of Nagaoka de coupling scheme l|HamaniJ . Il967t 
lNagaokallT96lll967|l . For J < the bound state splits 
off the band edge and was found to move towards an 
asymptotic value 



A 



1+5VS'(5+1^ 



-1/2 



(10.1) 



where g = AA^o, and A is the superconducting coupling in 
the BCS weak coupling theory. Since g <C 1 the bound 
state remains close to the gap edge for all values of J < 0. 
This qualitat ive result was late r confirmed by N RG cal- 
culations (Sa kai et"all 1993; S atori ei"all Il992l) . which 
showed that the binding energy is well approximated by 
e ~ 1 — TT^ where the renormalized exchange con- 

stant 



Jeff - 



2\J\/W 



1 + {2\J\W) ln(T4^/A)' 



(10.2) 



Therefore the ferromagnetic case corresponds to weak 
coupling and small phase shift of scattering at low tem- 
peratures. 

The ground state of th i s sys t em was arg u ed to 
be a doublet llSakai et all Il993t ISatori et all Il992t 
ISoda et all Il967^ . since the ferromagnetic interaction 
renormalizes to weak coupling and the impurity spin re- 
mains essentially free. Recently it was suggested that, 



since superconducting interaction is relevant in this sys- 
tem, and therefore above a critical coupling, Jq that 
depends on A (Jc is larger for smaller A), the ground 
state of the coupled superconductor-impurity system is a 
triplet (rn^ = 0,±1) (,Yoshioka and Ohashi .1998,) . This 
suggestion is worth exploring further. 



2. Antiferromagnetic coupling 

If the coupling is antiferromagnetic, in a normal metal 
Kondo screening corresponds to J — > oo and to phase 
shift of scatt ering, 5 7r/2. The Hartree-Fock analysis 
l|ShibaL 1197.^) is insufficient to fully describe this effect. 

The limit Tj^ <C A was considered following the 
work of Shiba (^ iillcr-Hartmann, 1973: Soda et al, 1963 
IZittartz and Miiller-HartmannL |l970), and the position 
of the localized excited state was found with various de- 
grees of accuracy. Notice that in this regime the localized 
state lies close to the gap edge as it does for ferromagnetic 
coupling. In the opposite limit, 3> A approximate 
solution for the position and the residue of the bound 
state was o btained in Refs. (Miillcr-Ha rtmannl Il973t 
IZittartz and Miiller-Hartmann. ,1970), however, the re- 
sults were inexact due to the nature of their approxima- 
tion. Later, within the local Fermi liquid appraoch, the 
energy of the bo und state in this limit was found to be 
ijMatsuural Il977|) 



where 



a : 



ttA , AeTx 

In — . 

ttA 



(10.3) 



(10.4) 



This result clearly shows that the phase shift of scattering 
depends on the ratio T^/Tk- 

The properties of the bound state, including its posi- 
tion and spectral weight, for arbitrar y values of Ty/Tp . 
were obtaine d with the help of NRG ijSakai et al\. Il993t 
ISatori et g^J . 11992) . They found level crossing similar to 
the quantum phase transition of the previous section at 
T/f/A - 0.3. For Tj^/A > 0.3 the impurity moment 
is largely quenched by the time the depletion of states 
caused by superconductivity affects screening. In that 
case the ground state is a Kondo-screened singlet, while 
the excited intra-gap state is a doublet with the spectral 
weight a « 2 for Tr: A ^ 1, corresponding to a single- 
particle state. Here a is defined from 



--lmG{uj + i6)/T: = — 

TT 2 



d{LU-Ea)+S{Lu + Eo) 



(10.5) 



On the other hand, for T^A < 0.3 the Kondo effect is 
suppressed by by the opening of the superconducting gap, 
the ground state is a doublet corresponding to a free spin 
state, while the bound excited state is a Kondo singlet. 
The spectral weight, a w 0.5 for Tr- <C A, and changes 
discontinuously at the phase transition point. 



28 



Level crossing means that the bound state is at zero 
energy for T^/ ^ « 0.3, while it is close to the gap edge 
for both Tk ^ A and Tk <C A. Numerical results show 
that the energy of the bound state is not symmetric with 
respect to the crossing point: Eq/A < 0.5 for 0.03 < 



respect to tne crossmg pomt: 
Tk/A < 1 ^Satori etallhm^. 



3. Anisotropic exchange and orbital effects 

Several more complicated aspects of Kondo screening 
in superconductors attracted attention in recent years, 
and we review them briefly, referring the reader to the 
original papers for further information. Anisotropic ex- 
change interaction, Jz J±, allows the investigation of 
the crossover between the Ising regime, J± = 0, when 
the spin-flip is disallowed and there is no Kondo screen- 
ing, and the isotropic exchange considered so far. The 
main features of the phas e diagram are discussed by 
Ijoshioka and OhashilllQQSl and new phases occur on the 
ferromagnetic side. In particular, these authors find an 
extended regime of Ising-dominated ground state even 
for J± ^ 0. In addition, they find small regions of the 
phase diagram around isotropic ferromagnetic and Ising 
antiferromagnetic lines, where there exist two localized 
intra-gap states. They also obtain a perturbative ana- 
lytic expression for the shift of the bound state energy 
due to anisotropy of the interaction. 

Using the numerical RG approach to analyse An- 
derson's model allows to interpolate between asym- 
metric magnetic scattering, Kondo problem, and 
non-magne tic scattering, includ i ng re sonance U = 
hmit l|Yoshioka and Ohashl l200(t . In par- 
ticular, the crossover from magnetically induced 
bound state to the resonance non -magnetic scattering 
regime ()Machida and ShibataL fl972|) was studied. 

Finally, so far we only discussed purely s-wave super- 
conductors. Fully gapped systems include also materials 
with a complex order parameter combining two (or more) 
out of phase unconventional gaps, such as dx^-y^ +idxy, 
or px + ipy In both of these cases Cooper pairs have or- 
bital degrees of freedom that also couple to the impurity 
spins, leading to multichannel Kondo effect. In addition, 
for p-wave pairing, the total spin of Cooper pairs is s = 1 , 
so that non-trivial changes in screening occur depending 
on whether the impurity spin 5 = 1/2 or = 1. The 
NRG analysis of the Kon do problem in this system was 
carried out very recen t lv l|Koga and Matsumotol l2002at 
iMatsumoto and Kogai |200^. They found that the two 
order parameters are indistinguishable when only I — 
impurity scattering partial wave is taken into account, 
i.e. only the depletion of the density of states due to 
the gap, rather the spin structure of the Cooper pair dic- 
tated the Kondo screening. In that case the moment of 
the ground state is determined by the orbital structure of 
the Cooper pair. However, inclusion of higher harmonics 
with I 7^ for scattering (extended impurity potential), 
leads to some novel dependencies of the screening and 



ground states on the exchange couplings. 



B. Kondo effect in gapless superconductors 

Gapless systems like d-wave superconductors with 
the low-energy quasiparticle density of states follow- 
ing the power law, N{E) cx with the exponent 
r > 0, constitute a marginal situation. The Kondo 
effect in these systems where the host single particle 
density of sta tes follows a power law, h as been studied 
intensively llBorkqwski and HirschfeldL Il992l Il994|: 
Bulla et ali. l2000l Il997l: ICassanello and Fradkin , 



19961 Il997t IChen and Javaprakashl Il99£ 



Gonzalcz-Buxton and Ingcrsent, 1998': Han et 



2002. 2004: Ingcrsent. 1996, Ingersent and SI 199S ; 
hLll993HLogan and GIossodI |200f]HPolkovnikovl 12002 : 



ito: 



Polkovnikov et al\ '2001'; 'Voita', '2001': 'Voita and Bulla, 
2001 ; Withoff and Fradkin 1990; 7M&m- et oL 
feoOll: IZhu and Tind. l2001albl). ^ Frad kin and 
co-workers ("Cassanello a nd Fradkinl Il996l Il997t 
IWithoff and Fradkin, 1990) first employed a combi- 
nation of the poor man's scaling argument and the 
large- TV approach to the case of spin-i (impurity degen- 
eracy = 2) for < r < 1, an showed that a Kondo 
effect takes place only when the electron-impurity 
exchange J exceeds a critical value. Otherwise the 
impurity decouples from the band. However, the study 

based on a nonpcrturbativc rcnormalization group 

appro ach l|Ch en and Jayaprakash, 1995; Ingerscna, 
to a spin-i identified particle-asymmetry as a key 
factor in determining the low-temperature physics. At 
small asymmetry, the critical coupling Jc above which 
the impurity moment is screened becomes so large for 
all r > i that the Kondo effect is suppressed. Away 
from the particle-hole symmetry, any quenching of the 
impurity moment is accompanied by a low-temperature 
decrease in the impurity resistivity, rather than the 
increase found in metals. The discrepancy between the 
two categories of work may stem from the mean-field 
nature of the large-A'^ method, or from the symmetry 
breaking that is implicit, for all A^ > 2, in the restriction 
that the impurity level be singly occupied. In real 
systems, the power-law variation of A'^(e) is restricted 
to an energy range |e| < Aq, with N{() ~ A^(A) for 
Aq < |e| < W . The NRG approach gave results entirely 
consistent with those known for gapped systems (The 
gap 2Ao is for r = oo limit). At the particle- hole 



^ There have also been discus sion about local moment formation in 
high-temperature cuprates (Khaluillin et al., 1997: Kilia n et oil 
J999; Nagaosa and Lee, 1997; Simo n and Var ma. 1999J. Dis- 
cussion on Kondo problem in s-wave superconductors, in 
unconventional superconductors with time-reversal-symmetry- 
broken p airing state and in insulators can be found in the lit- 
eratures_jKog^^n^^^tsumotol, l2002bL|Mateiumot^^n^j<oga , 



200i iQgura and Sasot 119 93 ; Saka i et all Il99a ISatorT et al. 



ll992HYoshioka and Ohashail99ai200i 
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symmetric case, an impurity in an insulator retains 
its moment, no matter how large J is; away from 
this symmetry, the spi n is screened provided that 
J> Jc~ 2W/\n{W/An) ^ Takegahara et oiJ.ll992D . 

Notice that considering the Kondo effect in a system 
with the power law dependence of the DOS is not the 
same as analysing the competition between supercon- 
ducting and Kondo correlations. Within the Kondo ex- 
change model, the Hamiltonian of the single magnetic 
impurity in a medium other than a superconductor, can 
be written as: 

/oo 1 J 

deN{e)eclc,„ + — ^[(L/q + ^H^c^'^ 
k,k' 

k,k' 

(10.6) 

where N{e) is the electron density of states, Nl is the 
lattice size, and all operators are for electrons. In con- 
trast, for a superconductor, we need to rewrite the in- 
teraction via the Bogoliubov quasiparticles, and enforce 
the self-consistency condition on the gap. The result- 
ing hamiltonian is rather lengthy, and follows straightfor- 
wardly from symmetric Anderson or Kondo hamiltonian, 
so that we do not give it here. We note that the formation 
and screening of the local moments in d-wave supercon- 
ductors has been i nvestigated using the vari ational wave 
function approach ijSimon and Varmalll999l) . 

The interest in the Kondo impurities in the high- 
temperature cuprates is motivated by the recent STM 
and NMR experiments around single impurities. Un- 
like most dopants, Zn, Ni are believed to substitute 
for Cu on the cooper-oxide plane and causes effective 
changes to the local electronic structure without much 
change of hole concentration. Simple valence count- 
ing suggests that if the Zn and Ni impurities main- 
tain a nominal Cu^+ charge, the Zn^+ would have a 
(Sd)^*^, S — configuration and acts as a nonmagnetic 
impurity while the Ni^+ would have a (3c?)^, S = 1 
configuration and acts a magnetic impurity. Although 
it is natural that the spin-1 impurities car ry an on- 
site magnetic moment l|Mendels et all Il999(l expected 
to give rise to the Kondo physics, the behavior associ- 
ated with the nonmagnetic impurities is completely un- 
expected. Nuclear magnetic resonance (NMR) experi- 
ments performed with nonmagnetic spin-0 (Zn,Li , Al) in 
doped cuprates l|Alloul et a^J. Il99ll: llshida et all Il993i 
119961: iMahaian et a,ll 1X994 120001: iMendels et all I1999D 
showed clearly that each impurity, itself carrying no mag- 
netic moment, induces a local -S* = ^ moment sitting on 
the nearest-neighbor Cu orbitals. It was also demon- 
strated that the magnetic properties associated with the 
substitution of these impurities strongly depends on the 
hole doping: In the underdoped regime, the moments 
retain their Curie law below the superconducting transi- 
tion temperature while near the optimal doping, the 
Kondo screening might persist even to T — > though 



strongly reduced. ^From the NMR, it is known that the 
induced moment is distributed around the impurity. We 
would like to emphasize that this moment is merely a 
particular bound state of conduction electrons near the 
impurity and the precise form of the interaction between 
the induced moment with other conduction electrons is 
a priori unknown. At this stage, we are unable to make 
a definitive conclusion about its importance for our un- 
derstanding of high-temperature cuprates. However, in a 
broader sense, the magnetic impurity embedded in super- 
conductors is a very-well defined theoretical issue. Gen- 
erally, the system Hamiltonian with a magnetic impu- 
rity consists of a d-wave BCS state Hbcs, a potential 
scattering term Hpot, and a magnetic term Timag- The 
magnetic term can described by the Anderson impurity 
model or Kondo spin exchnange model, and the impurity 
spin can be either coupled to a single site or be spatially 
distributed. For the Anderson model with the single-site 
coupling, the magnetic term is given by: 

Hmag = ^ [Vkdcl^da + H.C.] + Ed'^d^da + UdUd^ndl ■ 
ka 

(10.7) 

In the strong Ud limit, the Anderson model can be 
mapped onto a Kondo s-d exchange model through the 
Schrieffer- Wolff transformation: 

Urr.ag = ' S , (10.8) 

where Sq = \ c^^cTccnCticr' is the spin operator for 
the conduction electron at the impurity site. The cor- 
responding models for the multi-site coupling case are 
given by: 

Uniag = ^[V/dclcrO'cr + H.c] + e<i ^ d\d„ + UdTld^ridl , 
la a 

(10.9) 

and 

/ 

respectively. We point out that in the large Ud limit, the 
Schrieffer- Wolff transformation will map the model as de- 
scribed by Eq. (|10.9(l onto to a Hamiltonian nonequiva- 
lent to the Kondo model described by Eq. H10.10|l . There- 
fore, Eq. (|10.1U|I conies from different origin. 

The Anderson impurity model for a single-site cou- 
pling in d-wave supercon ductors. Eg. II10.7[I. was studied 
by Zhang, Hu, and Yu l|Zhang et all l200ll) . A sharp 
localized resonance above the Fermi energy, showing a 
marginal Fermi liquid behavior, was predicted for the 
impurity states. The same logarithmic dependence of 
self-energy and a linear frequency dependence of the re- 
laxation rate were also obtained, indicating a new uni- 
versality class for the strong coupling fixed point. Al- 
most at the same period of time, the multi-site cou- 
pling Anderson impurity model, Eq. (110.91) was consid- 
ered by Zhu and Ting (,Zhu and Tina . .2001bia) while 
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the multi-site coupling Kondo impu rity was stud i ed by 
Polkovnikov, Sachdev, and Vojta ijPolkovnikovl l2002t 
iPolkovnikov et all . l200lt) . All these works show the ex- 
istence of Kondo resonance. However, the low energy 
structure of spectral weight of the conduction electrons 
depends delicately on the local environment surround- 
ing the dynamic impurity. The on-site potential scatter- 
ing was ta ken to be either ze r o (Z han g et a i. 2001) or 
very weak (jPolkovnikovl 120021: IPolkovnikov et all l200l|) 
so that the resonance peak is lo cated very clos e to the 
Fermi energy. Zhu and Ting l|Zhu and Tind. l2001hl) 
took into account the quasiparticle scattering from a ge- 
ometrical hole, where electrons are allowed to hop onto 
the four neighbors of the impurity site, and obtained a 
double- peak structure around the Fermi energy. They 
further ijZhu and Tinel |2001(^ considered the potential 
scattering term to be in the unitary limit {U — > oo), 
and found that the Kondo resonance effect is weaved 
into that from the strong potential scattering to de- 
termine the low energy quasiparticle states. The deli- 
cate influence of the potential scattering on the Kondo 
physics as well as the local electronic structure in d-wave 
superc onductors has been re -emphasized by Vojta and 
Bulla l|Voita and BullaLl200l|) . 

To be concrete, we present a discussion based on the 
multi-site coupling Kondo impurity model, as given by 
Eq. HlO.lOfl . As demonstrated in previous sections, the 
problem of a single-site potential scattering can be ex- 
actly solved. In the Nambu space, the Green's function 
is given by 

(10.11) 

where the T-matrix due to the potential scatterer is 



T-i(zc^„) = T3/C/-g°(0,0;iw„) 



(10.12) 



and is the Green's function for the system in the ab- 
sence of impurities and has been given in early discus- 
sions. In the presence of both a potential scattering and 
a Kondo impurity, the system Green's function is found 
to be: 

g{i,j\iujn) = g{i,j;iujn) + ^fifi'g{iJ;ii^n)TK{iu!n) 



l.l' 



xg{l',j;iu;n) 



(10.13) 



Here I and I' label the sites neighboring to the impurity 
site at (0,0), and Tk is the T-matrix for the Kondo impu- 
rity. The variables (pi have different meaning depending 
on the approach to the Tk- In the large- A" approxima- 
tion (equivalent to the slave-boson mean-field approxi- 
mation), where 

^ = iujn - - ^ fifi'T3g{l, I'; iuJn)T3 , (10.14) 

Ll' 

where ipi are the complex Hubbard-Straonovich fields, 
and are determined, together with the Lagrange mul- 
tiplier A, by the saddle point solution. However, in 



the numerical renormalization group technique, when the 
strongest d-wave-like channel is considered, the variables 
are taken to be (pi = +(— )1 depending on the bond ori- 
entation. Note that this d-wave-like pattern is merely a 
band structure effect and has nothing to do with the d- 
wave symmetry of the superconducting order parameter 
of the host. The LDOS in the presence of both types of 
impurities is obtained as: 



1 



Im{Tr 



gii,i;Lu + iO+) 



I + T3 



}. (10.15) 



Figures ^1 through [2] show the LDOS for a four- 
site Kondo impurity different various strength of the 
potent ial scattering, calculat ed using the NRG tech- 
nique (|Voita and BullaL l200lD . which removes some ar- 
tifacts of the slave-boson method. It is shown clearly 
that the Kondo effect is very sensitive to the strength of 
the potential scattering. In the absence of the potential 
scattering, sharp resonance peak show up directly on the 
impurity site, and on its next-nearest neighbors with re- 
duced intens ity, whic h is consistent with the experimental 
observation ijPan et al. , 2000^) . For a moderate value of 
the potential scattering, as shown in Fig. 1201 the global 
particle-hole asymmetry changes its sign and the Kondo 
peak appears at the opposite side of the Fermi level com- 
pared to Fig. El For a strong potential scattering, the 
resonance peak directly from the impurity scattering be- 
comes dominant, and the Kondo effect is weaved into 
the overall structure of the LDOS. In this case, the in- 
tensity of the peak is strongly suppressed by the on-site 
potential scattering and a double-peak structure with en- 
hanced intensity is seen in the LDOS at the sites nearest- 
neighboring to the impur ity site. This result wa s also ob- 
tained by Zhu and Ting l|Zhu and Tind . l2001bj) based on 
the Anderson impurity model. It is expected that, in the 
unitary limit of the potential scattering, the LDOS has 
a zero intensity at the impurity site and a sharp single 
peak at its nearest neighbors. Consequently, the spatial 
shape of the resulting pattern is m ore compatible with 
the experiment af ter the filter effect (^Mar tin et alV l2002t 
IZhu et adboOOd) . as seen in Sec. IXIIII is taken into ac- 
count. It is also shown in this simple model that the large 
LDOS from the resonance state induced by the strong po- 
tential scatterer reduces dramatically the critical Kondo 
coupling, indicating the determination of the Kondo ef- 
fect by the local rather than global environment in which 
the magnetic impurity is embedded. 



XI. DYNAMICAL IMPURITIES 

A. Inelastic scattering from a single spin in li-wave 
superconductors 

We will address the inelastic tunneling features due to 
the scattering off a local spin impurity. Assume that we 
have localized magnetic atom with spin 5 on a surface of 
a d-wav e superconductor. We will treat the prob l em fol - 
lowing l)Balatskv et all 120031: iMorr and Nvberd. l2003(l . 
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Electrons in a superconductor interact with the locaHzed 
spin via point-hke exchange interaction at one site JS cr: 

H = ^e(k)cLck. + ^[A(k)4^clkx + h.c] 

k k 

Ckcr'''<T<T'Ck'o-' + ffAifiS • B , (11.1) 

k,k',CT,fj' 

where Cko- is annihilation operator for the conduction 
electron of spin a, ^(k) is the energy of the electrons, 
A(k) = Y (cos — cos ky) is the d-wave superconducting 
gap of magnitude A ~ SOmeV in typical high-Tc mate- 
rials. The local spin S is a \S\ = 1/2. We focus here 
on the effect of the Zeeman splitting of the otherwise de- 
generate local spin state in the external magnetic field B 
with splitting energy ujo = ul = g^sB. Below we use 
a mean field description of superconducting state at low 
temperatures T <^ Tc- Assuming field B <C Hc2 we will 
ignore the orbital and Zeeman effect of the field on the 
conduction electrons 

We are interested in a local effect of inelastic scattering 
of electrons. Thus only local properties will determine 
the conduction electron self-energy. Results we obtain 
will also hold for a normal state with linearly vanishing 
DOS, such as a pseudogap state of high-Tc superconduc- 
tors. In the case of a normal state one would model nor- 
mal pseudogap state with a single particle Hamiltoninan 
■^0 = Sk^WcLckfT with N{uj) uj. 

Because of the vanishing DOS in a d-wave 
superconducting state Kondo singlet formation 
occurs only for a coupl ing constant exceed- 
ing some critical value llCassanello and Fradkinl 

1996t iGonzalez-Buxton and Ineersend . Il998t 

Withoff and Fradkinl 'l990). Fbr a particle- hole 
symmetric spectrum Kondo singlet is not formed for 
arbitrarily large values of J. Another situation where 
Kondo effect is irrelevant is the case of ferromagnetic 
coupling J. This allows us, quite generally, to consider 
a single spin in a d-wave superconductor that is not 
screened and we ignore the Kondo effect. 

In the presence of magnetic field B||z spin degener- 
acy is lifted and components of the spin S||z and S _L B 
will have different propagators. It is obvious that only 
transverse components of the spin will contain informa- 
tion about level splitting at ojq = ojl- We have therefore 
focused on , S~ components only. The propagator in 
imaginary time r is x{t) — {TtS^ {t)S^ {0)) with Fourier 
transform and continuing to real frequency Xo{^) — 
For free spin we have {Sz) = tanh(a;o/2T)/2. 



For more general case of magnetic anisotropy this does 
not have to be the case. To be general we will keep {Sz)- 



To minimize the orbital effect of magnetic field one can apply it 
parallel to the surface of superconductor. The magnetic field is 
penetrating the surface sheath on the scale of penetration depth 
so that its effect on the conduction electrons for d-wave SC is 
small. 



We begin with evaluation of the DOS correction due 
to coupling to localized spin. Self-energy correction is: 



l](a;0 = J'T G(k,c^, ~ aOx+"(aO , (11-2) 

k,0„ 

where GO(k, uji) = [lui - mW^^i? " C'(k) - A^k)]-^ 
is the particle Green's function in d-wave superconduc- 
tor, G-i = G(0)~i - E, F°ik,LUi) = [A(k)] [(ic^,)' " 
£_^{k) - A2(k)]-i; rii = 2'kIT is the bosonic Matsub- 
ara frequency and uji — (21 + 1)ttT;1 = 0,1,2... is the 
fermionic frequency. Using spectral representation and 
analytical continuation onto real axis iujn —^L0 + i6 we 
find for imaginary part of self energy S(ijj) : 

Iml](ijj) = —J^ {Sz)lniG{Lj—tJo)[nF{tJ—Lj())—nB{i^o)^^] , 

(11.3) 

where n_F(w) — l/[l + exp{(3uj)],nB{ui) = l/[exp(/3w) — 1] 
are Fermi and Bose distribution functions. This local 
self-energy leads to the modifications of the DOS. In 
this solution we treat the self-energy effects in G to 
all orders, i.e. G in Eq. H11.3|l is full Green's function 
— Gq^ — S(a;) and solution for E is found self- 
consistently for a local vibrational mode. The mod- 
ifications of the superconducting order parameter and 
bosonic propagator were ignored in this calculation. Re- 
sults are presented in Fig. [221 To proceed with analytic 
treatment, unless stated otherwise, we limit ourselves be- 
low to second order scattering in S. Difference between 
self-consistent solution and second order calculation are 
only quantitative and small for small coupling. Correc- 
tions to the Green's function G{r,r',uj) = G'^(r, r',w) -|- 
G"(r, 0, w)S(tj)G°(0, r', lu) + F"(r, 0, uj)J:{uj)F*°{0, r, lu). 
For simplicity we define K{T,uj,ujo) = — [ri^? (w — ) ~ 
nsi^o) ~ 1 — Q{lo — loq)] which becomes a step function 
at low T <^ ojQ, the limit we will focus on hereafter. Cor- 
rection to the local density of states as a function of posi- 
tion comes from the correction to the bare Green's func- 
tion G° ; SNir,uj) = 1 / Trlm[G° (r, 0, lu)^ (uj) G° (0, r, lu) ± 
F°{r,0,Lo)Y,{uj)F*'^{0,r,iu)], where keeping it general, 
the plus sign corresponds to the coupling to the local 
vibrational mode and minus - to the spin scattering re- 
spectively. The strongest effect will be at the impurity 
site. For on-site density of states we have: 



dN{r^O,uj) TT 



2'-^ ~ 1-^0 



No 
'2a;, /A 



,iJSNo)^ A 



K{T, LO, LOo) 

(11.4) 



■ = 27r (JbNo) K [1 , UJ , ujo) m 



No 



X In 



4A 



4A 



+ (luo -too), - |A| , (11.5) 



\uj + ujo — A\ 

where we used for on-site Green's function G°(0, 0, w) 
■^ln(^)-f zttH; 



7Vo[¥ln( — ) + ittH], for w < a and we retained 
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only dominant real part of G". In opposite limit 
~ A we retained only dominant imaginary part of 

GO(0,0,w) = iTTN{uj) = -2i7Voln(J^^). At r = we 
have F°(0,0,w) = 0. Complete DOS N{uj) and deriva- 



tive 



dui 



are shown on Fig. 1221 For arbitrary position 
A^(r, uj) we would have to add a Friedel oscillation factor 



that 



A(r) = [|G°(r,^)P ± |F«(r,.)P] ^ ,,^^^^y.^^^^/,y. 
describes the real space dependence of the Green's func- 
tion on distance for small <C A. Here r_L||ki?j_ is the 
component of r = (r^,r||) that is along the Fermi sur- 
face near the nodal point of the gap and r||||kj^[| is the 
component perpendicular to the Fermi surface at nodal 
point. Existence of the nodes in d-wave case results in 
the power law decay of A(r) in all directions and it has a 
four fold modu lation due to gap anisotropy (See detailed 
discussions in (jSalkola et ai, 1997)), also Sec.lWl). The 
final result is our Eqs. H11.4I11.5|) . 
It follows immediately that 



f dl I dl 
" dvl dV 



5N{r 



(JS'iVo)2^^e(F 



"dV^ 



0,V)/No ^ 



(11 



Here we have used the fact that the derivative of {uj — 
Ci;o)6('^ — wq) with respect to lu yields 0(w — wq)- Thus in 
a d-wave superconductor and in a metal with vanishing 
DOS N{ll!) = Nq^ one should expect a step discontinu- 
ity in d?I/dV^ at the energy of a local mode with the 
strength J^Nq (see Fig. I22|l. This result is qualitatively 
different from the case of conventional metal. For metal 
with energy independent DOS we have from Eq. (|11.4|) 
for T < Wo 

^^SN{r^O,V)^J^N^e{V-LOo), (11.7) 

and the second derivative will reveal a delta function 
d'^I/dV'^ ~ J'^N^5{ijj - Wo) The effect in d-wave super- 
conductor is clearly smaller than correction to DOS in a 
normal metal with the same coupling strength. 

For completeness we also give the result for inelastic 
scattering in a metal with the more general DOS N{lo) = 
l/7rImG°(0,0,w) = {u;/A)'^No with power 7 > that is 
determined by the microscopic properties of the material. 
Then, from Eqs. (|1 1.311 1.4|l we have for w < A: 

^w/iv - ^^(^ = 0' ^)/^o ~ - ^oVeiv - wo) , 

(11.8) 



Taking imaginary part of G(r, w) we arrive for T ^ wo 
at: 

Iml](w) = ^(S,)[\u^{r ^ 0)\HiLj - cvo - E^) 
2wo 

-H|i;„(r = 0)|2j(w- Wo + £;„)] ,w > C(.11.10) 

At negative w < one has to replace wo — > — wo in 
Ea. (|ll.lU|l . For example, consider a magnetic impu- 
rity resonance in d-wave supercond uctor at energy Wim,n, 
such as a Ni indu ced resonance ijHudson et all l200lt 
ISalkola 6^^119971) . Then only the term with resonance 
level Ea = Eimp will dominate the sum over eigenstates 
a in the vicinity of impurity site. Inelastic scattering 
off this impurity induced resonance will produce addi- 
tional satellite split away from the impurity level by wo, 
see Fig. Sharp coherence peaks will also produce 
split satellites. Again, for a local phonon mode one gets 
a similar splitting of impurity level with wo now being 
the phonon energy. 
, These results suggest the possibility of single spin de- 
pJ,ection as one monitors the feature in d^I/dV^ as a 
function of position and external magnetic field. If we 
take experimentally seen DOS Ng 1/e V with JNq ~ 
0.14, A = SOmeV l)Hudson et all . l200l|) and assuming 
the field of ~ lOT we have wo = ImeV (corresponding 
to the Zeeman splitting of ^ ImeV in a magnetic field, 
we have from Eqs. H11.4I11.6|I 



6N{r = O,w)/A^o - 10" 



- Wo 



e(w-wo). (11.11) 



Result is expressed in terms of the relative change of 
DOS of a metal A'o. For observation of this effect one 
would have to sample DOS in the vicinity of eV = wq oc 



B. Assuming w 

^4t7/4^ 10~^. Expressed as a relative change of DOS 



dV I dV 

of a superconductor N{uj 



Wo = Wo we have from Eq. 
as a relative chai 
Nquj/ A effect is: 



5—1— 
dV' dV 



6N{r = 0,w)/7V(wo) 



io-2i^i^e(w-wo) 

Wo 



"dV^ 



(i/-wo)T-ie(F-wo) 



Depending on the value, we get divergent singularity at 
Wo for 7 < 1, or a power law rise for 7 > 1. In case of 
7 = 1 we recover the result for d-wave superconductor 
and for a pseudogap normal state. 

Quite generally one can express the results in terms of 
the spectrum of superconductor. We can write Iml](w) 
using spectral representation for G(r,w). In supercon- 
ducting case, using Bogoliubov Ua{r), Va (r) for eigenstate 
a, we have 



(11.12) 

It is of the same order of magnitude as the observed vibra- 
tional modes of localized molecules in inelas tic electron 
tunneling spectro scopy STM, lETS-STM (,Hahn and Hoi 
I2OOII: IStipe et all [1998). The satellites at A wo pro- 
duce the effect on the scale of unity and clearly seen even 
for small coupling. The important difference is that for 
localized spin the kink in DOS is tunable with magnetic 
field and this should make its detection easier. 



B. Localized vibrational modes in d-wave superconductors 

When a localized vibrational mode is coupled to elec- 
trons in a superconductor, the Hamiltonian is given by: 

H = 5l6cLcfe. + ^[Afc4^cl,^+/i.c.] 



G(r,w) 



E 



- Ec+i6 



ka 



[11.9) 



-gJ2{b^ + b)ci 



"Ocr'^Oo- , 



(11.13) 
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Up to a second order in coupling constant g calcula- 
tions are very similar if not identical to the ones in 
the previous section. F or more detailed analysis see 
l|Morr and NvbereL I200I . 

The proposed extension of the inelastic tunneling spec- 
troscopy on the strongly correlated electrons states, such 
as a d-wave superconductor and pseudogap normal state 
would open up possibilities to study the dynamics of local 
spin and vibrational excitations. The DOS in these sys- 
tems has a nontrivial energy dependence of general form 
N{uj) ~ U!^,'J > 0. This technique could allow for a Zee- 
man level spectroscopy of a single magnetic center, thus, 
in principle, allowing a single spin detection. There is a 
feature in dl/dV ~ {uj—loo)'^^^Q{uj—u!o) near the thresh- 
old energy loq that indicates the inelastic scattering. One 
also finds strong satellite features near the gap edge due 
to coherence peak for a superconducting case. The singu- 
larity is a power law and qualitatively different from the 
results for a simp le metallic DOS (Hah n and Ho, 200,1,; 
IStipe et a^J . ll998(l . For the relevant values of parameters 
for high-Tc the feature is on the order of several percents 
and makes the feature observable in these materials. Sim- 
ilar predictions are also applicable to the local vibrational 
modes, where ujo becomes a vibrational mode frequency. 



XII. INTERPLAY BETWEEN COLLECTIVE MODES AND 
IMPURITIES IN d-WAVE SUPERCONDUCTORS 

Conventional wisdom dictates that local probes have 
poor momentum resolution, and therefore cannot identify 
collective modes that exist at specific wave vectors. How- 
ever, recent p rogress in relat ing the STM observations 
with ARPES l|McElrov et all [2003) by using a Fourier 
transform of the image taken over a large area opens 
up the possibility to study connections of STM measure- 
ments with other spectroscopies. Since the ability of the 
STS to identify the Fermi wave vector is most naturally 
explained by the sensitivity of the technique to Friedel 
oscillations in the electron density due to impurities, we 
now look into the DOS features arising from the interplay 
between dynamic scattering of of the collective mode and 
static disorder. 

Neutron scattering is one of the spectroscopic measure- 
ments which revealed a resonance excitation, the so called 
42 meV peak that is commonly present in a number of 
materials. It has been pr oposed tha t STM may be used 
to detect neutron mode l|Zhu et all l2004a(l . The main 
tool for such a measurement is the Inelastic Electronic 
Tunneling Spectroscopy (lETS), see Sec. lXIl Specifically 
the proposal is to use Fourier transform (FT) tunneling 
maps and search for features that represent an lETS sig- 
nature. 

We limit our consideration to the example of a 
well defined mode at wave vector (tt, tt) with energy 
ujQ. The ultimate goal is to detect the bosonic spec- 
tral function, be it magnetic spin mode or some lat- 
tice modes. Recent efforts indicate possible connec- 



tion between the kink in ARPES data on quasipar- 
ticle dispersion and phonon mod e s and , possibly, su- 
perconductivity llDamascelli et all l2003t iGweon et all 
\2004 iLanzara et all |2004|) . It was suggested that 
the full Eliashberg function in frequency and mo- 
mentum space may be ex tracted from ARPES data 
l|Vekhter and Varmal 1200^ . and the challenge is to de- 
sign a similar procedure for its determination from the 
FT lETS STM. Efi'orts to relate the data from ARPES, 
STM, and transpor t m easure ments in cu p rates have 
recent ly intensified (|Scalapino e t all 120041: IZhu et all . 
I2004cl) . 

We need the impurity scattering to produce features 
that can be analyzed using Fourier transform. We con- 
sider weak Born scattering from distributed scalar poten- 
tial U (vi) = Ui, i is the lattice site index. We find indeed 
that disorder potential can strongly modify the features 
as seen in local DOS. One of the interesting findings is 
that characteristic wave vectors of the impurity potential 
Uq — 2 • Ui exp(iq • r^) do play a crucial role in defining 
characteristic wave vectors of the DOS modulation. 

To explain the method we point out that FT STM 
data do contain momentum information. However all 
the analysis up to date on FT STM data was done 
within the framework of elastic scattering that considers 
the natural quasipartic le excitations at the Fermi surface 
llHoffman eit qiJ . l2002H) . No inelastic scattering processes 
off the collective mode were included in the analysis. To 
consider the scattering of quasiarticles off the collective 
mode one has to explicitly keep track of the self-energy 
effects that result from scattering. In this case simple 
noninteracting quasiparticle picture is not adequate. In- 
elastic scattering of quasiparticles requires us to consider 
excitations off shell, for example to consider excitations 
at energies that are typically A -|- fio TOmeV. At these 
energies the combination of the Fermi surface effects, 
typical wavevectors of the collective mode and typical 
wavevectors of the random potential all enter in deter- 
mining the momentum of the inelastic tunneling features 
as seen in FT STM. 

We limit ourselves to the second order scattering be- 
tween carriers and bosonic excitations and at this level 
there is no conceptual difference in the method as applied 
to spin or phonon bosonic mode. 

We start with a model Hamiltonian describing two- 
dimensional electrons coupled to a collective spin mode 
and in the presence of inhomogeneity: 

Ti- = Hbcs +'Hsp +'Himp ■ (12.1) 
Here the BCS-type Hamiltonian is given by Tiscs — 

Ek,^(£k - M)cLcka + Ek(^k4^cLk| + AJC-kiCkt), 

where cj^^ (cko-) creates (annihilates) a conduction elec- 
tron of spin a and wavevector k, ek is the normal state 
energy dispersion for the conduction electrons, /x the 
chemical potential, and Ak = -^(cosfc^; — cosfcj,) the 
d-wave superconducting energy gap. The coupling be- 
tween the electrons and the resonance mode is modeled 
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by an interaction term Hsp = gj^i^i ' where the 
quantities s^, and Si are the couphng strength, the 
electron spin operator at site i, and the operator for the 
collective spin degrees of freedom, respectively. The dy- 
namics of the collective mode will be specified below. The 
quasiparticle scattering from impurities in the Hamilto- 
nian is given by, iJ^^p = J^icr Uicj^Cia ,where is the 
strength of the impurity scattering potential. For sim- 
plicity, only the case of nonmagnetic scattering is con- 
sidered here. By introducing a two-component Nambu 
spinor operator, '^i = (ci|, c|j^)'^, one can define the ma- 
trix Green's function for the full Hamiltonian system, 
G{i,j;T,T') = -{Tr[^i{T) (g) ^'](t')]). An algebra leads 
to the full electron Green's function with impurity scat- 
tering: 

j' 

-g[l\i,f;ito,,)g2iij\j;iLOn)] ■ (12.2) 

Here Q'^^^ is the dressed Green's function, with its Fourier 
component given by: 

[^(")]-i(k;zc.„)= ' 



-Ak - 1121 iu)n + ?k - S22 

(12.3) 

where Ck = £k ^ = (2n + 1)'!tT is the fermionic 

Matsubara frequency for fermions. When the inelastic 
scattering of quasiparticles from the collective mode oc- 
curs, the self-energy is obtained to the second order in 
the coupling constant as: 



ig'T 



q 



(12.4) 

where x(q;*^i) is the dynamical spin susceptibility 
Xij{T) = {Tr{Sf{T)SJ{Oj)) and ili = 21ttT the bosonic 
Matsubara frequency, and Gq is the bare superconduct- 
ing Green's function. We have also assumed that the 
d-wave pair potential is real. For a single-site impurity, 
the equation of motion for the full Green's function can 
be exactly solved. For the case of multiple impurities and 
especially the inhomogeneous situation, some approxima- 
tion needs to be taken. We consider here the case of the 
weak impurity scattering limit so that the Born approx- 
imation can be used. With this ansatz, we arrive at: 

C/ii(«,j;ia;„) ^ g^f^\i,j;iu!n) + 6gii{i,j;iu!n) (12.5) 
with 

5gii{i,j;iuJn) = Uj'[g["j^ {ij';iujn)gn W J;iuJn) 

(12.6) 



^l2^(i,/;*^^)^2i^(j',j;«'^«)] 



The LDOS at the i-th site, summed over two spin com- 
ponents, is 



p{Yi,E) = Im5ii(i,i;£'-|-i7) 

TT 



where 7 = 0+. We are especially interested in the cor- 
rection to the LDOS from the impurity scattering: 

5p{v,,E) ^ --liiY5gii{i,i-E + i-i) . (12.8) 

TT 

Its Fourier transform is: 
<5p(q,i?) = ^Jp(^,i?)e-^'^•■■• 

i 

u, 



Nni 



^[^^';Hk + q;£; + ^7)^^?^(k;^; + ^7) 



-g['i^*{k~<l;E + ^J)g[Tik■,E + ^^) 

-g[l\k + q;E + zj)gi°\k;E + ij) 

+g[T{k - q; + ^l)GiTi^■, E + *7)] (12.9) 

where the Fourier transform of the scattering potential 
is 



f^q = E ' 



(12.10) 



The corresponding Fourier (wave- vector) spectrum is de- 
fined as 



P(q,i?) = |<5p(q,i?)| 



(12.11) 



Up to now discussion and formulation are quite gen- 
eral and can be used to study the effects of any dynamic 
mode once the susceptibility x is known. Consider now 
specific case of magnetic mode with susceptibility taking 
a phenomenological form (based o n the inelastic neutron 
scatt ering observations), see also (jEschrig and NormanL 

l2onn|) : 



x(q; i^i) 



1 



+ r^o 



(12.12) 



where we denote the wavevector Q = (vr, tt) and 
the mode energy by fio- This form captures the es- 
sential feature of resonant peak observed by neutron 
scattering experiments in t he superconducting state of 
cuprates ()Zhu et alV l2004a|) . For the normal-state en- 
ergy dispersion, we use £k = —2t{coskx + cosfcj,) — 
At' cos cos ky , where t and t' are the nearest and next- 
nearest neighbor hopping integral. Unless specified ex- 
plicitly, the energy is measured in units of t. We choose 
t' — —0.2 to model the band structure of the hole-doped 
cuprates. Since the maximum energy gap for most of the 
cuprates at the optimal doping is about 30 meV while 
the resonance mode energy is in the range between 35 
and 47 meV, we take Aq = 0.1 and ilo = 0.15 (i.e., 
I.5A0). The chemical potential (/x « —1.15) is tuned 
to give an optimal doping value 0.16. To mimic the in- 
trinsic life time broadening, in our numerical calculation 
we take 7 of Eq. H12.7|l to be O.OSAq. A system size of 
N = N^, X A'^,, — 256 X 256 is taken in the numerical 



(12.7) calculation. 
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Figure plots the density of states in a clean system 
for various values of coupling constant g. Without the 
electron-mode coupling {g — 0), the density of states is 
peaked only at the maximum gap edges ±Ao. When 
there exists the electron-mode coupling, for example 
g — 0.2 as shown in the figure, several new effects emerge. 
As a result of the additional anomalous self-energy intro- 
duced through the coupling, the maximum gap edge is 
renormalized to Arem which is larger than Aq. More im- 
portantly, the singularity in the quasiparticle self-energy 
causes additional poles in the Green's function, and new 
peaks show up at the energy zLEr = ±(Aj. -I- fio). A 
strong implication of this result is that regardless of 
the renormalization of the energy gap, the position of 
the new peaks relative to the superconducting coherent 
peak is shifted by fio- In addition, with the appearance 
of the peaks away from the gap edge, the intensity of 
the superconducting coherent peaks is reduced such that 
the sum rule is obeyed. The intensity of the peaks at 
negative energies is stronger than that at positive ener- 
gies since the van Hove singularity is below the Fermi 
energy. These results, for the clean case, are consisten t 
with ear lier studies of th e ARPE S llAbanov et al. 
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2002 ; ICampuzano et'oH. Il999l: iDessau et 
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Norman and Ding, "1998; "Shen and Schrieffeil Il997|) and 
bos (Abanov and Chubukov. 2000,) . The shift of states 
due to inelastic sc attering is also expect ed for scattering 
off of local mode l|Balatskv et a/.l.l2003(l . 

We now turn to the Fourier spectrum in the presence 
of disorder. An accurate description of this problem re- 
quires an extremely high energy and spatial resolution. 
Therefore, a very large system size should be considered. 
For the quasiparticle scattering off a single impurity, a s 
considered in the work by Zhu et al. l|Zhu et a^J . l2004a(l . 
one can first calculate the LDOS within a small window 
aroimd the impurity site in a very big system size (e.g., 
N = 1000 X 1000), then perform a Fourier transform over 
the window size and even with masking of sites. This pro- 
cedure and flexibility does not exist in the case of disor- 
der and inhomogeneity with multiple scattering centers. 
Summation over the wave-vector in the Brillouin zone 
constrains us to consider a moderate system size. 

In Fig. 1251 we show the Fourier spectrum at the energy 
—Er for the coupling constant g — 0.2 with a structurless 
scattering potential Uq — Uq. This structureless Uq cor- 
responds to a single-site impurity in real space. In the ab- 
sence of the electron-mode coupling [g = 0), the Fourier 
spectrum has a strongest intensity at q = (0, 0) and its 
equivalent points, and a moderately strong weight along 
the edges of the square around q ~ (tt, tt). When there 
exists the electron-mode coupling, as shown in Fig. 1251 
the spectrum has the strongest intensity point at four 
corners of the square, and a moderately strong intensity 
at the four ridges of the square. This implies that as 
one takes line cuts along the diagonal the first feature at 
±Er always will be at the wavevector {tt — 6,tt — 6) be- 
fore the (tt, tt) point. Regardless of whether the electrons 
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are coupled to the collective mode, the spectrum has a 
minimum in intensity at q = (tt, tt), which is different 
from the results when the pre-dominant Friedel oscilla- 
tion is filtered out of the Fourier transform. When the 
scattering potential has a clear stucture, which might be 
relevant to the inhomogeneity in high- Tc cuprates, the 
pattern of the Fourier spectrum changes dramatically. 
As an ansatz, we propose the following stucture for the 
scattering potential: 



ql + 4[cos2(9^/2) + cos2(9j,/2)] 



(12.13) 



where the parameter go describes the extent of the peak 
at q = (tt, tt). This structure factor has a highly non- 
trivial consequence on the Fourier spectrum of the lo- 
cal density of states. This is because overall modula- 
tion (5/9(q, uo) (X J7q and for C/q peaked at (tt, tt) FT DOS 
5p{i\) will also be peaked at this wavevector. As shown in 
Fig. 1261 the Fourier spectrum now has a strongest inten- 
sity at the four ridges, which are located along the diag- 
onals of the first Brillouin zone, and a secondly strongest 
intensity at q = (tt, tt). Except in the four small lobes 
around the corners, the spectrum has moderate intensity 
inside the square. 

Here we have considered the situation where all impu- 
rities have identical potential scattering strength. If one 
dopes a strong scattering impurity, for example Zn sub- 
stituted for Cu, into a high-Tc cuprate which has inhomo- 
geneity coming from weak potential scatterers, the pat- 
tern of the Fourier spectrum of the local density of states 
is mostly determined by the scattering off the strong im- 
purity. 

To summarize, the main results of this section are : 
1. the energy of the inelastic feature is at Er — Aq -I- 
r^o '^-^ lOmeV for optimal doping. Given that the gap 
is position dependent in observed spectra, this energy 
will be position dependent. 2. the typical wave vectors 
along diagonals where the inelastic features are seen are 
determined by a number of factors: the momentum of 
the disorder potential U{q), the doping and positions of 
the "diamonds" seen in Figs. IJSl and 3. the first 
feature that is seen in FT lETS STM in our calculation 
is occurring at wavevectors that are inside (tt, tt) vectors. 



XIII. SCANNING TUNNELING MICROSCOPY RESULTS 
A. STM results around a single impurity 

The STM has established itself as a remarkably power- 
ful and versatile tool for studying the electronic proper- 
ties of solids. Its remarkable energy and spatial resolution 
makes it particularly well suited for materials character- 
ized by small energy and short length scales. It measures 
the tunneling current varying with the voltage bias and 
the tip positions. In the tunneling Hamiltonian formal- 
ism, the differential tunneling conductance — the deriva- 
tive of the current with respect to the voltage bias, is 
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given by 

^cx-l dc.^|TkpA.(k,c.)/;^(c.-ey), (13.1) 

where fpo is the Fermi distribution function, and 
Acr(k,uj) is the electron spectral function of the 
sample. The tunneling matrix clement, |Tkp — 
X]q l-^^kqP^tip(q, w), where Afkq is the matrix element 
representing the overlap of the electronic states on the tip 
and sample. Using a tip with featureless DOS around 
the Fermi energy, we can assume |Tkp is energy inde- 
pendent. If we further assume a k-independent tun- 
neling matrix element, one can find that the tunnel- 
ing conductance is proportional to the local density of 
state at the tip position, which we have chosen to be 
the origin, p{eV) = - / dujJ2k^^MK^)f'i^^ ~ eV). At 
zero temperature, it is simply given by the imaginary 
part of the electronic Green's function we have used 
heavily for the discussion throughout the work, that is 
p{v^eV) = — i ImGo-(r, r; = eV). Here we have 
labeled the tip position by r. 

The experimental attempts to detect and accurately 
resolve the sub-gap features in the density of states in su- 
perconductors with impurities have a long history. This 
feature was observed in the planar junctions doped with 
magnetic impurities in earlier years. However, a direct 
observation of the sub-gap states induced by a magnetic 
impurity did no t occur until late 1990 's. In 1997, Yazdani 
and co-workers l) Yazdani et aZ.Lll997|) deposited adatoms, 
Mn, Gd and Ag, on the (llO)-oriented surface of a super- 
conducting Nb sample, and examined the local electronic 
structure around them. Figure f?7l shows their STM tun- 
neling spectrum measurement. The main fundings are: 
(1) The local density of states has no much difference 
for the tunneling through Ag impurity atoms and far 
away from them, since Ag impurity atoms are believed 
to be non-magnetic in nature; (2) The LDOS at energies 
less than the Nb's gap is enhanced when the tunneling 
is through Mn and Gd magnetic atoms. The enhance- 
ment happens at the length scale of lOA, indicating the 
bound nature of the impurity states; (3) The LDOS spec- 
tra are asymmetric about the Fermi energy. Within the 
framework of the BdG theory, the authors used a two- 
parameter magnetic impurity model, where the electrons 
are coupled with the impurity through an magnetic ex- 
change interaction J and experience a nonmagnetic po- 
tential scattering U . The obtained results, consistent 
with the Yu-Shiba-Rusinov prediction and more recent 
theoretical works, fit the experimental data. However, 
the model calculation required the value of J of the or- 
der of 4 eV, in the strong coupling limit, and failed to 
capture the detailed spatial dependence of the spectra 
around Gd site. 

The pioneering STM experimental research on the lo- 
cal electronic structure around single defects and im- 
purities in high- Tc cuprates was carried out by two 
grou ps led by Eigl er at I BM Almaden Research Cen- 
ter ijYazdani et alV Il999(l and Davis at UC Berke- 



ley ijHudson et alV Il999j) . Since the high- Tc cuprates 
have a d-wave pairing symmetry, even nonmagnetic im- 
purity scattering would affe ct the superconduct ivity. By- 
ers, Flatte, and Scalapino (i Bvers et al\ . Il993l) were the 
first to suggest the use of STM to study the spatial vari- 
ations of the tunneling conductance nea r impurities. In 
particular, it was theoretically predicted ijBalatskv et alV . 
1995; Salkola et ai, 1996) that quasiparticle resonance 
states are induced around a nonmagnetic impurity in 
a c?-wave superconductor, in striking contrast to s- 
wave systems. The sample used by Yazdani and co- 
workers is overdoped Bi2Sr2CaCu208 with a supercon- 
ducting transition temperature of 74K and a transition 
width of 3K. The sample used by the Berkeley group 
is Bi2Sr2CaCu208+5 with a transition temperature of 
87K and a transition width of 5K. The STM experi- 
ments were operated at 5K and 4.2K, respectively. The 
STM spectroscopy on these samples, which were nom- 
inally undoped with known impurities, shows clearly 
the enhancement of the local density of states near the 
zero voltage bias in regions where the chemically in- 
duced defects in the sample are located. The experi- 
ments provided a strong evidence for the existence of 
low-energy quasiparticle resonance states around single 
nonmagnetic impurities, as predicted theoretically. The 
asymmetric or splitting of the measured resonance near 
the zero bias may come from the fact that the particle- 
hole symmetry may be broken by impurities and de- 
fects locally or the underlying realistic band structure 
of BSCCO (F latte and B vcrs, 1998; Zhu et al. ^OOOa). 
However, in these two experiments, the location in the 
crystal and the identity of these scattering centers are 
unknown. Moreover, since the enhancement of the lo- 
cal density of states at these scattering centers is not 
dramatically large and the coherence of high-Tc super- 
conductors is so short, it is very difficult to investigate 
in detail the local electronic structure around them at an 
atomic scale. 

New STM study on the impuri ty effects in BSCC O was 
reported by the Berkeley group ijPan et al l l2nonH) . The 
samples were Bi2Sr2Ca(Cui^2;Zn2;)208+5 single crystals 
with intentionally doped with x = 0.6% Zn. The crystals 
have the transition temperature of 84K and a width of 
4K. To search for low-energy quasiparticle states associ- 
ated with the Zn atoms, the authors first mapped the 
differential tunneling conductance at zero sample bias in 
a larger window of the surface and found a number of 
randomly distributed bright sites corresponding to the 
areas of high LDOS. Then they measured the tunneling 
spectroscopy exactly at the center of a bright scattering 
site. As shown in Fig. |5S1 the spectrum showed a very 
strong DOS peak at the energy Q. — —1.5 ± 0.5 meV. 
The peak intensity can be up to six times greater than 
the normal-state conductance. At the same time, the in- 
tensity at the superconducting coherence peak is strongly 
suppressed, indicating the almost complete local destruc- 
tion of superconductivity. These phenomena is consistent 
with the theoretically predicted characteristics of quasi- 
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particle scattering off a nonmagnetic unitary impurity in 
a d-wave superconductor. The strong intensity of the 
near-zero-bias peak allows the authors to give a close 
inspection of the electronic structure around the Zn im- 
purity. As shown in Fig. 1201 the STM differential con- 
ductance imaging at Q = —l.bmeV exhibits two novel 
features: Firstly, it has strongest intensity directly at the 
impurity site and local maxima at the sites belonging 
to the sublattice containing the impurity site, while lo- 
cal minima at the sites belonging to the other sublattice. 
Secondly, the intensity decays much faster along the gap 
nodal direction than along the bond direction. These new 
features are totally unexpected. The theory based on a 
potential scattering model would predict a vanishingly 
small intensity at the impurity in the unitary limit. The 
first feature motivated theorists to study the electronic 
structure a round a Kondo impurity in a d-wave super- 
conductor fPolkovniko v et al 1 1200 H IZhang et all l200lt 
|Zhu and Ting, 2001b), and consider the importance of 
the BiO layer which i s the exposed surface llMartin et all . 
12002; .Zhu and Tin'gll2?mcl IZhu et o4l200ncD . 

The STM study on the local electronic structure 
around a magnetic Ni atom in BSCCO was also r eported 
by the same Berkeley group ijHudson et a^J . l200Hl . It was 
found that there two spin-resolved resonance states in- 
duced by the Ni atom, in contrast to the case of Zn atom 
in previous experiment where only a spin-degenerate res- 
onance state is induced. The energy of four resonance 
peaks in the tunneling spectrum are, ±f2i and ±f22, with 
= 9.2 ± 1.1 mev and ^2 = 18.6 ± 0.7 meV. The ex- 
perimental result is reasonable agreement with a theo- 
retical m odel with bo th nonmagnetic and magnetic scat- 
tering fS alkola et all . [1997: Tsuchiura et at, 2000). By 
substituting the values of f2i^2 and the maximum super- 
conducting e nergy gap An = 28 meV into the theoretical 
formula (,Salkola aij . ll997|) : 

" ~2Nf{Uo ± J) In \8Nf{U ± W)\ ^^^'^^ 

with Np the normal-state density of states at the Fermi 
energy, U and W the strength of nonmagnetic and mag- 
netic scattering, one can find NpU = —0.67 and NpW = 
0.14. This result indicates that, despite Ni atoms possess 
a magnetic moment, the scattering off them is dominated 
by potential interactions. In addition, the experiment 
also showed that the intensity at the gap edge in the 
tunneling conductance directly at the Ni impurity site 
is almost unaffected, in comparison with that far away 
from the impurity, supporting the scenario that th e high- 
Tc su perconductivity is magnetically mediated ijPined . 
Il997|) . 



B. Spatial distribution of particle and hole components 

Spatial distribution of tunneling intensity clearly ex- 
hibits alternation between positive and negative bias, see 



for example. Fig. 1291 It appears as a rotation of an im- 
purity induced cross upon changing the sign of the bias. 
Since the effect is so explicit in the images we will address 
it here. 

Apparent rotation of the impurity intensity can be un- 
derstood as a result of interplay between particle and hole 
components of the Bogoliubov quasiparticle. This effect 
is a general property of superconductivity and is seen in 
both s-wave and d-wave superconductor s llHudson et al\ . 
I200lt iPan et all l2000at lYazdani et aH Il997|) . see also 
Sec. IIXI The Bogoliubov quasiparticles, that are na- 
tive excitations in superconductor, have both particle and 
hole component. The sites where there is a large parti- 
cle components will have large intensity on positive bias 
site and hence will be show up as bright sites on positive 
bias. Sites with large hole component will be bright on 
negative bias. Fig. 

Let us define the respective amplitudes of particle and 
hole amplitudes of the Bogoliubov quasiparticle, Un{i) 
and Vn{i) for site i and for particular eigenstate n. They 
obey the normalization condition |u„(i)p-|-|v„(z)p = 
1 for any fixed site i. Consider now a site where, say, 
Un{i) is large and close to 1 for particular eigenvalue. It 
follows therefore that for the same site the w„(i) would 
have to be small, since the normalization condition is 
almost exausted by |m„(i)P term alone. Similarly, for 
the sites where has large magnitude, u„(i) would 

have to be small. Large M„(i) component would mean 
that quasiparticle has a large electron component on this 
site. Hence the electron will have large probability to 
tunnel into superconductor on this site and the tunnel- 
ing intensity for electrons positive sample bias will be 
large. Conversely, for those sites the hole amplitude is 
small |wn(«)||Mn(i)| and the hole intensity negative sam- 
ple bias will be small. Similarly, for sites with large hole 
amplitudes |wn(«)| |u„(i)| the electron amplitude will be 
suppressed and this site will be bright on the hole bias. 
Therefore if there is a particular pattern for the large 
particle amplitude (sampled on positive bias) on cer- 
tain sites i, the complimentary pattern of bright sites 
for hole tunneling (on negative bias) will develop as a 
consequence of the inherent particle-hole mixture in su- 
perconductor. This is the physics behind what appears as 
the cross rotat i on upon bias switch, see n in experiments 
ijHudson et al\ . 120011: iPan et adl2000a(l . see Fig. ^ 



C. Fourier-transformed STM Measurement 

The ingredient of the Fourier-transform STM tech- 
nique is to collect a large set of tunneling conduc- 
tance data (at a fixed voltage bias) in the real space, 
and then to perform a Fourier transform. Th is tech- 
nique was first applied bv lHoffman et a/.l . l2002al to study 
the quasiparticle states generated by a quantized mag- 
netic vortex in the mixed state of slightly overdoped 
high-Tc superconductor, Bi2Sr2CaCu208+5. A Cu-0 
bond-oriented "checkboard" pattern with 4ao period- 
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icity. The 4ao modulation periodicity is one half of 
that (8ao) of the field-induce d SDW modulati o n ob- 
served in neutron sca ttering l|Khavkovich et oil l2002t 
iLake et "0/1120011 12002') on other cuprate materials. This 
field-induced "checkerboard" pattern has been inter- 
preted as the induction of two-dimensional spin den- 
sity wave around the vortex core where the supercon- 
ductivity is suppressed (Anderse n and Hedceard. 200^ 
iTakieawa eTai] . 120031 Izhu et all l2002j) . the nucleation 
of the antiferromagnetic order bro ught about by loca l 
quantum fluctuations of a vortex I Franz et all l2002|) . 



and the frozen of d-wave hole pairs into a crystal by 
the magnetic field l)Chen et "^■ 120021) . Most of the the- 
ories rely on the proximity of the system to a quan- 
tum critical point so that it is very sensitive to external 
perturbations. The same kind of checkerboard pattern 
has also been predicted around a single strong impu- 
rity with induced local moment in the optimally doped 
cuprates (Chen and Ting. .20031. .2004 .Liang and Lea 
l2002HZhu et^l2002^ 

The challenge comes from the experimental obser- 
vation of a similar checkerboard patte rn even at zero 
field in the same dop i ng regim e ( Hoffm an et all l2002lj: 
iHowald et d\ . 2003; McEhov et at, 2003*) . In the 
exper iments by Hoffman et aL . ■2002b' .McElrov et all 
I2OO3I the Fourier analysis of the images of the energy- 
dependent modulations yields the dispersion of wavevec- 
tors. Instead, Howard and co-workers observed the exis- 
tence of static striped density of electronic states, i.e., 
the four-period peaks in the Fourier transform of the 
data are present at all energies, including very low en- 
ergies. One can understand these two effects separably. 
Those peaks showing energy dispersion co mes from the 
a quasiparticle scattering from impurities J Bvers et all 



119931: IWang and LeeL 12003: IZhang and TingLi2003ll2004^ . 

A heuristic model b ased upon the e lectro nic band struc- 
ture is as follows llMcElrov et all |200l: In BSCCO, 
four nodes exist in the superconducting gap A^. Be- 
low the gap maximum Aq, the contours in fc-space along 
which quasiparticle exist at a given energy are banana- 
shaped, as shown in Fig. The quasiparticle density 
of states at energy E = lu, p{E = uj) is proportional 
^'^ Isk^ui l^fe-^fel^^'^'*'' where the integral is performed 
over the contour Ek = lo. Each 'banana' exhibits its 
largest rate of increase with energy, |Vfci?fe|~-^, near its 
two ends. Therefore, the primary contributions to p{E) 
are from the octet of momentum-space regions centered 
around kj(i?), j — 1,2,..., 8, at the end of the 'ba- 
nana's. (Red circles in Fig. In the presence of im- 
purities, quasiparticles will be elastically scattered. A 
quasiparticle located in momentum-space near one ele- 
ment of the octet is highly likely to be scattered to the 
vicinity of another element of the octet, because of the 
large density of final states there. For each in a rep- 
resentative octet, there are seven characteristic scatter- 
ing wavevectors. This octet-model then predicts a to- 
tal of 56 scattering wavevectors. Of these, 32 constitute 
a complete set of inequivalent wavevectors and there- 



fore 16 distinct ±q pairs can be detected by Fourier- 
transformed sca nning tunneling spec troscopy. The ex- 
perimental data ijMcElrov et 0/1120031) are in good agree- 
ment with this model. From the material point of view, 
although no external impurities were introduced in a con- 
trolled manner into the sample of these experiments, the 
source of quasiparticle scattering may be closely related 
to t he experiment a lly observed nano s cale inhomogene- 
itv iHowald et all l200lt iLang et all l2002t iPan et all 

\mM- 

In cont rast, the obser v ed no n-dispersive LDOS mod- 
ulations ijHowald et all l2003|) should be interpreted 
by invoking a static (or fluctuating) charge- or spin- 
ordered state llKivelson fit oil l2003t iPodolskv et all 
I2OO3I: iPolkovnikov et all 120031) . The emergence of a 
competing ordering is due to the quantum criticality 
with or without the aid of the inhomogeneity in the 
sample. An even stronger evidence of the compet- 
ing ordering is provided by recent observation that the 
electronic states at low energies within the pseudogap 
state in Bi2Sr2CaCu208+(5 exhibit spatial modulations 
wit h an energy-inde pendent incommensurate periodic- 
ity ijVershinin et o/.l . I20b4^ . Theoretically, a complete 
microscopic model with all these elements has not yet 
been developed. 



D. Filter 

We point out that for unitary scattering impurity in 
any model it is difficult if not impossible to produce large 
intensity on the impurity site. Unitary scattering pro- 
duces a node in the wave f unction. Y et, experimentally, 
the impurity site is bright (|Pan et ai . 2000b). One ex- 
planation is that the image seen by STM is not the real 
intensity of the impurity state that is buried below the 
exposed layer in STM experiments. One needs to have a 
model on how intensity is transmitted to the top layer. 
Th e idea of a filter the n comes in naturally. Martin et 
al. ijMartin et all\206^ proposed an idea of filter that in- 
tensity as seen at the top layer by STM is a convolution 
of initial intensity due to impurity scattering and filter 
function that comes from the effective hopping matrix 
element between CuO planes, tk oc \ cosk^ — cosfc^p. 

The reasoning goes as follows. In order to tunnel be- 
tween layers it is advantageous to involve tunneling be- 
tween s-wave orbitals that extend out of the Cu-0 plane. 
These orbitals are off the chemical potential and virtual 
hopping on these orbitals would bring large energy de- 
nominators in any perturbation scheme. Still it pays to 
engage s-wave orbitals of Cu because one gains on the ex- 
ponential overlap factors between s-orbitals in adjacent 
planes. The electronic orbitals near chemical potential 
are essentially d^-2_yi orbitals of Cu (hybridized with p- 
orbitals but we ignore this hybridization). The d^2_y2 or- 
bitals on the impurity site are orthogonal to the s-orbital 
of the impurity site. The next available s orbitals are on 
the nearest Cu sites. Therefore electron hops virtually on 
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to Px or py orbitals of nearest O and then onto s-orbital. 



The amphtude for the hops Cu d,j.2-y2 —> O p^ 



Cu 



would be different for hops along horizontal and vertical 
directions, as one can verify from Fig. H34I) . 

For example the hopping to the Cu site on the right 
one would get for and amplitude Ai i^^iy) '■ 



{d.i\px){Px\3i+x) 



(-l)exp(ifca;a) 



A. 



i,i~\-y ^ 



[Ep — Ed] [Eg — Ep] [Ep — Ed] [Eg — Ep] 

{di\Py){Py\Sl + y) ( + 1 ) CXp ( ^ fcj; ) 



[Ep - Ed] [E, ~ Ep] [Ep - Ed] [Es - Ep] 



(13.3) 



in the second equation we considered a plane waves that 
describe the states without impurity scattering. One im- 
mediately can see that the signature for the horizontal 
and vertical amplitudes is opposite in sign regardless of 
the phase assignment of p-orbitals for pure case and for 
general amplitudes of the states produced by impurity 
scattering. For a quantum mechanical process to hop 
from one site to nearest neighbor s-orbitals one would 
have to add the amplitudes: 



Afot = A. 



A,, 



A,. 



cos{kxa) — cos(fcj,a) . 



A 



i,i-y 



(13.4) 



Again, second equation refers to the pure plane wave 
analysis to make contact with the bands structure calcu- 
lations for the tunneling matrix element ijAndersen et al\ . 
IT995il . Upon hopping on the s-orbitals electron hops to 
the next layer and retraces its path exactly in reversed 
sequence as described above. Therefore the amplitude for 
the hopping will be proportional to the square of the Atot- 
The net hopping matrix element has the from consistent 
with |(ir2_,,2| modulations: 



[Atot\^ cx [Ai.i+x 



A 



iS+y 



Ai^i-y I 



(13.5) 



This particular filter is directly connected to the inter- 
plane hopping matrix element obtaine d with in the band 
structure calculation llAndersen p.t all 1199,^ . However 
the idea that one has to involve the s-orbitals of the Cu- 
O plane is relevant also for an exposed Cu-0 layer as one 
would need to tunnel from the s-orb itals of the t i p ont o 
relevant s-orbitals of the Cu-0 plane l|Misra et oiJ . l200^ . 

Alternative filter due to blocking of certain hop- 
pin g matrix e l ements has been considered by Zhu et 
al. ijZhu et all . l2000c^) . For an analysis of the local ef- 
fects of impurity one need to consider a local tunneling 
matrix elements that has to connect impurity orbitals to 
s-orbitals on neighboring Cu atoms that have a great- 
est overlap between Cu-0 layers. The net effect of the 
filter is to produce large spectral intensity on an impu- 
rity site and nearest neighbor sites to be dark. More 
recently, the measured quasi-continuous data has been 
convert ed to a set of LDOS d efined on a two-dimensional 
lattice dWang and Hi A 120041) . which is suitable for a rig- 
orous comparison between the tight-binding model stud- 
ies and the STM experimental data. 



Another important observation one can make by com- 
paring STM and local N MR results available in Li doped 
YBCO superconductor ijBobroff et ad 1200 1|) . In case of 
Li impurity NMR revealed that maximum intensity in 
NMR signal comes from four nearest neighbor Cu sites 
and is quite localized near impurity. This observation 
would be consistent with the notion that strongly scat- 
tering impurity produces large density of states on near- 
est sites. The crucial difference between NMR and STM 
is that NMR observation does not require electronic tun- 
neling. Magnetic field is measured instead. Hence there 
is no filter to apply to native electronic states in Cu- 
O plane to obtain NMR real space distribution. There- 
fore, depending on the type of measurement one might 
need or need not to use the filters. Details depend on the 
nature of the measurement. 



XIV. AVERAGE DENSITY OF STATES IN 
SUPERCONDUCTORS WITH IMPURITIES 

The Green's function formalism is well suited to 
the analysis of the combined effect of many uncorre- 
lated impurities in the bulk of a superconductor. The 
first treatment of the superconducting properties us- 
ing this technique was given by Abrikosov and Gor'kov 
(^brikosov and Gorkov, 1963) ^ pioneering paper. 
The basic assumptions underlying such calculations were 
given in Sec. 1111. lj| After averaging over different impu- 
rity distributions following Eq. (|3.14|) . the translational 
symmetry in the system is restored, and therefore the 
Green's function takes the general form 

G"i(k,w) = iw„ - e(k)r3 - Aoa2T2 - E (14.1) 
= iuj - e(k)T3 - Aa2T2. (14.2) 

Here the second line explicitly takes into account the ma- 
trix structure of the self-energy, E. The superconducting 
gap in the presence of impurities is determined by the 
self-consistency condition, Eq. H2.21|l . which reads here 



cut, 



(14-3) 

The transition temperature is the temperature at which 
a non-trivial solution of the self-consistency equation first 
appears. Together with the recipe for computing the 
self-energy these equations form a general basis for treat- 
ing ensembles of impurities in superconductors. We note 
here that we always ignore the contribution of E3, which 
is equivalent to the renormalization of the chemical po- 
tential. This is always allowed in computing the den- 
sity of states, although the corrections may need to be 
taken into account in e valuating the response functions 
l)Hirschfeld et a?J . ll988|) . The basic assumption for com- 
puting the self-energy is that, in addition to neglect- 
ing the interaction between spins on different impurity 
sites, see Sec. IIII.BI we can neglecting the interference 
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effects of scattering on different impurities (which the 
order (p_fO~^; where I is the mean free path). 



where Aqo is the gap in the pure material at T = 0. 
The transition temperature is determined fr'om 



A. s-wave 

1. Born approximation and tiie AG Theory 

In a seminal paper Abrikosov and Gorkov analysed the 
effect of the impurity scattering on superconductivity in 
the Born approximation. We briefly review this analysis 
to compare its outcome with the results of theories going 
beyond Bor n approximation. We follow the treatment of 
iMakiL 119691 

Consider a general impurity potential combining the 
potential and the magnetic scattering, 

t/,™p(k - k') = Upotik - k')T3 + J(k - k')^ • a, (14.4) 

where a is defined in Eq. (|3.5(l . AG considered the self 
energy in the Born approximation, 



S(cj,k) = riirnp 



(27r)3 

Integrating over k' we find 



C/„„p(k-k')G(k',c^)[/,„,p(k'-k). 

(14.5) 
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(14.7) 



where the potential (rp) and spin-flip (ts) scattering 
times are given by 

- = n,rapNo [ dh\Upot{k~k')\'', (14.8) 

Tp J 



- = n™piVo5(5 + l) / dn\J{k-k')\\ (14.9) 

Ts J 

Here we averaged over all possible directions of the im- 
purity spin. 

In the absence of spin-flip scattering both A and uj are 
renormalized identically, and it follows from Eq. (|14.3() 
that the gap remains unchanged compared to the pure 
case. This is in accordance with Anderson's theorem. 
The spin flip scattering time (which violates the time- 
reversal symmetry) enters the equations for uj and A with 
the opposite sign. Therefore, introducing u = w/A, we 
find 



(14.10) 



It follows that the gap in the single particle spectrum is 
Egap = A(l - (Ars)~2/3)3/2 f-Qj. ^ 2, and vanishes 
for Ats < 1. This gapless region starts at the value of 
pairbreaking parameter a 



T, 



In^, (14.12) 



where ipix) is the digamma function and Tco is the tran- 
sition temperature of the pure material. Consequently, 
superconductivity is destroyed {Tc — 0) when 



= t;^ = ^Tco/27 = Aoo/2 > a', 



(14.13) 



where 7 k, 1.78. As a' « 0.912q;c AG predicted that a 
regime of gapless su perconductivity exists for a ra nge of 
impurity scattering ijAbrikosov and Gorkovl IT960|) . This 
was first confirmed in experiments by Woolf and Reif 
iWooff and ReiA 1196,^1) . 

The evolution of the density of states with 
increasing disorder w a s in v estigated in deta. il 
l)Ambegaokar a nd Griffi^ Il965t ISkalski et all. Il964() . 
and is shown in Fig. 1351 For a < a' a. hard gap in 
the single particle spectrum persists up to the critical 
impurity concentration, as shown in Fig. 1361 This result 
is clearly at odds with our discussion in Sec. IVII which 
shows that even a single magnetic impurity creates a 
localized state in the superconducting gap. 



2. Shiba impurity bands 

In the AG theory the impurity concentration and the 
strength of the exchange coupling contribute to the sup- 
pression of superconductivity as a single pairbreaking pa- 
rameter, a = = (2nimp/7rA^o) sin^ (Sq cx riimpJ^ S{S + 
1) for isotropic exchange, see Eq. (|14.9() . This is a result 
of the Born approximation; in general, the phase shift 
So and the concentration of impurities n-i^p are separate 
variables that control different aspects of impurity scat- 
tering. For example, in the limit of dilute concentration 
of strong magnetic impurities, the AG approach yields a 
small scattering rate, and a single-particle spectral gap 
virtually identical to that in a pure limit. On the other 
hand, we have learned that in this regime each impurity 
is accompanied by a bound state with the energy below 
the gap, and therefore we expect a finite number of these 
sub-states to exist in a superconductor. This section ad- 
dresses this dichotomy. 

Analysis of the strong scattering regime requires go- 
ing beyond the Born approximation, and here we use 
the self-consistent T-matr i x app roach ijHirschfeld et all 
Il986t ISchmitt-Rink et all Il98q) . where the self-energy 
S(p,a;) = n, 



Hmp'E'p^P'f and 



fp,p' = Up,p' + J dpiC/p,piG(pi,cj)fp,,p- 



a' = ^ = Aoo exp(-7r/4). 



(14.11) 



(14.14) 



Following the treatment described in Sec. IVII we anal- 
yse the pairbreaking in different angular momentum 
channels. The effective pairbreaking parameter in the 
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Z-th channel is ai 



5(1 — ef)/(27rA^o), where e; is the 



position of the corresponding bound state, see Eq. H6.1U|I . 
In analogy with the AG treatment, we find that the ratio 
Un = pn/A(u;„) sat isfies the equation ijChaba and Nasi 
ll972HRusinOTl[T96ai 



/=0 ' 



K 

7/2 



(14.15) 



where the gap is determined self-consistently from 
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The width of the gapless regime now also depends on 
the details of scattering. For / — channel only the 
:ap vanishes \yhen the pair breaking exceeds the value 
Rusinovi Il969t IShibaL ITqbI 



= 2e2 exp[-7re2/2(l 

Q^cr 



eo)]. 



(14.20) 



This equation should be contrasted with Eq. p4.10|l . The 
pairbreaking parameter, ai now depends separately on 
the position of the single-impurity resonance state, e/ and 
the impurity concentration, in contrast to the AG theory. 

The growth of the impurity band has been investigated 
for the sp herically symmetric case of purely magnetic 
scattering (IChaba and NagiL Il972t Ir usinovl Il969t IShibaL 
[r96&) . The critical concentration of impurities at which 
the transition temperature vanishes is obtained by set- 
ting Tc = in the gap equation, 



ln^=^(l/2 + a/2^r,)-V'(l/2), 

c 

where now llGinzberd . lT97l 

a ^ ^{21 + l)ai. 



(14.17) 



(14.18) 



Since the gap equation is identical to that considered by 
AG, the critical pairbreaking, a^r = Ao/2. However, now 
the critical concentration of impurities depends on the 
phase shift of scattering by individual impurities, and on 
the position of the single impurity resonance, see Fig.lSTI 



In the Born approximation (weak scattering) the bound 
state moves to the gap edge, — 1, and we regain the 
result of Abrikosov and Gorkov. For stronger scatter- 
ing, eo < Ij the realm of gapless superconductivity is 
enhanced compared to the AG theory. As higher order 
harmonics are included, the threshold at which the den- 
sity of state s at the Fermi en ergy becomes non-zero shifts 
even lower (iGinzbersi lTpTfll). 

For Z = in the limit ao ^ A the width of the impurity 
band around Eq is estimated to he W = (8aoA)^/^(l — 
eo)^^''', and therefore varies as n 

™p isElHilllHS)- There- 
fore if the resonance state at Eq is sufficiently close to the 
gap edge, the concentration, cq, at which the top of the 
impurity band merges with the continuum above A is 
smaller than the critical concentration, c', at which the 
bottom of the impurity band reaches the Fermi surface 
and the superconductor becomes gapless (Shiba, 1968), 
see Fig.|2Hl In fact, the AG result is an extreme example 
of this behavior when the states due to individual impu- 
rities are infinitely close to the gap edge, and therefore 
the effect of increasing impurity concentration is an ap- 
parent decrease of the gap until the onset of the gapless 
behavior. 



J 



3. Quantum spins and density of states 

In the fully quantum treatment of the impurity spin. 
Sec. 13 we discussed the competition between gapping 
the density of states due to superconductivity, and the 
onset of the Kondo screening of the impurity moment. 
The main conclusion was that, in contrast to classical 
impurity spin, the position of the bound state is not sim- 
ply given by the value of the bare exchange coupling but 
depends sensitively on the ratio Tk/Tc- Once the po- 
sition of the bound state is established, in the limit of 
independent impurities one can consider the growth of 
the impurity band in analogy with the previous section. 
As discussed previously, for ferromagnetic coupling of the 
impurity to the conduction electrons, the bound state is 
always close to the gap edge, the scattering is weak, and 
we can expect that the Abrikosov-Gor'kov theory gives 



correct results. 



When Kondo screening is effective, for antiferromag- 
netic coupling, the behavior of the density of state and 
the transition temperature was studie d a series of papers 



bv Miiller-Hartmann and co- 


workers (IMiiller-Hartmann 


19731: iMiiller-Hartmann and Zittartzl Il971 


Schuh and Miiller-Hartmann 




19781 IZittartz et al 



The main new result was the prediction of the 



re-entrant behavior for small T^jTc ^ 1- In that case 
the phase shift of the scattering increases upon lowering 
temperature, but remains moderate at enabling the 
transition to the superconducting state. Upon further 
decrease in temperature, scattering becomes stronger 
and suppresses superconductivity in a range of phase di- 
agram of Fig. 1211 Finally, at lowest temperatures below 
Tr-, the system re-enters local Fermi liquid regime with 
weak scattering and superconductivity may re-appear. 
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While further work l)Jarrell Il990t iMatsuura et adll977|) 

cast doubt on the existence of the third transition, 
region of two solutions for Tc{nimp) was confirmed by 
theoretical studies. In particular, a combination of quan- 
tum Monte Carlo technique with Eliashberg equations 
gave the dependence of the re-entrance transition on 
the electron-phonon coupling constant, while accounting 
non-perturbatively for the Kondo effect ( JarreH , 1990[l . 
see Fig. I^Hl Moreover, the initial decrea se of with 
increasing impurity conce ntrat i on is fast l|,Tarrel]L Il99flt 
iMiillcr-Hartmann and Zitt art j. Il97ll). and depends on 
the coupling strangth IjJarrelll |l99(]f) . The behavior 
of the den sity of states in this l ir nit was investigated 
in d etail l|Bickers and ZwicknaglL Il987t I.Tarrell et oiJ . 

The overall shape of the transition temperature 
as a function of impurity concentration with re- entrant 
transition was observed in (LaCe)Al2 alloy series (jMaolel 
Il97,'^ . 



B. d-wave 

For completeness we briefly consider the growth of 
the impurity band with finite concentration of impuri- 
ties. As was mentioned above, scalar (non-magnetic) 
impurities are pair-breakers for any nonconventional su- 
perconductor, and substantially change the low-energy 
spectrum of superconducting quasiparticles. This prob- 
lem has been addressed in great detail in the framework 
of the self- consistent T-matrix approximation (for ex- 
ample, see llBalatskv e t oL. 1994; Gorkov and Kalugin, 
1985 [ iHirschfeld and aoldenfeldL .199,1- .Hirschfeld et all 



T9M n988t lT.eet 119931 iPethirk a"nd PmeA 
Schmitt-Rink et al\ . \l9H (f)) . which leads to the finite den- 



sity of states at the Fermi level. Here we briefly review 
the main steps and give results for the quasiparticle scat- 
tering rate and low-energy density of states for complete- 
ness. 

For finite impurity concentration, the self-consistent 
Green's function, averaged over impurity positions, was 
given in Eq. H3.21|l as 



G-i(k,c^) = Go-i(k,c^)-S(w). 



(14.21) 



with S(a;) 



-iT(u;). In the case of particle-hole sym- 



metry l|Hirschfeld et all Il988|) , and unconventional gap 
(defined by us as having a zero average over the Fermi 
surface, see Sec. ^ the only non-vanishing component of 
the T-matrix is proportional to tq. 



Toito) 



(14.22) 



The T-matrix has to be determined self-consistently with 
5oM) = 2%EkTrG(k,..)To. 

Solution of this equation leads to a finite density of 
states at the Fermi level. This result was first ob- 
tained for Born scat tering ijGorkov and Kaluginl 1198,'Tt 
lUd ea and Rici[l98l, leading to an exponentially small 



N{0)/No « 4T^AQexp(-2Aor), where r is the normal 
state scattering rate. The results are much more dra- 
matic for unitarity scatte r ing (c = 0) ijHirschfeld et all 
Il986t ISchmitt-Rink et all Il986|) . when straightforward 
algebra yields 



7 ~ ^71imp(Ao/7riVo), 



(14.23) 



where 7 = —lmY,{u! 0) is the scattering rate for low- 
energy quasiparticles. For w < 7, the density of states is 
determined by impurities and is finite: A^imp(0)/Afo — 
2'j / IT Aq. The characteristic width of the impurity- 
dominated region is ~ 7 oc .yriimp. 

The origin of the finite density of states is the impurity 
band, grown from the impurity-induced states (consider 
c = 0). Scaling of the impurity bandwidth 7 (x y^nhnp 
has been obtained earlier for the case o f paramagnet ic im- 
purities in an s-wave superconductor l)Shibal . [l968|) . The 
fact that 7 oc y/n\^ is obeyed in the case of a d-wave su- 
perconductor with scalar impurities as well is consistent 
with the claim that the low-energy states in a disordered 
d-wave superconductor are indeed formed from the bound 
states at finite concentration. Many questions about the 
exact nature of the interference between impurity sites in 
unconventional superconductors remain as of now unan- 
swered. Wc briefly reviewed some of the relevant work in 
the introduction, but do not discuss it in depth here. 

Notice that the results above are for isotropic impu- 
rity scattering. Anisotropic impurities may preferentially 
scatter electrons between regions with the same, or close 
values of the gap, so that the scattering is inefficient 
in suppressing Tc- For gen eral impurit y phase shifts 
this has been considered by IChoi Il999l while for the 
model wit h dominant s mall angle scattering in cuprates 
lAbra hams and Varmal i200Q„') the effect was considered 
by Kee \Ke A\2o7)W . 



XV. OPTIMAL FLUCTUATION 
A. Introduction 

So far we concentrated on discussing the effect of a 
single impurity on its immediate surrounding and on the 
combined effect of an ensemble of scattering centers on 
the spatially averaged properties of a superconductor. In 
the case of a single pairbreaking impurity the charac- 
teristic length is simply the superconducting coherence 
length, ^0- In the Abrikosov-Gorkov approach the gap 
is assumed to be uniformly suppressed. If the coher- 
ence length is short, this assumption breaks down as the 
energy cost of local suppression of the order parameter 
becomes smaller than the cost of uniform reduction of 
the gap. In this case again the length scale of this sup- 
pression is of the order of ^o- These results are obtained 
by carrying out a standard impurity averaging procedure 
at the mean field level, i.e. averaging over all the pos- 
sible configurations of impurity atoms ijAbrikosov et all 
I1963D . 
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It is clear, however, that some physics is missing in such 
an approach. Among all the realizations of the impurity 
distribution in a sample of size Lq there exist regions 
where the local impurity concentration, on some charac- 
teristic scale L Lq, differs significantly from the aver- 
age concentration, n. If the local impurity concentration 
is sufficiently high, for L > superconductivity may be 
locally destroyed of sufficiently suppressed to generate a 
bound quasiparticle state at an energy E <^ Aq. 

Of course, such regions are rare. There is a high en- 
tropy cost in creating an impurity droplet with the con- 
centration significantly different from the average and 
hence the probability of encountering these regions is 
small. However, the states localized in these droplets 
make a non-perturbative contribution to the density of 
states averaged over the entire sample, N{E), and quali- 
tatively modify its behavior compared to the mean field 
(Abrikosov-Gorkov and Shiba) treatment. Quite dramat- 
ically, they make any s-wave superconductor with a small 
concentration of magnet ic impurities (Ar^ 3> 1) gapless 
l)Balatskv and TrugmanL ITqQTi) . It is due to such a dra- 
matic modification that the interest in these "tail" states 
stretching below the mean field gap edge has peaked in 
recent years. 

The problem of tail states did not originate in the 
study of superconductivity. The contribution of regions 
of anomalous impurity concentration to the net density 
of states below the gap edge was first considered in d oped 
semiconductors by Lifshitz ijLifshitzl . IT964albl Il967j) . He 
was the first to show that such rare impurity configura- 
tions create a local profile in the Coulomb potential that 
can have bound states, and therefore gives rise to the 
non-vanishing density of states below the bottom of the 
band. Eg. Henceforth the states localized in the droplets 
of impurities have become known as 'Lifshitz tails', and 
have been extensive l y studied { Halpcrin and j^axi Il966t 
IVan MieghCT3.ll992tlZittartz and Langer..l96e 



While, in retrospect, it is natural that inhomogeneities 
lead to a low-energy tail in the density of states in su- 
perconductors in much the same way, little attention has 
been paid to this problem until the paper b y Balatsky and 
Trugman l|Balatskv and Truemanl Il997j) . Their study 
was stimulated by the experimental observations that 
the tunneling density of states in s-wave superconductors 
with magnetic impurities is far greater at low energies 



than the Abrikosov-G orkov theory sugg ests ( Bader etai\ . 
ll975tlEdelstehlll967tlWoolf and Beif 119651) . A number 



of theoretical studies of the tail states followed, and this 
topic is now a subject of active interest. 

Below we first briefly review the physical picture of the 
tail states in semiconductors, and then describe how it 
is applied to the subgap states in superconductors with 
impurities. 



B. Tail states in semiconductors and optimal fluctuation 

In a semiconductor there are two distinct situations: 
a) heavily doped, and b) lightly doped with impurity 
atoms. In the former case a localized tail state with en- 
ergy E < Eg forms in the impurity-rich region, and the 
extent of its wave function greatly exceeds the average 
distance between individual shallow sites. Therefore the 
exact impurity potential can be replaced by a smooth 
function, averaged over regions containing many impuri- 
ties. The probability of realization of the potential with 
the "right" energy of the bound state among all the pos- 
sible impurity distributions determines its contribution 
to the DOS. In the latter case the number of impurity 
sites needed to form a bound state depends on how deep 
below the band edge the energy of such a state is. For 
example, if each impurity binds an electron at energy Ei , 
while E2 is the energy of the state bound by two impu- 
rities on neighboring lattice sites, to obtain a localized 
state below Ei but above E2, one simply needs to find 
a region where the two impurities are at a particular fi- 
nite distance from each other. The probability of finding 
such an i mpurity pair determines the density of states 
llLifshitzllT964hlll967D . As we go to energies below E2 
we need to position three impurities etc. 

For energy, E, the most probable (albeit still very rare) 
configuration of impurities that creates a potential U, 
with a bound state from the solution of Schrodinger's 
equation [Hband + U]4> = ElU]^), such that £\U] = E, 
and therefore contributes the most to N{E) is called the 
optimal fluctuation. Given the probability density for 
the impurity potential, P[U], and the density of states in 
this potential. 



N{E)^ VUP[U]S{E - £[U]), 



(15.1) 



the optimal fluctuation is obtained by using the saddle 
point approximation and minimizing the resulting func- 
tional with respect to U. This approach finds the cheap- 
est (from the entropy consideration) impurity potential 
that creates a bound state at E. Therefore it optimizes 
the non-uniform impurity distribution (fiuctuation from 
the uniform average) to the given energy, hence the name 
"optimal fluctuation" . The general technical difficulty of 
minimization lies in its essential nonlincarity: the opti- 
mal potential depends on the wave function of the parti- 
cle in this potential. 

Let us consider the example of many uncorrelated shal- 
low impurity centers forming an extended potential. It 
is described by the Gaussian probability density. 



P[U] (X exp 



Saddle point approximation for Ea. ll5.1() gives 
N{E) 



In- 



^0 



-S[Uopt], 



(15.2) 



(15.3) 
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where the optimal fluctuation is obtained by minimizing 
the functional 



(15.4) 



with respect to the potential U and the Lagrange mul- 
tiplier A. At the simplest level it is sufficient to con- 
sider only the potentials where S\U] = E is the lowest 
energy state in the potential U; fluctuations where E co- 
incides with the higher eigenstates are exponentially less 
probable. In a semiconductor the kinetic energy of the 
quasiparticles is /2m* , where m* is the effective mass. 
Consequently, in a potential well of depth U (all energies 
are measured from the band edge) and size L the energy 
of the localized state is of the order of [/ + l/(mL^) = E 
{h = 1). In the optimal fluctuation E ^ U ~L~^,so that 
the action for such fluctuation is tSft/j » L'^U'^ /Un, or 
\n\N{E)INo\ ~ -\E\^-'^I'^ jUl llHalperin and LaxlllQBfit 
[Lifshitz, 1964lJ. Importantly, the characteristic size of 
the optimal fluctuation, L (x increases as the en- 

ergy approaches the band edge, while its depth, \U\ ~ \E\ 
decreases: the potential becomes more shallow and ex- 
tended. 

More formally, note that the energy of the bound state 
is the expectation value of the hamiltonian over the wave 
function of the bound state, V'(r), is equal to E, 



£[t/] = (i/) = (V^|-^+C/|7/^)=i? 

2m 



(15.5) 



Minimization of the action in Ea. H15.4|l with respect to 
IJ dictates that 

U{x) = -A[/2(^| — K-) = -AC/oV(a:), (15.6) 

while minimization with respect to A dictates that, in this 
potential, the bound state is at energy E, i.e. (setting 
m* = 1 for simplicity) 



A[/oV'(r) 



V'(r) = E^{r). 



(15.7) 



In one dimension this eq uation is exactly solved to give 



in one dmiension tnis eq i 
l|Halperin and La^lligGGj) 



2 

Af/2 



sechKX, 



8k, 



(15.8) 
(15.9) 



with E — — K^/2. Therefore the "optimal action" 
S{Uopt) ^ kVC^o ~ as expected. 

In higher dimensions the corresponding equation is 
not solvable, however, one can extract the energy de- 
pendence of the action by assuming a spherically sym- 
metric optimal fluctuation and an exponentially decay- 
ing at large distances bou nd state to find t he Lifshitz tail 
N(E ) oc exp(-|i;|2-''/2) llLifshitzL Il964bl: iLifshitz et all 
Il988|) . To obtain the pre-exponential factor one needs to 
consider all the wave functions in the potential, and the 
correspondin g analysis has only been carried out in low 
dimensions llHalperin and La: 



nl y been 
i|l96i). 



C. S-wave superconductors 

1. Magnetic and non-magnetic disorder 

Since the effect of the tails is most dramatic for fully 
gapped superconductors, most studies focused on conven- 
tional, s-wave superconductors with pairbreaking mag- 
netic impurities. The general route followed in all the in- 
vestigations is similar to the approach described above: 
given the probability density of different impurity con- 
figurations, and the hamiltonian of the system with the 
potential of each impurity distribution, we find the most 
probable configuration of impurities that gives rise to a 
state at a given energy within the gap. Technical im- 
plementations of this algorithm vary depending on the 
specifics of the problem at hand, see below. 

There are important differences between the physics of 
the optimal fluctuation in such a superconductor and an 
optimal potential well for quasiparticles below the band 
gap discussed in the previous section. First, since the 
superconducting quasiparticles consist of electron pairs 
close to the Fermi surface, their kinetic energy is not 
simply that of a band particle, but is given instead by 
the Hamiltonian 



H = ^T3 + A(r)ria2, 



(15.10) 



where and ct; are the Pauli matrices in the particle- 
hole and the spin space respectively, so that Tiaj is a 4x4 
direct product. Therefore, while the envelope of the tail 
state wave function still varies smoothly over the length 
scale of inhomogeneities in the impurity distribution, the 
rapid oscillations on the atomic scale associated with the 
Fermi surface have to be taken into account. As will be 
seen below, these considerations substantially modify the 
behavior of the tail states. 

Second, the scattering potential is a matrix in particle- 
hole and spin space. In general, an impurity site acts 
both as a potential and a magnetic scatterer, so that the 
total scattering potential is 



U(r) = E 



C/oT3<5(r - r,) + J(r - r,)S, • a] ,(15.11) 



using the Nambu notations. The potential part of the 
scattering, Uq, is not pairbreaking in accordance with An- 
derson's theorem. However, since the size of the optimal 
fluctuation is large compared to the correlation length, it 
is necessary to distinguish between the cases where the 
motion of quasiparticles within the optimal fluctuation is 
diffusive (strong potential scattering, Ar <C 1, r <C Ts, 
where r is the transport lifetime) and ballistic (weak po- 
tential scattering, r ^ t^). Moreover, we should also 
distinguish between strong and weak magnetic scattering: 
if the magnetic scattering is strong there are resonance 
(Shiba-Rusinov) states in the gap, and the tails stretch 
not from the mean-fleld gap edge, but from the localized 
impurity band. If the magnetic scattering can be treated 
in the self-consistent Born approximation, the tail states 
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emerge below the Abrikosov-Gorkov renormalized single 
particle spectral gap, Ag = A(l — (Ar^)"^/^)^/^, where A 
is the self-consistent value of the superconducting order 
parameter. In the AG limit the probability density for 
the magnetic impurity potential is gaussian, as it is aver- 
aged over a large number of impurity sites. In contrast, 
in the unitarity scattering limit there are subgap states 
localized on one or a few impurities; consequently, we 
deal with the Poisson density distribution. These vari- 
ous possibilities provide for a rich variety of behavior that 
is still a subject of active interest. 

All models of tail states due to magnetic impurities 
studied so far ignore interacti ons between the impurity 
spins : it was shown in Ref. ijLarkin and Ovchinnikovl 
Il972|) that the RKKY interaction and glassy behavior of 
impurity spins modify the AG results very weakly. The 
models also treat impurity spins as classical Heisenberg 
spins, and therefore cannot account for the Kondo effect. 
This is justified either when the Kondo temperature of 
individual impurity sites is much smaller than the super- 
conducting transition temperature, Tk (and deple- 
tion of states at the Fermi level prevents screening of the 
local moment), or in the opposite limit, Tk ^ T^, when 
the mom ents are quenched already in the normal state 
iMlillcr- Hartmann and ZittartzL Il97l|) . 

To our knowledge, the first paper discussing the influ- 
ence of non-uniform impurity distribution on the tran- 
sition temperat ure in analogy with Lifshit z's work ap- 
peared in 1968 f Kulik and Itskovichl [1968). These au- 
thors found that, in the limit of average impuritity con- 
centration n <^ Her of the Abrikosov-Gorkov theory, 
there are localized regions that become superconduct- 
ing at a temperature > Tc{n), where Tc{n) is the 
corresponding AG transition temperature. The differ- 
ence between the two was evaluated for parabolic one- 
dimensional variation of th e effective impurity potential. 
iKulik and Itskovichl Il968l also noted that their results 
will be modified if there is non-magnetic as well as mag- 
netic scattering, but did not address this question fur- 
ther. 



2. Diffusive limit, weak magnetic scattering 

If the scattering on individual magnetic impurities is 
weak, the optimal fluctuation is created by large droplets 
of these scattering centers. Since the impurities are un- 
correlated, the probability density for the impurity po- 
tential is Gaussian, which greatly simplifies the analysis. 

Historically, most of the studies have been carried 
out in the diffusive limit. One of the first papers in- 
vestigated the smearing of the gap edge due to local 
fluctu ations in the effective interact ion between elec- 
trons ijLarkin and Ovchinnikovl Il972|) . If the correla- 
tion length of the inhomogeneities, Tc S> where ^ is 
the coherence length of the dirty superconductor, ^ ^ 
(ZJ/A)^/^, and D is the diffusion constant, the order pa- 
rameter simply locally adjusts to the value of the inter- 



action and the density of state is determined by the local 
gap amplitude, 

/>oo 

N{E)= N{E,A)W{A)dA, (15.12) 
Jo 

where W{A) is the probability density of the gap. 

In the opposite limit of short-range correlations in the 
pairing interaction, the finite density of states below the 
mean field gap edge is due to the states spatially localized 
in correlated droplets of size tq ^ C[(^o ~ E)/A]~^^^ 
(increasing rapidly as E —> Aq as in a semiconductor), 
which leads to N{E) oc exp(-[(Ao - E)/A]^^'^ in d = 3. 
As in semiconductors, the high entropy cost of a large 
droplet is offset by the lowering of the kinetic energy of 
the bound state. Indeed, in a clean system with At^ 3> 1, 
and therefore Sq ~ A, we fi nd the characteristic kinetic 
energy, D/r^ ~ V^o - 

Recently it was argued that the above result is flawed 
since it does not account properly for the rapid os- 
cillations of the wave function of t he bound state on 
the s cale of the Fermi wavelength ijMever and Simon^ 
1200 ll) . These authors used a field-theoretical approach 
that maps the disordered superconducti ng system onto 
a no n-linear tr-model (for a review, see lAltland et al\ . 
|2000() to show that, while the droplet size for the optimal 
fluctuation is identical to that obtained by Larkin and 
Ovchinnikov, the subgap density of states is N{E) oc 
exp{ — [(Aq — i?) / A] which gives the exponent 

3/4, rather than 5/4, for d = 3. 

The paper that brought the investigation of the sub- 
gap states in superconductors into the limelight after a 
quarter-century-long hiatus was the study of the den- 
sity of states due to regions where the impurity concen- 
tration is sufficie nt to lo cally d estroy superconductivity 
ijBalatskv and Trugmanl Il997(l . In that case the spec- 
trum of the fluctuation region is similar to that of a dis- 
ordered metallic grain of the same size, L, and depends 
on the mean level spacing of the grain, Sl- The aver- 
age density of states was obtained in two steps. First, 
an average over all realizations of disorder for grains of 
size L yielded Nl{E) ~ Sj^^. Second, the probability 
of finding a fiuctuation region of size L with the criti- 
cal concentration of impurities, ric, for a given average 
impurity concentration, n, Pl^Ucju) was used to de- 
fine the average density of states over the entire sample, 
N{E) ~ J dVPL{nc;n)NLiE). This integral was esti- 
mated to find 

N{E) - exp[-L^{na ln(nc/n) - + n)], (15.13) 

as £' — > 0. Here Lq = (^oO^^^ is of the order of the 
coherence length in a dirty superconductor with / <C Co- 
At energies closer to the gap edge, in the spirit of op- 
timal fiuctuation, it is not necessary to destroy super- 
conductivity completely to generate the tail states. Us- 
ing the instanton approach for the nonlinear a-model, 
Lamacraft and Simons demonstrated how these states 
arise out of inhomogeneous instanton configurations for 
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the action ijLamacraft and Simon 3. 1200(1 l20ni|) . The re- 
sulting optimal action reads 

/A _ P\ 

So = a,(Aor,)2/3(l_Ao^,)-V3)-(2+.)/8 j^^^o_fc^j 

(15.14) 

density of states, varies as N{E) ~ 
exp[-47rg(e/L)'*-2So] exp{-[(Ao - E)/A](^~^y^}. 
Here f/ is the bare conductance and ^ 1 

That approach appears sufficiently general to anal- 
yse nucleation of domains in a variety of sys- 
tems. It was used to re-derive within this frame- 
work ijLamacraft and Simonsl l200l[l the universal gap 
fluctuations in small meta llic grains, first obt ained us- 
ing random-matrix theory dVavil ov et qiJ . l200lt) . namely 
N{E) ~ exp[-(Ao - Ef^^], valid for Aq - E <s: Aq. In 
this regime the spatial extent of the optimal fluctuation 
is greater than the size of the grain, so that effectively 
we are in zero dimensions, d — 0, and the exponent 3/2 
agrees with the general result of Lamacraft and Simons, 
(6 — ci)/4. In the same zero-dimensional limit, but at 
lower energies, E <C Aq, the random matrix theory gives 
N{E) ~ (|£;|/(5'^/^AV^)expf-7rr,(A n - EY/S], where S 
is the mean level spacing in the grain l)Beloborodov et o,l\ . 
I2OOO}. 

3. Diffusive limit, strong scattering 

Recently the field theoretical treatment has 
been extended to t h e cas e of strong scatterers 
llMarchetti and Simonsl I2002D . In that case the 
probability distribution of scattering strength is Pois- 
sonian rather than Gaussian. In the field theoretical 
language this implies that the action cannot be expanded 
to second order in the magnetic potential, as it was for 
the weak potential. Marchetti and Simons circumvented 
this difficulty by considering the dominant contribution 
of droplets densely populated by magnetic impurities, so 
that S, ^ Is I- As we saw above, an impurity band 
emerges within superconducting gap in the limit of near- 
unitary scattering already at the level of the mean field 
theory. Consequently, the tail states extend from the 
edge of the continuum above Ag as well as from the top 
and bottom of the impurity band, see Fig. 1411 According 
to Marchetti and Simons in all these cases the density 
of states varies as N{E) cx exp[-(|£; - E.l/A)^^-'^^^, 
where Ei is the appropriate band edge. Therefore the 
exponent of the action is identical to that found in other 
systems in the diffusive limit. 

4. Ballistic limit, weak scattering 

It was noticed early on that in some systems the mag- 
netic scattering is dominant: upon increasing the con- 
centration of impurities the increase in residual resistivity 
ratio correlates with the sup pression of the su perconduct- 
ing transition temperature llEdelsteinL 1 1 967(l . Since both 



magnetic and nonmagnetic scattering contribute to the 
resistivity, but only the magnetic part suppresses Tc, this 
is an indication of almost pu rely spin-dependent scatter- 
ing. Shytov and co-workers ('Shvtov et all \200^ consid- 
ered the subgap states in this limit in clean (I ^ ^0 or 
Ats 3> 1) limit, when the spectral gap obtained in the 
self-consistent Born approximation nearly coincides with 
the order parameter, Aq ~ A. 

Once again, since the impurities are assumed to be 
weak, the optimal fluctuation for states not too far from 
the gap edge is large and shallow, and the spin-dependent 
potential has the Gaussian probability density. When 
the size of the optimal fluctuation is much greater than 
the superconducting coherence length, Z 3> i S> ^0, the 
motion of the quasiparticles in this potential is ballistic. 
As a result, direct mapping on the non-linear u- model is 
not feasible, and the problem requires quantum mechan- 
ical treatment akin to that in a semiconductor discussed 
above. 

As in that case, we first consider the one-dimensional 
problem. An important assumption (discussed below) 
is that a ferromagnetic fiuctuation maximizes the effect 
of the impurity potential. Choosing the direction of the 
impurity spins along the y axis, performing rotation 172 
CT3 , we remove the vector character of the slowly varying 
potential U, and consider the hamiltonian 

H± = ^^3 ± Aon ± C/(r). (15.15) 

The hamiltonian, however, still remains a matrix in the 
particle-hole space, and the wave functions of the optimal 
fluctuation are the Nambu spinors 

Let us again discuss the physical behavior of the op- 
timal fluctuation qualitatively. We linearize the kinetic 
energy near the Fermi surface, so that typical kinetic en- 
ergy in an OF of size L is ^ ~ vp/L. Then the energy of a 
quasiparticle in the optimal fluctuation (measured from 
the Fermi energy) is E ~U+ a/Aq -f Vp/L^. For the en- 
ergies close to the superconducting gap, (Ao — E)/Ao <C 
1, the OF is large (L ^ = vp/Ao) and shallow 
{\U\/Ao < 1), so that E- Ao U + vj,/{AoL^). In- 
troducing the dimensionless energy e = E/ Aq, we ob- 
tain, in analogy with the arguments above, |J7|/Ao — 
^g/L^ ~ 1 — e. Notice that the size of the fluctua- 
tion is L ~ Co/Vl — e 3> Co- As a result, we find 
S[U] « LU'^/Ul = AI^q{1 - efl'^lUl. From the defi- 
nition of C/q it follows that 

- In « S\Voyt\ ^ (AoT,)(l - efl\ (15.16) 

iVo 

The energy dependence in Eq. H15.16|l is identical to the 
result of Lifshits in rf = 1, despite the linear, rather than 
quadratic, dependence of the kinetic energy on the size 
of the droplet. This follows from the smallness of this 
energy compared to the gap: even though ^ oc 1/L, the 
expansion is in ^ . 

The minimization of the saddle-point action proceeds 
exactly following the steps in section IXV.BI For spin 
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"up" particles £+[U] 
respect to U gives 



Minimization with 



U{x) = -\U^{^\ 



5H+ 
~SU 



(15.17) 



In principle this variational derivative includes the effect 
of the self-consistent suppression o f the gap. However , 
it can be explicitly demonstrated (|Shvtov et all 1200,'^ 
that the effect of self-consistency is small. Then, in 
exact analogy to the semiconductor problem, U{x) — 
-XUq{^*{x)^>{x)), where {'9*^) denotes the scalar 
product in particle-hole space. In turn, Schrodinger 
equation takes the form 



d 

-j-fF^Ta + Aon - AC/o^(«'**) 



(15.18) 



This equation is solved by introducing the bilinear forms 
'if* {x)Ti'i> (x) . These bilinears play the role of the 
Halperin-Lax wave function in the Nambu space, and 
yield 



Ro 



1 



1 



,(15.19) 



Co arccos e e + cosh(2a;Vl - e^/Co) 
Ri = i?o(e + Co ^0 arccos e), (15.20) 

R2 - 



(15.21) 



and i?3 = ijShvtov et all l2004l). The phys ical potential 
of the optimal fluctuation is ijShvtov et all 12003) 



Uix) 



1 



2 Ao e + cosh(2a;\/r^ /Co) ' 
which corresponds to the value of the action 



S[U] ~ 87r(Aors) \/l — e'^ — e arccos e 



(15.22) 



(15.23) 



For e w 1 the length scale of the optimal fluctuation is 
Co/vl — ?, its depth is U ^ Ao(l — e^), and the density 
of states N{E) ^ exp[— (1 — e^)'^/^], in complete agree- 
ment with qualitative estimates. 

The most important observation of lShvtov et all l2003l 
is that in higher dimensions the optimal fluctuation is 
strongly anisotropic, in contrast to both the conven- 
tional semiconductors and superconductors in the diffu- 
sive limit. This is a direct consequence of the composite 
nature of superconducting quasiparticles: they are made 
out of objects that move with the Fermi velocity. The 
wave function of the subgap state is concentrated along 
the quasiclassical trajectory, which is a chord in a po- 
tential of any shape. Consequently, there is little energy 
cost in reducing the size of the OF in the "transverse" 
direction, while the smaller volume makes such fluctua- 
tions more probable, see Fig. 021 As a result, the opti- 
mal fluctuation is strongly elongated in one (x) direction. 
The wave function of the bound state can be written as 



\E'(a;,y) = exp(?fci?a;)$(a;, y), where y denotes the trans- 
verse d—1 coordinates, and $ is a slowly varying function. 
The kinetic energy of the quasiparticle is 



C*^ 



Akpx 




(15.24) 

The transverse size of the fluctuation can therefore be 
reduced until the second term becomes comparable to 
the first, i.e. Ly ~ [XpL^Y^^ , where Xp ~ /cj^^ is the 
Fermi wavelength. Consequently, |C/|/Ao ~ 1 — e and 



Co/ VI - e, and 



(1-6)^ 



(15.25) 

where Ep is the Fermi energy. Consequently, the density 
of states, N{E) ~ exp[-(l - e)^^-'^)/^]. 

The action for this anisotropic fluctuation is smaller 
than that for an isotropic droplet with the same energy 
of the bound state, by a factor of {Ep / - 
g)-(£'-i)/4^ SO that the corresponding DOS is exponen- 
tially higher. 

Since the optimal fluctuation is a result of a sad- 
dle point approximation for the functional integral, 
Eq. H15.1|l . it is only vahd when S\Uopt\ ^ 1, or 



2(d-l) 



(15.26) 



For d = \ this condition becomes 1 — e ^ (AqTs) 
while for d = 3 it does not depend on the gap amplitude, 
l-e> {kpl)-^. 

It is possible to compare the densities of states given 
by different approaches at the cross over scale be tween the 
diffusive and the ballistic regimes l|Vekhter et at. 2003). 
A transition to the diffusive regime occurs when the size 
of O F L > vpT.,, or 1 - e < (AoTs)-^. The result of 
Ref. l|Lamacraft and Simon^ 2000) for AqTs >• 1 is S/j = 
(AoT,)^/3(£;^/Ao)'^-i(l-e)(^-'')/4. C onsequ ently, at the 
crossover point the action from Eq. (|15.25|l is smaller, 
Sd/So 0^ iEp/Ao)^'^-'^^HAoT,y^'^ > 1, and the OF 
found by Shytov et al. corresponds to a greater DOS. 
Therefore the structure of the OF near the crossover be- 
tween the ballistic and diffusive regimes still resembles 
closely that given above. As the size of the OF increases 
even further, the anisotropic fluctuation becomes insup- 
portable due to diffusive motion. 

Balatsky and Trugman l)Balatskv and TrugmanLfl997j) 
considered only the DOS at -E = due to the suppres- 
sion of superconductivity by paramagnetic impurity po- 
tential. They needed a large volume ffuctuation, > C^j 
which is less probable and yields l ower DOS than that of 
Eq. 115.25(1 . IVekhter et dl l2003l checked whether local 
suppression of the gap from Aq to E due to a large num- 
ber of impurities with uncorrelated spins (as opposed to a 
ferromagnetic OF above) is advantageous. For 1 — e ^ 1 
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the local pairbreaking rate, 7, needed to reduce the gap 
to E is « 1 + (1 — e)(AoTs)^/^, and the volume of the 
region has to be at least equal to that of the anisotropic 
OF to avoid high kinetic energy cost (this is an underesti- 
mate since it ignores proximity coupling to bulk) . In that 
case the optimal action Sbt/Sq ~ (AoTs)^/'^(i?_F/Ao)c, 
where c = Uimp^p is the atomic concentration of impu- 
rity atoms. As a result, for realistic values of c and clean 
samples Sbt 3> So, and the DOS given by the action in 
Eq. 115.251) is higher. 

Therefore the ballistic limit of the action obtained 
by IShvtov et all 120031 is expected to be valid near the 
crossover between the ballistic and the diffusive regimes. 



5. Ballistic regime, strong scattering 

As of today, we are not aware of any investigations 
of the structure of the optimal fluctuation in the bal- 
listic regime, when there exist bound states on individ- 
ual magnetic impurities. It is reasonable to assume that 
the result differs from the standard Lifshitz formula for 
the same reason as in the section above: the wave func- 
tions of the states localized on magnetic impurities in 
superconductors oscillate with the Fermi wavelength, see 
Sec. lVIl As a result, in the dilute impurity limit, the shift 
of the energy level localized on, for example, two impuri- 
ties located at distance R » will be supp ressed by 
the typical factor exp(— i?/^o) teu sinovL Il968|) . Conse- 
quently, the states significantly below the impurity band 
must be created by a large number of impurities or impu- 
rities located on neighboring lattice sites. This problem 
still awaits further investigation. 



6. Summary 

In s-wave superconductors with magnetic impurities 
the density of states does not vanish irrespective of the 
concentration and nature of the impurity scattering. The 
tails of the density of states extend into the mean field 
gap. Therefore all superconductors with magnetic impu- 
rities are gapless. This behavior is qualitatively illus- 
trated in Fig. 



XVI. SUMMARY AND OUTLOOK 

While considering the role of impurities in conventional 
and unconventional superconductors, this review focused 
on theoretical and experimental results that highlight the 
new physics beyond standard Abrikosov-Gor'kov theory, 
Anderson theorem and average lifetime effects. The stud- 
ies of disorder in s-wave superconductors were carried out 
in detail in the 1960's. We discussed more recent results 
in this field. Our main emphasis has been on how indi- 
vidual impurities influence local electronic states in their 
immediate vicinity, and on deviations from the standard 



Abrikosov-Gorkov theory on mesoscopic scales. This fo- 
cus is dictated both by the advances in experimental tech- 
niques, which can now use NMR methods and STS mea- 
surements to probe electronic states with atomic spatial 
resolution, at the scales where impu rities perturb their 
surrounding, ijFischer and et a 11. 120041) . and the concomi- 
tant development of new theoretical approaches. 

The stimulus for such extensive studies is that impuri- 
ties are markers that allow to reveal the nature of correla- 
tions and pairing of the state where impurities are placed. 
Indeed particular pattern of impurity-induced electronic 
states is closely connected to the symmetry of the super- 
conducting gap and helps us to understand the nature of 
superconducting pairing. If strong electronic correlations 
in the ground state are present, they also are reflected in 
details of impurity induced states. Therefore watching 
the waves created by throwing a pebble in the pond of 
correlated electrons helps us understand the properties 
of the underlying electronic liquid. 

We kept the discussion general to allow applications 
to other systems and materials. For instance, this was 
our rationale for employing the BCS state to describe 
superconductivity. We believe that it is a good approx- 
imation in heavy fermion, organic superconductors and 
SrRu04, at very low energy. At the same time, devia- 
tions from this mean field picture may provide additional 
details on the underlying physics of the particular mate- 
rial. Majority of the data at a moment are obtained in 
high- Tc materials. It is clear that similar local effects 
are present around impurities in other unconventional 
supe rconductors, e.g. in Na^ Co02 • 2/H2O superconduc- 
tors l|Wang and Wand. 12004') . although we are not aware 
of any data on single impurity states in these materials. 
Given the importance of the impurity states, this field 
will undoubtedly be extended to other systems by future 
experiments. 

Outlook for the future. New ideas and directions con- 
tinue to emerge in the studies of electronic properties 
induced by impurities. The suite of new experimen- 
tal tools that address local electronic effects, such as 
STM, will help to clarify the role of interference be- 
tween several impurities, and pave the way towards con- 
necting the microscopic local states with average prop- 
erties. Recent t heoretical work addressed some aspects 
of this subject llAndersenL l2003t I Atkinson et alL [ 200: " 



Morr a nd Balatskv. "2003': 'Morr and Stavropoulosl 12003 : 
Zhu et al 2003. 2004b.} . and is awaiting direct compari- 
son with experiment. 

Another promising avenue is combining the spatial 
resolution of STM-STS with the time resolution. The 
subject is still in its infancy, both theoretically and ex- 
perimentally, but hold immense promise for the future. 
Sec. IXII reviewed some of very recent work in this direc- 
tion. Temporal and spatial characterization of the states 
generated by dynamical impurities allow exploration of 
the correlations inside the electronic state in which im- 
purity is placed. One obvious example where such char- 
acterization is crucial is the Kondo effect in a supercon- 
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ducting state. It is desirable to have a time resolved 
measurements that allow to visualize the Kondo effect 
in a superconductor. Another interesting problem that 
needs further elaboration is a role of collective modes in 
impurity-induced states. We are only starting to investi- 
gate these questions, as discussed in Sec. IXIII 

Real progress on these problems will be made when we 
have real data. As usual, one should expect that the data 
will have surprises that were not anticipated in simple 
theoretical models. This will motivate further theoretical 
studies, stimulate more measurements, and therefore will 
lead to a rapid further development of the field. They 
can provide space (and time) resolved window into the 
intimate workings of the correlated electron matter. We 
have every reason to be enthusiastic and optimistic about 
the future the field of impurity states in superconductors, 
and in other correlated electron systems. 
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List of Symbols 



Quantity 


Explanation 


a 


Lattice parameter 


b[b<) 


Bosonic annihilation (creation) operators 




Fermionic alihihilation (creation) operators 


a 


Spatial dimension 


D 


Half energy bandwidth 


Ao 


Superconducting energy gap 




Momentum-dependent superconducting energy gap 


(t>n{r) 


Electron eigenfunction 


Ep 


Electron Fermi energy 



G(t, r'), G{t, r) Electron Green's function in coordinate space 

G(a;„,k), G(k, a;„) Electron Green's function in Matsubara frequency and momentum space 

H, H, Hint Hamiltonian 

J, Jo , Jc Exchange coupling 

L Linear dimension of a system 

H Chemical potential 

N{e) Electron density of states 

N{e, r), N{E, i) Electron local density of states 

tp{r{ip^ (r)) Fermionic field operators in continuum space 

|\t'o) BCS variational wavefunction 

l^^-i) Excited variational wavefunction with single particle excitation present 

r Spatial coordinates 

<T Pauli matrices in spin space 

T Pauli matrices in Nambu space 

Va0-yS, Va0 Superconducting pairing interaction 

S Local spin operator 

t, t' Electron hopping integral 

u Electron-like Bogoliubov quasiparticle wavefunction amplititude 

T{(j) T-matrix 

T Temperature 

V Hole-like Bogoliubov quasiparticle wavefunction amplititude 

U Hubbard on-site electron-electron interaction 

Uq Impurity scattering potential 

W Half energy bandwidth 

Wk _D-density-wavc order parameter 

^0 BCS superconducting coherence length at low tc^niperaturc^s 

^(T) BCS temperature dependent coherence k ugtli 

I 
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Figures 



FIG. 1 Multiple scattering on a single impurity. Thick (thin) 
line denotes full (bare) Green's function, and the dashed line 
denotes scattering process. The second line defines the T- 
matrbc according to Eq. 13.191 . 



FIG. 2 Impurity bound state in a metal at energy uq is formed 
as a result of a multiple scattering. 

FIG. 3 Graphic solution of the Eq. (17.511 for large Uo is 
shown. Only physically relevant intersections with small 
imaginary part f2o" <C are shown. In the region where 
the imaginary part of the local Green's function exceeds the 
real part, the resonance is broadened and merges with con- 
tinuum. Resonances below (or above, for a different sign of 
Uo) the fermionic band are the sharpest, with most of spec- 
tral weight, and the virtual bound state inside the gap is well 
resolved for large Uo (small c). 

FIG. 4 Illustration of the cross shaped nature of the impurity 
state. Shown is the spectral density Aa{v, ±r2o) as a function 
of position and spin in units of A'qAo for a) /i = and b) 
^ — —W, 2W is the bandwidth, in a two-dimensional d-wave 
superconductor as a function of position around a classical 
magnetic moment ( No Jo = 10 and Uo — 0) located at r = 
0; a is the lattice spacing. These results are computed self- 
consistently with ^ = 10a. At half filling, the spectral density 
obeys particle-hole symmetry: ^| (r, ilo) = ^j. (r, — f2o). The 
energies of the shown virtu al-bound states are a ) Qo = 0.05Ao 
and b) fio = 0.5Ao. From jSalkola et a^.lll997^ 

FIG. 5 Self-consistently determined gap function near scalar 
impurity in a 2D d-wave superconductor. Gap suppression is 
strongly localized nea r impurity sit e aside from weak oscillat- 
ing tails. From ( Franz et all Il996ri 

FIG. 6 The LDOS as a function of energy at the impurity site 
(left panels) and at one of its nearest neighbors (right panels) 
in a 2D lattice. The upper panels are for various values of 
repulsive potential, Uo ~ (black line) , 2 (blue line) , 5 (green 
line), and 10 (red line) while the lower panels are for various 
values of attractive potential, Uo = (black line), -2 (blue 
line), -5 (green line), and -10 (red line). The band structure 
parameter values are t = 1, t' = 0, and the chemical potential 
M = 0. 

FIG. 7 Same as Fig.|S|except that the band structure param- 
eter values are t = 1, t' — —0.2, and the chemical potential 
/i = -0.784. 
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FIG. 8 Same as Fig. El except that the band structure pa- 
rameter values t — 1, t' = —0.3 and the chemical potential 
/i = -1.0. 



FIG. 9 An impurity state in a high Tc superconductor: (a) 
The DOS in the pseudogap regime used in this article (see 
also [11]) and (b) the DOS in the superconducting state as 
was used in [1] . In both phases there is a resonant state. 



FIG. 10 (a) The density of states N{uj) = -go{uj)/TT. Around 
the pseudogap states are only partly depleted e.g. No{u>) — 
No\uj\/ApG, where No{lli) = A^o for Apg < \(^\ < W/2 with 
W the bandwidth, (b) The real part g'o{u}) of Green's function 
together with 1/Uo and Uo positive, il' is the real part of the 
solution of the equation (?o(f^) = 1/C^o close to zero and there- 
fore with sharp bandwidth, (c) The impurity induced reso- 
nance at Q' = — Apg/2(7oAo ln(2L'^oAo). Because the other 
three solutions of Eq. (18.111 have much broader bandwidth, 
they are not depicted here. All the figures are taken on the 
impu rity site. From Kruis, Martin and Balatsky iKruis et all 
1200111 



FIG. 11 DDW-DOS for the clean case (solid line) and in the 
presence of a non-magnetic impurity with Uo = 1 eV: (1) DOS 
on the impurity site, (2) DOS on the nearest-neighbor site, 
and (3) DOS on the next-nearest-neighbor site. The other 
parameter values are: t=_30<yrneV , Wo = 25meV, t' = 0, 
and ^ = 0. From Morr llMoril |2002|) . 



FIG. 12 (a) Fermi surface in the DDW state with t' = 
— 0.3t, ^ — — 0.91t (corresponding to a hole-doping of 10%) 
and Wo = 25 meV. The hole pockets are centered around 
(±7r/2,±7r/2). (b) DOS in the DDW state with the same 
band parameters as in (a), for the clean case (solid line) 
and in the presence of a non-magnetic impurity with Uo = 1 
eV: (1) DOS on the impurity site, (2) DOS on the nearest- 
neighbor site, and (3) DOS on the next-nearest-neighbor site. 
Inset: SC DOS for the same band parameters as in (a). From 
Morr llMoril. l2002fl 



FIG. 13 Local density of states with Ao = 0.68t, = —0.3t 
and Uo — lOOt. (a) A(r„„,Lj) versus co. Solid lines: n„ = 
10"'',10"*,10"^, and 5 x 10"^ with decreasing peaks. The 
dotted line is the LDOS at = and Uo = for comparison, 
(b) N{r, 0.05t) at n„ — 0. The impurity is at the center, (c) 
The same as (b) for — 5 x 10~^ . The gray scale is the 
same i n (b) and (c ). The other parameter t' — 0. From 
Wang llWanell2002t ). 



FIG. 14 The local effect of a magnetic moment on the low- 
energy spectral density in an s-wave superconductor. 



FIG. 15 Two variational states are shown schematically. I'I'o) 
is a standard BCS wave function that contains only paired 
particles and has unscreened impurity spin 5*. I^*!) is a varia- 
tional wave function that describes the formation of the bound 
state between particle with the spin opposite to the local spin 
( for antiferromagnetic coupling); this state is inherently a 
non BCS state and electonic spin quantum number differs by 
one unpaired spin compared to I'J'o). 



FIG. 16 Energies of two variational states are shown. |^'o) is 
a standard BCS state with energy Eo- is a variational 

state that describes the formation of the bound state between 
particle with the spin opposite to the local spin with energy 
El. Level crossing between states with different symmetry 
occurs at some critical value of the coupling Jcrit- This is 
an example of a first order quantum phase transition with no 
divergent length or time scale associated with it. 



FIG. 17 a) The bound-state energy Qo, b) the spectral 
weight of the pole for positive and negative energies in 
units of NoJo,No = Np, c) the spin polarization (s(r = 0)), 
and d) the gap function A(r = 0)/Ao at the impurity 
site r = as a function of Jo in the s-wave superconduc- 
tor. Lines denote the T matrix results for the uniform or- 
der parameter and symbols the self-consistent mean-field re- 
sults on a square lattice at half filling. The quantities of 
the impurity-induced intragap quasiparticle state belonging 
to the branch Jo < Jcrit are denoted by solid lines and 
solid symbols, whereas those ones belonging to the branch 
Jo > Jcrit are marked by dash ed lines and open symbols. 
Taken from llSalkola et a^.Ul997^ . 



FIG. 18 Cartoon of the intrinsic vr junction near magnetic 
impurity in s-wave superconductor. 



FIG. 19 Calculated tunneling density of states for the four- 
site Kondo impurity model at 15% hole doping with a realistic 
band structure {t = 0.15 eV, t' = -t/4, t" = t/12), Ao = 0.04 
eV, and m = -0.14 eV. The Kondo coupling is J = 0.09 eV, 
the potential scattering U = 0. Top: Local DOS vs. energy 
for the impurity site (red) and the nearest (blue) and second 
(green) neighbor sites. Bottom: Spatial dependence of the 
local DOS at w = -2 meV. Left: Local DOS in the Cu02 
plane. Right: Local DOS after applyin g the filter effect pro - 
posed by Martin, Balat sky, and Zaanen IIMartin et all\200^ . 
From Vojta and Bulla llVoita and Bullall200ll) 



FIG. 20 Same as Fig. 1191 but with potential scattering 
U = t = 0.15 eV. Here, J = 0.065 eV. The lower panel 
shows the local DOS at uj — +2 meV. From Vojta and 
Bulla JVoita and BullaLl200l|) . 



FIG. 21 Same as Fig. 1191 but with potential scattering 
U = 4t = 0.6 eV. Here, J = 0.04 eV. The lower panel 
shows the local DOS at ui = +3 meV. From Vojta and 
Bulla JVoita and Bullall200lD . 
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FIG. 22 The solution for the DOS (black line) and its en- 
ergy derivative (red line) are presented for a local boson 
mode scattering in a d-wave superconductor. The normal 
self-energy was treated self-consistently as a full solution of 
the Eq. 111.3L ignoring vertex corrections and gap modifica- 
tions. Apart from the feature ai uj — ujq we get also strong 
satellite peaks at A-|-aJo that are a consequence of a coherence 
peak in DOS of a d-wave superconductor. These satellites are 
a specific property of a superconducting state and will not be 
present in a pseudogap state. These features are best seen in 
Energy scale is given in units of A, the dimensionless 
coupling constant is taken to be 1. For comparison we plot 
the results for local mode frequency ujo/^ = 0.2,0.4,0.6 in 
the first three panels. The lower panel gives the results for 
the asymptotic analytic solution, that assumes <^ A using 
Eq. Ull.l^ . for LOO = 0.4. The overall features are similar for 
both cases, however the an alytic solution shows a somewhat 
larger feature. Taken from llBalatskv et aZ.Ll2003t) 



FIG. 23 Appearance of the satellite peaks for an impurity 
resonance uiimp at Uimp ± loq is shown schematically. The 
satellites will have different spectral weight. Imagine we inject 
into system an electron at energy coimp+i^o- To create a peak 
at one need to excite local mode and the energy of the 

electron will be equal to the difference between local state 
and local mode energies. Similarly, to obtain the peak at 
ujimp from injected electron at energy ujimp — i-jq one needs 
to add local mode energy to an electron. For this process 
to occur the local mode has to be excited to begin with and 
hence this process will have very low weight at low T. These 
two processes will also have different matrix elements. Overall 
relative weight of the side peaks is proportional to J^Nq which 
we assumed to be small. In case of magnetic scattering, when 
LJo = gu sB the splitting w ill be tunable by the field. Taken 
from ( Balatskv et all 1200311 



FIG. 24 The density of states for a clean system with g — 
and 0.2. 



FIG. 25 Fourier spectrum at i5 = —E,. for the constant g — 
0.2 and the structure factor Uq ~ Uo- The parameter Uo = 
0.3. 



FIG. 26 Fourier spectrum at iJ = —E,. for the constant g = 
0.2 and the structure factor Uq = Uogo/{qo + 4[cos'^(gi^/2) + 
cos^ (17^/2)]}. The parameter Uo = 0.3 and qo = 0.5. 



FIG. 27 Left panel: The dl/dV spectra measured near (A) 
Mn, (B) Gd, and (C) Ag atoms and far away from the im- 
purity where local density of states can be fit by the BCS 
theory. Right panel: Constant-current topographs and si- 
multaneously acquired dl /dV images show the spatial extent 
of the bound state near Mn and Gd adatoms. (A) Constant- 
current (32 A by 32 A) topograph of a Mn adatom. (B) Image 
of dl /dV near the Mn adatom, acquired simultaneously with 
the topograph in (A) by using an ac detection. The areas 
where dl /dV is reduced (dark) show the extent of the bound 
state. This reversed contrast comes about because a dc bias 
voltage was chosen well above the energy of the bound state, 
where the bound state affects dl/dV only indirectly by con- 
tributing to the total current I. (C) Constant-current (32 A 
by 32 A) topograph of a Gd adatom. (D) Image of dl/dV near 
the Gd adatom, acquir ed simultaneously w ith the topograph 
in (C). From Yazdani llYazdani et aLLIlQQTh . 

FIG. 28 Differential tunneling spectra taken at the Zn-atom 
site (open circles) and a location far away from the impurity 
(filled circles). Note that even on the impurity site one has 
peaks at both positive and negative bias albeit of very differ- 
ent magnitude that are reflection of the particle hole character 
of the impurity resonance. To fit the data one can use a simple 
potential scattering model with essentially unitary scattering 
phase shift 8 = 0.487r. Phase shift is related to a impurity 
potential Uo via simple formula: cot 9 = ^^^^^ . From Pan 
et al. llPan et a^.ll2000b^ . 

FIG. 29 Differential conductance spectra above the Ni atom 
and at several nearby locations. Differential conductance 
spectra obtained at four positions near the Ni atom showing 
the maxima at eV — ±f2i. Intensity as a function of po- 
sition relative to impurity site reverses upon change of the 
bias sign. This effect is explained as a result of particle 
and h ole components of t he impurity state. From Hudson 
et al. llHudson et aZ.ll200l[) . 

FIG. 30 High-spatial-resolution image of the differential 
tunneling conductance at a negative tip voltage bias eV = 
-~l.2mV at a 60 X 60A^ square. Also shown d-wave gap nodes 
orientation and lattice sites to indicat e that impurity sta te is 
registered to lattice. From Pan et al. llPan et aZ.ll2000bD . 

FIG. 31 Tunneling DOS for tunneling on Ni impurity site. 
Note that there are always states at opposite bias as well. 
The peak intensity is largest on either positive or negative 
bias depending on the position. To fit the data one need t o 
use both Uo and J. From Hudson et al. jHudson et ttZ.ll200lD . 



FIG. 32 The particle and hole components of the impurity 
wave function for a magnetic impurity in a s-wave supercon- 
ductor is shown. A) Impurity wave function 'I'B(r) and B) 
r^*I's(r) are shown. The maxima of particle and hole com- 
ponents occur at different positions. This results in the dif- 
ferent image of the impurity state, seen on positive and neg- 
ative bias. This effect is a general property of a supercon- 
ductor regardless of t he symmetry of the pairing state. From 
iYazdani ei a^.lll997^ . 
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FIG. 33 A representative set of seven scattering vectors 
cti(E) of the 'octet' mode l. Reproduced with permission 
from ijMcElrov et alli200^ . 



FIG. 34 The real space image of different orbitals on Cu, 
nearest O and nearest Cu sites are shown. Quantum me- 
chanical interference produces the filter effect that changes 
the d istribution of the impurity state intensity (jMartin et all 
l2002fl . 



FIG. 35 Density of states in the Abrikosov-Gorkov theory 
of magnetic impurities in superco nductors. Here F = t^^ . 
Reproduced with permission from JSkalski et a;.LIl964D . 



FIG. 36 Plot of the dependence of the order parameter, A, 
transition temperature, Tc, and the single particle spectral 
gap, ila here, on the scattering rate F = t^^ . Reproduced 
with permission from JSkalski et all\l964) . 



FIG. 39 Reduced transition temperature normalized to pure 
system as a function of the impurity concentration for differ- 
ent eletcron-phonon coupling, Ap. Th e impurity concentration 
c = n,rap/{2TvfNoTco- From iJarrelk .igQaV 



FIG. 40 The spatial structure of the optimal fluctuation in 
the ballistic and the diffusive limits. 

FIG. 41 Qualitative picture of the density of states in an s- 
wave superconductor with magnetic impurities. Blue shade 
denotes regions where the mean field density of states is finite. 
Red shading signifies the finite DOS induced by the deviations 
of the local impurity distribution from the average. The den- 
sity of states in these tails is exponentially small but finite. 
If the impurities are weak, the well-defined impurity band is 
absent and the tail extend from the mean field gap edge. 



FIG. 37 Evolution of the spectral gaps and density of states 
for strong magnetic impurities (eo <C Ao). Left panel shows 
the available states (shaded) as a function of the impurity con- 
centration. Right panel shows the qualitative features of the 
density of states for different values of the impurity concentra- 
tion, labeled by vertical lines A,B,C,D on the left. Critical 
concentration of impurities corresponds to line B, when the 
impurity band first touches lj = 0. At the same time, spectral 
gap between the top of the impurity band and the bottom of 
the continuum states persists to higher impurity concentra- 
tion (line D). 



FIG. 38 Evolution of the spectral gaps and density of states 
for weak magnetic impurities (eo ^o). Left panel shows 
the available states (shaded) as a function of the impurity 
concentration. Right panel shows the qualitative features of 
the density of states for different values of the impurity con- 
centration, labeled by vertical lines A,B,C on the left. The 
impurity band and the continuum above the gap merge at a 
low impurity concentration,see line B, and further evolution 
of the density of states is very close to the predictions of the 
AG theory. Critical concentration (line C) marks the onset 
of gapless superconductivity. 



Tables 





5-WAVE 


P-WAVE 


D-WAVE 


POTENTIAL SCATTERING 




+ 


+ 


MAGNETIC SCATTERING 


+ 


+ 


+ 



TABLE I Effects of potential and magnetic impurity scat- 
tering on the S-, p- and d-wave superconductors is shown 
qualitatively. indicate that impurity scattering is a pair- 
breaker and "— " is that impurity scattering is not a pair- 
breaker. There is a qualitative difference between the poten- 
tial scattering in s-wave superconductors and any other case. 
Potential scattering impurities are not pairbreakers in s-wave 
case due to Anderson theorem. This is an exceptional case. 
For any other case any impurity scattering will suppress su- 
perconductivity. Obviously the details depend on scattering 
strength and other details. At high enough concentrations 
both magnetic and nonmagnetic impurities will suppress su- 
perconductivity regardless of symmetry. 
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